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This  report  preaetrs  a  sesi-empirical  (non-Gauss ian)  model  for  the 
obscurants  chat  say  be  generated  by  military  activity.  ”hese  in¬ 
clude  smoke  and  obscuring  aerosol  munitions,  dust  clouds  generated 
by  non-nuclear  explosive  munitions,  and  dust  generated  by  military- 
like  vehicles.  The  field  data  and  some  laboratory  experiments  are 
discussed,  and  numerous  validations  are  provided. 
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SUMMARY 

This  report  presents  a  semi -empirical  (non-Gauss ian)  model  for  the 
obscurants  that  may  be  generated  by  military  activity.  These  include 
smoke  and  obscuring  aerosol  munitions,  dust  clouds  generated  by  non-nuclear 
explosive  munitions,  and  dust  generated  by  military-like  vehicles.  The 
field  data  and  some  laboratory  experiments  are  discussed. 

Basically  the  approach  adopted  here  is  to  generate  a  predictive  model 
that  will  permit  the  analyst,  as  well  as  the  field  operator,  to  establish 
the  thermal  and  space-time  history  of  the  obscurant  cloud  and  the  concentra¬ 
tion  of  tne  obscurant  therein.  The  input  data  for  the  predictions  are  readi¬ 
ly  available  to  the  user.  In  a  similar  vein,  the  explosive  dust  model  sub¬ 
sumes  the  smoke  dispersion  model. 

The  basis,  justification,  and  validation  for  the  various  algorithms 
are  presented,  and  the  mathematical  expressions  to  find  the  mean  values  and 
the  deviations  of  Che  obscurant  concentration  spatial-temporal  behavior. 
Directions  for  model  improvement  are  discussed. 

Lastly,  an  attempt  is  made  here  to  present,  at  the  ends  of  sections, 
a  calculac tonal  schema  to  compute  the  necessary  parameters  to  establish 
the  effect  of  the  various  obscurants  on  battlefield  electro-optical 
systems. 
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1.1  BACKGROUND 


This  is  a  final  report  on  the  General  Research  Corporation  (GRC) 
effort  dealing  with  the  formulation  of  a  model  to  describe  obscurants 
and  their  behavior  in  a  military  environment.  The  effort,  during  the 
past  two  years,  has  been  sponsored  by  the  Night  Vision  &  Electrooptics 
Laboratory  under  Contract  DAAK02-74-C-0366 ,  and  ably  monitored  by  L.  Obert. 

The  model  initially  addressed  the  dispersion  of  smokes  and  other 
obscuring  aerosols.^  Emphasis  was  given  to  a  model  that  can  readily  be 
used  in  combat  simulation  models;  however,  the  physical  aspects  of  the 
model  were  sufficiently  broad  and  accurate  that  it  could  be  used  for  a 
variety  of  applications.  Subsequently  the  model  has  been  extended  to  in¬ 
clude  explosive  and  vehicle-generated  dust. 

The  explosive  dust  model  has  two  stages.  The  initial  dust  cloud  is 
determine-  --y  -.he  hydrodynamics  of  the  explosion  (high-explosive)  process, 
and  the  dispersion  of  this  cloud  is  described  by  the  smoke  dispersion  model. 
Hence,  the  overall  utility  of  the  model  becomes  broader.  It  should  be 
noted  that  the  explosively  dispersed  bulk  whice  phosphorus  shell  (155  mm) 
can,  if  desired,  be  modeled  like  explosive  dust. 

During  the  past  several  years  the  US  Army  has  conducted  extensive 
field  test  programs.  Data  derived  from  these  efforts  have  been  employed 
here  to  some  extent  to  assist  in  validating  the  concepts  used.  However, 
full  access  and  use  could  have  furthered  the  formulation  of  a  more  thorough, 
flexible  model.  The  reasons  for  this  statement  are  as  follows: 

1.  The  field  test  data  represent  a  large  body  of  unique 
data  on  the  generation  and  turbulent  dispersion  of 
various  aerosol  generators,  particularly  quasi- 
instantaneous  types. 


^Zirkind,  R.  ,  An  Obscuring  Aerosol  Dispersion  iModel,  GRC  Report  231,  Vols. 
I  and  II,  December  1978. 


2. 


Verification  of  the  laboratory  experiments  and  theoretical 
results  would  have  been  possible. 


Unfortunately,  the  dat3  have  been  "shoehorned"  in  a  preconceived  dispersion 
model.  Hence  new  model  development  by  the  US  Army  is  a  remote  possibility 
and  many  Important  questions  may  remain  unanswered. 

1.2  APPROACH 

The  objective  of  the  sponsored  effort  has  been  the  development  of  a 
model  that  can  readily  predict  the  properties  of  a  dispersing  aerosol  cloud 
produced  by  military  munitions  and  smoke  generators.  To  meet  this  objective 
a  semi -empirical  approach  was  adopted  to  be  consistent  with  conventional 
military  procedures,  field  conditions,  and  physical  phenomena.  Specifically, 
the  chemical  officer  has  a  limited  set  of  data  for  planning  purposes  and 
therefore  the  inputs  should  be  consistent  with  this  fact. 

Of  greater  importance  is  the  fact  that  existing  representations  of 
diffusion  based  upon  Gaussian  distribution  functions  appear  too  inadequate 
for  aerosol  releases  at  or  near  ground  level  where  the  entire  dispersal 
process  occurs  within,  at  most,  a  several  minute  period.  Examination  of 

gust  data  for  l-sec  sampling  and  5-min  and  10-min  averaging  periods  exhibit 

2 

large  dirferenceu.  This  implies  that  a  very  large  sample  of  instantaneous 
concentration  values  would  be  required  to  approach  a  Gaussian  distribution 
and  obtain  the  mean  value.  On  the  ocher  hand,  the  fluctuations  per  se  are 
of  equal  import  since  perils  of  low  concentration  along  the  line  of  sight 
of  a  sensor,  comparable  to  acquisition  time,  would  defeat  the  usage  of  the 
obscurant.  Further  the  several  o’oseuraats  used  are  initially  buoyant 
plumes;  i.e. ,  they  release  significant  thermal  energy,  and  therefore  a 
fraction  of  their  temporal  behavior  is  dominated  by  heat  exchange  and  not 
momentum  exchange. 

For  the  above  reasons,  the  approach  taken  was  to  establish  the  mean 
concentrations  on  the  basis  of  experimental  data  using  Pasquiii  stability 

V "  "  " 

'"Exploring  the  Atmosphere 'a  First  Mile,"  (H.  H.  letcau  and  D.  Davidson, 
eds.),  Vol.  U,  Pergamcm  Press  1957.  (Proceedings  of  the  Great  Plains 
Turbulence  Field  Program. ) 
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parameters  and  establish  the  variations  therefrom,  again  from  experimental 
observations. 

In  the  main  text  we  will  provide  the  background  material  to  justify 
the  various  rules  adopted  for  the  proposed  models  and  the  models  proper. 

Although  we  have  discussed  the  obscuring  aerosol,  the  remarks  are 
equally  applicable  to  explosive  dust;  that  is,  we  have  developed  a  semi- 
empirical  model  to  describe  the  cratering  phenomena  aud  subsequent  dis¬ 
persion  of  the  dust  cloud.  Once  this  cloud  is  "stabilized, "  its  dispersion 
is  handled  like  that  of  any  aerosol  cloud. 

the  report  is  divided  into  three  parts:  Part  L  -  Smokes  and  Obscuring 
•  irosol  Model;  Part  2  -  Explosive  Dust  Obscuring  Model;  and  Part  3  - 
Vehicular  Dust  Obscuring  Model.  Collectively,  the  three  models  constitute 
a  "dirty"  battlefield  atmospheric  model. 

this  writer  recognizes  chat  improvements  can  be  made  in  these  models; 
however,  additional  experiments  and  analyses  are  required  to  establish  the 
b-ases  for  the  improvements. 
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SMOKES  ANO  OBSCURING  AEROSQl  MQOEL 
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MODEL  DEVELOPMENT  SYNOPSIS 


1.1  INTRODUCTION 

la  Ref.  1  we  presented  a  semi-empirical  model  chat  can  predict  the 
spatial-temporal  characteristics  of  the  aerosol  cloud  produced  by  smoke 
and  other  obscuring  aerosol  sources  (artillery  shells,  mortars,  smoke  pots, 
and  aerosol  generators).  These  characteristics  are:  (1)  mean  concentra¬ 
tion,  (2)  path  length,  and  (3)  radiation  transfer.  In  addition,  Ref.  1 
provided  several  validations  of  the  proposed  model  where  the  detonated 
munitions  were  either  bulk  white  phosphorus  (WP)  or  phosphorus  suomunitions 
(VP  wicks) . 

Ref.  1  addressed  the  following  topics:  source  size  data;  fluctua¬ 
tions  effects;  and  the  rno.^el  structure  for  quasi-exo thermic  and  non- 
exothermic  obscurant  sources.  However,  topics  not  discussed  therein: 

(a)  A  treatment  of  the  still  air  case  (mean  wind  speed,  u  =■  0) 

(b)  Field  test  data 

(c)  Extinction  coefficient  data 

are  treated  in  this  report;  we  extend  here  the  model  presented  in  Ref.  1. 

In  the  remainder  of  this  section  we  summarize  the  discussion;  and  in  Sec.  2 
(2.1  to  2.5)  the  reader  will  find  a  complete  discussion  and  several  addition¬ 
al  validations. 

1.2  MODEL  EXTENSION  AND  MODIFICATION 

1.2.1  Concentration  Fluctuations 

The  experimental  concentration  data  exhibit  small  and  large  scale 
fluctuations  about  a  mean  value.  In  Ref.  1  we  indicated  that  horizontal 
wind-speed  variations  can  aecounc  for  the  overall  concentration  fluctua¬ 
tions;  however,  the  available  sampling  data  could,  at  best,  provide  an 
estimate  of  the  standard  deviation. 

^Zirkind,  R. ,  An  Obscuring  Aerosol  Dispersion  Model,  GRC  Report  231,  Vols. 

I  and  II,  December  1978. 
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Since  che  atmospheric  stability  descriptor  used  in  our  model  is 
the  Pasquill  Category  (Ret-  1),  we  have  selected  fluctuation  data  chat 
have  been  correlated  with  Pasquill  Category;  see  values  given  below. 

Pasquill  Category  A  B  C  D  b  F  G 
J(u)  1.13  1.2  1.6  1.33  1.08  0.89  0.51 

u'  /u  0.2  0.18  0.20  0.2  0.13  0.17  0.14 

cwlo 

where  U*  ^  is  che  root  mean  square  of  the  horizontal  wind-speed  fluctua¬ 
tion  u',  and  J ( u)  is  the  standard  deviation  of  u.  For  use  in  che  model, 

che  u'  value  should  be  used  to  obtain  an  estimate  of  che  small-scale 
Kilo 

fluctuations. 

The  larger  scale  fluctuations  of  interest  here  occur  with  a  period 
on  the  order  of  40  to  60  sec;  however,  data  available  from  field  tests 
have  not  been  analyzed.  These  fluctuations  are  responsible  for  che 
phenomenon  known  as  intermittency  (when  the  concentration  approaches  a 
near-zero  value  for  a  period  of  about  40  to  60  sec).  The  available  data 
indicate  that  the  effect  occurs  at  random  times  and,  therefore,  a  statis¬ 
tical  model  is  required  if  the  impact  on  sensors  is  to  be  treated  accurately. 

1.2.2  Aerosol  Dispersion 

The  basic  algorithms  to  find  the  cloud  growth  for  u  >  0  have  remained 

unchanged;  however,  che  rise  of  the  buoyant  cloud  has  been  clarified  and 

updated.  For  example,  when  the  atmospheric  stability  is  unstable  or  neutral 

then  the  final  rise,  AH  is 

max 

uH  — 1.61  (3.ox^)  v  u 

where  x*  =  14F^^  for  F  £  55  m^/sec^ 

xA  ■  34?^'**  for  F  >  55  m'Vsec'* 

in  which  x^  «  transition  distance  and  F  =  3.6  ;<  10  5  (heat  release/sec) 
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Similarly,  the  values  of  y  and  2,  the  Lateral  and  vertical  coordinates 
respectively,  have  been  given  as  continuous  functions  for  unstable  to  staole 
atmospheric  conditions;  see  Appendix  C. 


Lastly,  we  suggest  the  following  rules  be  used  for  still  air,  u  0  0: 
(1)  For  an  Instantaneous  source  the  maximum  rise  is 


H  «  1.89Q0*23  (m) 

oax 


(2)  For  a  maintained  source 


H  -  5F °-i5/;3/9 


Wne re  s  *  g((TT/.iZ)  +  ?]  i  T 


The  definitions  of  the  above  are  as  follows: 
g  s  gravitation  acceleration  (m/sec) 

=»  adiabatic  lapse  rate 
-T/2.Z  a  actual  lapse  rate 

T  »  ambient  absolute  temperature  ( *K) 


The  shape  of  the  aerosol  cloud  is  assumed  to  be  a  right  circular  cone  with 
a  vertex  angle  of  32*. 


In  Ref.  1  ve  omitted  the  effect  of  a  ground  or  elevated  inversion 
layer  on  the  growth  of  the  aerosol  cloud.  In  Sec.  2  the  rules  are  provided 
for  various  atmospheric  stability  conditions  and  are  incorporated  into  the 
model  structure.  Sec.  3. 


1.2.3  Model  Issues 

We  will  now  discuss  several  topics  that  require  iurtr.er  clarification 
and  probably  additional  research  and  development. 
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— -  - — " 


Source  Size 


The  experinenc.il  data  defining  the  initial  cloud  size  generated  by 
munitions  filled  with  bul-c  phosphorus  are  limited.  These  data  would  in¬ 
dicate  that  the  initial  cloud  dimensions  are  large;  e.g.,  the  IT  value  of 
the  radius  and  height  for  an  artillery-delivered  155-mm  shell  are  —3  m  and 
3  m,  respectively;  and  within  several  seconds  the  downwind  dimension,  x, 
can  be  an  order  of  magnitude  greater  than  che  calculated  value,  where 
x  -  wind-speed  *  time.  The  importance  of  this  parameter  is  twofold: 

(1)  it  gives  us  the  distribution  of  tne  source  material  for  buoyancy  cal¬ 
culations;  and  (2)  it  is  the  mathematical  representation  of  the  source, 
i.e.,  whether  Che  source  is  a  quasi-point  or  a  distributed  type.  If  che 
bulk  material  is  spread  over  a  large  area  then  che  buoyancy  will  be  less 
chan  if,  as  assumed  in  the  modei,  che  material  is  released  within  a  small 
area.  Wit"  respect  co  the  dispersion  modeling,  if  we  assume  a  distributed 
source,  then  che  cloud  development  would  be  similar  to  che  treatment  used 
for  submunlcions  like  the  white  phosphorus  wicks.  Since  che  modei  repro¬ 
duces  the  observed  cloud  geometry  quite  accurately  from  stacicaliy  fired 
munitions,  we  assert  chat  che  bulk  material  is  concentrated  in  the  vicinity 
of  che  detonation.  Also,  che  observed  cloud  at  early  times  (£  several 
seconds)  is  a  visible  phenomenon  generated  by  small  quantities  of  che  total 
fill  weight  and  che  coupling  of  the  high-explosive  burster  with  che  VP. 

This  topic  is  discussed  in  greater  detail  in  See.  2. 

This  issue  and  the  burn  time  require  additional  experimentation  and 
research  co  provide  a  sore  solid  base  to  this  and  ocher  models. 

For  che  white  phosphorus  subauni cions,  the  situation  is  similar, 
particularly  the  ground  dispersion  pattern  of  the  wicks  or  wedges  for  live 
firings.  Although  the  proposed  dispersion  pattern  simulation,  three  linear 
arrays  in  our  model,  yields  good  agreement  with  experiments,  we  strongly 
recommend  extensive  further  experimentation  to  obtain  statistical  data  of 
the  ground  dispersion  pattern  and  additional  burn  time  and  rate  data. 

* 

For  the  smoke  Obscuration  Model  this  parameter  is  important  since  the 
distribution  for  the  source  is  assumed  to  be  Gaussian  and  is  an 
initial  condition  for  the  diffusion  equation. 


*  * 


£  fflc lency 

The  literature  specifies  'he  fill  weight  of  che  active  basic 
chemicals;  however,  che  percentage  of  conversion  to  smoke  factor,  effi¬ 
ciency  is  not  available  except  in  a  qualitative  sense.  For  example,  we 
find  such  statements  as  "any  material  entering  a  pool  of  water  is  not  con¬ 
verted  to  smoke"  or  "fragments  will  be  buried  in  sofc  dirt",  evidently, 
this  should  be  quantified  by  experimentation.  In  che  interim  we  have 
assumed  chat  the  efficiency  for  hard  dirt  is  1007,;  however,  in  soft  dirt 
or  "muck"  the  values  are  352  and  50%  respectively. 

Extinction  Coefficient 

The  speccrai  extinction  coefficient  is  a  critical  parameter  to  estab¬ 
lish  the  optical  properties  of  the  obscuring  cloud;  however,  che  values  for 
the  different  obscurants  are  not  firmly  established.  For  example,  che 
values  are  dependent  on  the  particle  size  distributions.  Field  data  indicate 
that  the  size  distribution  varies  during  the  lifetime  of  che  obscuring  cloud, 
and  chat  che  distribution  descriptors,  mass  mean  diameter  (MHO)  and  geossetri: 
standard  deviation  differ  from  che  values  used  in  laboratory  experiments 

‘a 

co  obcs.n  the  speccrai  extinction  coef r'ic rents.  It  should  be  noted  chat 
these  and  related  topics  are  discussed  in  See.  1. 

Hence  an  in-depth  analysts  of  the  field  test  data  is  necessary  to 
establish  the  spectral  extinction  coefficient;  that  Is,  correlate  the 
laboratory  lata,  independent  concentration  and  particle  size  data,  and 
transaissometer  data. 

Lastly,  we  strongly  recossmend  that  an  analytical  and  measurement 
program  be  conducted  to  establish  the  iatermittency  properties  of  ailitarv 
smoke  clouds,  the  interaitteacy  rs  important  to  reliably  establish  electro- 
optical  sensor  performance,  especially  laser  systems. 
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BASIC  CONCEPTS 


2.1  GENERAL 

The  basic  objective  of  che  obscurant  model  is  to  determine,  to  some 
degree  of  accuracy,  sensor  performance  when  an  obscurant  cloud  exists  in 
che  field  of  view  of  the  sensor.  This  may  imply  che  temporal  history  of 

the  path  transmission,  cloud  radiance  (self  or  reflected),  or  other  param- 

★ 

ecers  that  affect  sensor  performance.  To  calculate  these  parameters  for 
different  obscuranc  sources  several  factors  are  required,  and  these  are  as 
r  o i lows : 

1.  Source  characteristics 

2.  Description  of  the  aerosols  produced  by  che  source 

3.  Spacial-cemporai  history  of  obscuranc  aerosol 
concentration 

4.  Environmental  interaction  wich  che  aerosol  cloud 
production  and  development 

With  che  exception  of  the  first,  the  remaining  factors  are  characterized 
by  che  turbulent  diffusion  of  source  material  in  the  environment.  In  che 
next  several  sections  we  present  a  discussion  or  each  factor. 

1.2  SOURCE  CHARACTERISTICS 

the  nomenclature,  fill  weight,  burst  radius,  and  ocher  technical 
data  far  inventory  munitions  are  >j?ven  in  U'S  Army  publications  iFM^-137, 
33,  etc.).  Several  of  these  data  sheets  are  given  in  Appendix  A.  Recent 
developmental  items  include  subauaicions  such  as  white  phosphorus  wedges 
and  wicks. 

Hence  we  now  have  a  spectrum  or  source?  varying  from  near  inscanca- 
eous  point  types  to  linear,  continuous  arrays,  the  importance  of  this 
lies  in  che  variety  of  mathematical  descriptions  needed  to  describe  the 
dispersion  history.  This  is  consid*.  -i  in  Sec.  2.5. 

a 

the  radiation  transfer  methodolog-w. 
referred  to  Ref.  1. 
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.  r  t  omitted  here;  the  reader  is 


2.2.1  Source  Size 


I 

I 

I 
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The  character  os  the  smoke  obscurant  is  of  interest  to  the  cloud  dis¬ 
persion,  We  consider  smoke  pots  and  fog  oil  generators  as  poinc  emitters. 
Until  further  data  are  available,  we  consider  that  munitions  filled  with 
wedges,  wicks,  or  plasticized  phosphorus  produce  ground  patterns  chat  can 
be  simulated  as  linear  arrays.  This  gap  in  our  knowledge  needs  to  be 
filled  by  further  explicit  tests. 

In  our  model,  we  assume  chac  artillery  shells  and  mortars  delivering 
the  bulk  white  phosphorus  are  quasi-point  sources;  the  source  is  in¬ 
variant  with  the  shell  and  mortar  size.  Since  this  macerial  is  a  major 
smoke  cloud  generator,  we  will  consider  it  here  in  greater  detail  chan  in 
Ref.  i  and  include  the  initial  cloud  development. 

Dolce  and  Mec2^  have  attempted  to  characterize  the  smoke  cloud  forma¬ 
tion  and  development  for  static  firings  and  artillery-delivered  bulk  white 
phosphorus.  We  will  first  examine  the  source  size  data.  Figure  1  presents 
the  Dolce  and  Metz  data  along  with  older  data,  the  half-width  of  the  source 
radius,  and  the  source  height,  r^.  I:  should  be  noted  chac  these 

values  refer  to  the  cloud  dimensions  and  not  the  ground  mass  distribution. 
Several  additional  values  are  added  from  the  Fall  197?  Dugvay  tests. 

The  artillery  curve  is  higher  than  the  mortar  cur/e.  We  will  now 
attempt  to  explain  this  and  ocher  data. 


In  Fig.  2  we  have  plotted  the  initial  and  subsequent  cloud  rise  for 
a  sO-aa  mortar  and  a  155-aaa  artillery  shell.  We  observe  that  the  plume 
centroid  has  a  different  time  history  than  the  initial  cloud  rise  and  sub¬ 
sequent  “source1*  neighc.  In  F.ef.  3  the  authors  identify  a  “weapon  phase" 
tise  it  ) .  This  time  corresponds  to  the  time  interval  from  initial  detona¬ 
tion  of  the  high  explosive  to  the  tise  the  cloud  becomes  white,  i.e..  loses 


Du  1 :  c ,  T .  2. ,  and  D.  ? .  Metz.  An  Ana l y s i s _ : he  $ soke  £  loud  3  a  1 4  frea 
August  237$  Jefferson  Proving  drouaS  $soWe'~Tea t ,  AMfCU  tec&Bie al'Hepo rt 
2Q1.  September  1977. 
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ids  orange  glow.  The  mean  values  of  the  weapon  phase  time  for  several 
munitions  are  given  below,  along  with  the  weights  cf  high  explosives  used 
to  disperse  the  white  phosphorus.  An  examination  of  the  times  clearly 


Munition 

HE  (lb) 

WP  (lb) 

cw  (St^ 

60  mm 

0.025 

0.75 

0.44 

81  mm 

0.08 

4.0  (1.75) 

0.56 

4.2  xn 

0.73 

.0.75 

0.87 

105  mm 

0.51 

3.83 

-0.6 

155  mm 

0.83 

15.6 

1.1 

indicates  that  they  correlate  with  the  mass  of  HE;  that  is,  the  mean  cloud 
temperatures  would  be  at  mean  ambient  for  times  equivalent  to  t^  (see 
Part  2).  Further,  the  power  released  by  the  HE  is  approximately  100  times 
greater  than  that  from  the  oxidation  of  the  white  phosphorus  within  the 
initial  25  msec.  At  later  times  (~  1  sec  intervals),  the  WP  material,  with 
the  exception  of  the  rare  large  piece,  will  be  oxidized  and  converted  to 
phosphoric  acid. 

Hence,  the  initial  rapid  rise  for  t  <,  5  sec  is  primarily  due  to  the 
HE  detonation  and  some  entrained  phosphorus  particles;  however,  the  effect 
on  the  vertical  cloud  rise  behavior  at  later  times  is  negligible.  This  may 
be  justified  on  energetic  considerations:  the  heat  release/unit  weight  is 
approximately  the  same  for  HE  and  WP,  whereas  the  ratio  WP/HE  varies  from 
8:1  to  30:1;  therefore,  the  total  impact  on  overall  buoyancy  is  small. 
However,  the  turbulent  flow  generated  by  the  HE  fireball  may  be  responsible 
for  the  early  history  of  a  W?  smoke  cloud. 

The  effect  of  the  large  difference  between  the  quantity  of  WP  and  HE 
in  the  two  weapons  is  now  clear.  For  60-mm,  the  small  quantity  of  HE,  the 
effect  is  small  and  the  centroid  is  not  different  than  the  "source"; 
however,  for  the  155-mm  the  explosive  effect  and  certainly  the  quantity  of 
WP  is  larger.  The  latter  is  responsible  for  the  significant  centroid  rise. 
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In  Ref.  3,  che  initial  dispersion  of  cloud  in  the  downwind  direc¬ 
tion  (x-axis)  is  not  addressed  except  that  o  =  for  wind  speeds 
<4  mps  (see  Fig.  1)  .  The  latter  statement  is  inconsequential  since 
the  data  indicate  that  che  cloud  length  at  the  end  of  the  weapon  phase  may 
be  in  excess  of  100  m;  see  Fig.  3.  Beyond  this  time  the  axial  growth  is 
characteristic  of  mean  wind-speed  magnitude;  i.e.,  x  =  u  *  t.  These 
results  are  not  dissimilar  for  static  and  dynamic  detonations;  e.g. ,  a 
static  firing  of  a  105-mm  had  a  ground  dispersion  of  —  61  m.  If  a  signif¬ 
icant  amount  of  bulk  phosphorus  is  dispersed,  these  data  are  difficult  to 
understand. 

A  simple  explanation  is  possible  of  the  overall  dispersion  observa¬ 
tions  discussed  in  Ref.  3  and  similar  tests;  that  is,  a  small  quantity  of 
material  can  generate  a  cloud  "apparent"  to  the  observer  or  camera.  Act¬ 
ually,  this  cloud  has  a  low  concentration  but  with  sufficient  glare  to  be 
interpreted  as  a  significant  cloud. 

If  che  phosphorus  is  dispersed  over  a  large  area,  then  the  engulfed 

air  volume  woulu  quickly  reduce  the  buoyancy  since  the  volume  increase  a 

2 

surface  area  =  (dispersion  width)  .  Hence  there  should  be  a  marked  etfect 
on  plume  rise;  however,  the  plume  rise  data  agree  with  experimental  and 
analytical  models  of  quasi-point  buoyant  plumes  to  be  described  later  and 
therefore  the  quantity  of  dispersed  material  is  assumed  to  be  small. 

In  view  of  the  limited  data  in  Ref.  1,  e.g.,  extremely  unstable  con¬ 
ditions  for  the  static  firings  and  other  similar  uncertainties  in  ocher 

test  results  on  cloud  geometry  and  dispersion,  we  cannot  resolve  the  above 

* 

problem.  A  definitive  experiment  is  required. 

2.2.2  Burn  Time 

4 

Salomon  and  Peterson  reported  their  findings  on  burn  rate  and  time 
for  several  types  of  smoke  munitions.  The  burn  rates  for  HC  and  VP  (wick/ 


It  should  be  noted  that  che  validations  performed  here  show  excellent  agree¬ 
ment  with  single  and  multiple  rounds  that  explicit,  valid  data  are  required 
prior  to  making  major  model  modification. 

4 

Salomon,  L.  L.,  and  E.  Peterson,  Testing  of  Smoke  Munitions  anf  Submunitions, 
Paper  H,  Proceedings  of  the  Smoke  Symposium  II,  April  1978,  DRCPM-SMK-T-003- 
787,  June  1978. 
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wedge)  are  shown  in  Figs.  4  and  5.  We  observe  that  the  rates 
tiveiy  unii’ora  tor  HC  and  WP  (wick)  but  nonunirorm  for  the  RP 
6.  Since  the  sample  is  small,  a  single  characteristic  may  be 
tor  phosphorus.  We  assume  Linear  burning  for  WP  and  the  rata 
given  in  Fig.  4. 


are  rala- 
wedge ,  Fig. 
appropriate 
for  HC  as 


The  burn  times  for  the  various  submunitions  a3  reported  by  Salomon 
4 

and  Peterson  are  given  in  Table  1.  It  should  be  noced  that  '!>  the  R? 
munitions  are  wedges  or  wicks  and  (2)  the  second  155ML  should  read  155H2. 


TA8LE  1 

SUBMUNITION  BURN  TIME 


Munlcion 

I ms. 

1 5  5  M 1 

HC 

155M1 

HC 

105  cm 

HC 

6"  Wick 

w? 

3"  Wick 

UP 

2.75"  Wick 

UP 

81  mm  (Navy) 

?P 

31  mm  (German) 

RP 

155  nan  (Navy) 

RP 

Time  (min) 

2.3 

1.3 

2.0 

6 . 5 

7.3 

4.3 

4.3 

3.6 

6.4 


The  burn  time  for  several  foreign  rounds  are  given  in  Table  2 
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FRACTION  OF  BURN  TIME 


Figure  6. 


Submuni  cior.  Burn  Rate,  Navy  155  mm  Wedge  (RP) 
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Salomon  and  Peteraon,  op.  cit. 
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First,  if  the  burn  times  in  Table  2  are  for  bulk  white  phosphorus 
(WP)  and  not  plasticized  white  phosphorus  (PWP)  then  they  are  not  in 
agreement  with  the  foreign  literature.  Second,  the  weight  of  WP  is  re¬ 
latively  low  to  burn  for  such  long  periods;  that  is,  on  the  basis  of 
firing  replacement  rates  the  estimate  burn  times  are  extremely  long.  The 
results  for  the  82-mm  mortar  are  longer  than  for  che  US  81-tnm  mortar  where 
80%  is  oxidized  in  £2  sec  and  the  balance  in  £10  sec.3 


Lastly,  the  experimental  data  on  concentration  histories  do  not 
agree  with  Salomon  and  Peterson's  data. 


In  the  Dugway  tests  held  in  the  fall  of  1977,  where  inventory  smoke 
munitions  were  tested,  the  4.2"  WP  mortar  like  the  81  mm  is  primarily 
oxidized  rapidly  (£2  sec).  To  support  this  we  note  that  the  concentration 

3 

is  -0.8  g/m  at  1  sec,  55  m  from  the  origin;  within  10  sec  the  entire 
quantity  is  oxidized. 

The  impact  of  the  explosive-generated  thermal  environment  on  the 
kinetics  of  bulk  phosphorus  does  not  appear  in  the  literature. 


With  respect  co  HC  canisters  the  following  burn  time  data  were 
provided  by  the  Program  Manager  (PM) /smoke-obscurants : 


Munition 
155  mm  Ml 
155  mm  M2 
105  mm  Ml 


Burn  Time 
120-240  sec 
60-240  sec 
180  sec 


It  should  be  noted  that  the  M2  has  —  50%  the  fill  weight  of  the  Ml.  The 
above  do  not  differ  from  the  earlier  results,  but  the  obscuration  of  smoke 
pot  emission  is  simple  compared  to  that  of  an  explosive  munitions  like  a 
mortar  or  artillery  round. 

^Private  communication,  CSL  Edgewood  Arsenal,  1  December  1978. 
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The  estimates  made  by  Salomon  and  Peterson  are  qualified  by  them 
and  relate  to  mass  loss.  In  the  neat  paragraph  we  discuss  the  burn  time. 

The  fill  weight  specified  in  a  standard  munitions  manual  gives  r.he 
basic  weight  of  basic  chemicals  chat  can  generate  an  obscuring  aerosol 
cloud,  cor  munitions  delivered  by  artillery,  mortars,  or  rockecs,  the  fuze 
altitude  (or  type)  is  important  since  for  a  near  or  on-ground  burst  the 
ground  hardness  plays  a  role  in  conversion  efficiency.  If  the  ground  is 
soft  or  adjacent  to  water  then  part  of  the  fill  material  may  be  embedded 
in  the  subsurface  material  and  not  available  to  the  aerosol  production 
process.  Also  an  efficiency  factor  appears  to  exist  for  smoke  pots;  thac 
is,  the  encire  material  is  not  converted  into  an  obscuring  aerosol  cloud. 

Soviet  pots  appear  to  have  efficiencies  as  low  as  402,  a  low  value. 
Observations  by  US  sources  are  unknown  to  this  writer;  however,  highly 
efficient  smoke  pets  should  not  be  difficult  to  produce. 

Thus  we  assume  thac  for  phosphorus  munitions  che  fill  efficiency 
is  tc  1  for  hard  pack  dirt,  832  for  loose  dirt,  e.g.,  plowed  field,  and 
£502  for  "muck."  For  smoke  pots,  it  is  again  assumed  that  if  the  pot 
ignites  with  che  outlets  exposed  to  che  atmosphere  then  che  efficiency 
m;l;  actually,  che  value,  on  che  average,  is  less.  Further,  the  per¬ 
centage  of  upright  cauisters  is  <100%;  the  exact  number  is  not  available. 

2.2.3  Yield 

The  yield  factor  is  defined  as  che  final  mass  to  che  initial  mass 
available  for  smoke  generation.  Johnson  and  Forney0  published  yield 
results  for  various  smoke  materials.  These  are  summarized  in  algebraic 
form,  where  &H  »  relative  humidity; 

Y (WP)  -  3.8  ♦  0.003  <RH%  -  1G)1'67 
Y(ZnCU)  -  1  +  0.051  (RK%  -  5)°‘85 
Y(A1G13/H.,5Q4)  -  2.3  *  0.016  (RH%  -  20)1*25 
Y(Fog  Oil)  -  1 


Johnson,  M.  C.  and  P.  0.  Forney,  The  Effectiveness  of  Obscuring  Smokes. 
0RG  tdgewood  Arsenal,  1972  (unpublished). 
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Recently  Rubel 1  published  theoretical  yield  results  for  white 
phosphorus  which  agree  wich  che  Johnson  and  Forney  theoretical  results. 

On  the  other  hand,  field  measurements  of  concentrations  would  indicate 
that  lower  values  give  better  agreement:  Y(VP)  a  2.9  +  0.CO3  (RH.»  - 
iG)1'*3'.  Ac  -»02  relative  humidity  che  difference  is  —1.74,  a  major 
difference  chat  could  noc  be  accounted  for  in  observed  and  calculated 
values  of  aerosol  concentration  unless  several  factors  compensate  for 
che  difference. 

For  completeness,  we  cite  cwo  recent  studies  that  consider  the  physico- 

8  9 

chemical  problem  of  phosphorus/phosphoric  acid  ’  ,  where  the  former  con¬ 
siders  particle  size  and  kinetic  times. 

With  respect  to  kinetic  times,  the  various  stages  are  modeled  with 
the  exception  of  the  oxidation  (combustion  process),  which  is  responsible 
for  the  heat  release  and  subsequent  conversion  to  phosphoric  acid.  Clearly, 
the  large  variations  in  burning  times  should  have  an  effect  upon  che  over¬ 
all  particle  size  distribution,  etc.  Farther,  che  environment  created  by 
che  high-explosive  bursters  should  alter  the  initial  amount  of  oxygen  and 
aoiscure  available  for  chemical  reactions  and  che  state  variable  of 
pressure  and  temperature. 

2.3  data  problems 

2.3.1  Extinction  Coefficients 

One  of  the  agonizing  problems  for  the  modeler  (and  analyst;  la  the 
magnitude  of  the  spectral  extinction  coefficients  (absorption  and  scattering; 


‘ Rubel,  3.  0.,  Predicting  the  Dranlet  Size  ned  Yield  Factors  of  a  Phosphorus 
Smoke  as  a  Function  of  Droplet  C-.  osition  and  Asfeient  Relative  Humidity 
under  Tact  leal  Conditions.  AR.CSL-it-7305?,  Chemical.  Systems  laboratory, 

APG,  November  19?B. 


Rubel,  C.  0..  An  Aerosol  Kinetic  Model  for  the  Condensational  Growth  of  a 
Phosphorus  smoke,  Proe.  of  Smoke  Symposium  til.  1979. 

^Friekel,  S.  H. .  C.  0.  Rubel.  and  g.  V.  Stuehiag,  Relative  Humidity  Dependence 
of  the  Infrared  Extinction  by  Aerosol  Clouds  of  Phosphoric  Ac .4.  Proc -  of 
Smoke  Symposium  III.  1979. 


for  the  different  obscuring  aerosols.  For  example  Mil ham ^  gives  the 

I 

extinction  coefficient  sc  1.06  urn  for  HC/RP  as  1.5  m“/g,  whereas  derived 
values  from  the  field  data  are  0.93  (HC)  end  0.59  (UP) ;  see  Table  3. 

The  differences  between  laboratory  and  field  trial  data  are  difficult  to 
rationalize  for  several  reasons: 

1.  The  concentration  history  at  the  transtais soma ter  is  dif¬ 
ferent  chan  che  values  at  the  sampler,  therefore,  cor¬ 
relation  is  noc  probable. 

I.  The  particle  si2e  distribution  in  che  field  changes  with 
time,  and,  therefore,  che  distributions  characteristic  of 
the  extinction  coefficient  are  unknown. 

J.  droadband  devices  are  utilized;  a.g.,  the  9.75-wn  trans- 
missomecer  actually  covers  from  —  9  to  11  urn;  see  Table 
3.  Vow  che  integrated  extinction  coefficient  is  a 
function  of  che  detector  spectral  response  and  che  spectral 
extinction  coefficient  so  chat  che  "effective”  extinction 
coefficient  depends  on  che  C  •  L  value,  hence,  che  sug¬ 
gested  field  values  for  "a"  is  not  obvious. 

Conversely  che  ability  to  validate  che  cransmissomecer  measurements  from 
a  predicted  C  •  L  vs.  time  cur/e  will  be  difficult  due  to  the  repeated 
deviations  in  che  experimentally  determined  extinction  coefficients. 

We  will  now  examine  the  laboratory  and  field  data,  and  che  impact 
of  particle  size  on  the  accuracy  of  spectral  extinction  coefficients 
when  a  transaissometer  is  used. 

Rubel  (Ref.  7)  reports  that  at  36S  relative  humidity  he  found  for 
the  phosphorus  particles  that  the  mass  median  diameter  -  1.14  ya  and 
geometric  standard  deviation,  't  *  1.34,  where  the  particle  size  is 

‘^Milhaa.  M.  ,  A  Catalog  of  Optical  Sstinecion  Data  for  Various  Aerosols/ 
Smokes.  Sdgewood  Arsenal  Report  EC-SF-77001.  June  1976. 
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defined  by  a  log-normal  distribution.  An  analysis  of  several  Dugway  crisis 
was  performed  here  for  particle  size  characteristics  where  the  relacive 
humidity  varied  between  32%  and  36%.  The  results  -’re  given  in  Table  4  for 
the  average  values.  Clearly  the  spread  is  greater  chan  can  be  attributed 


TABLE  4 

BULK  PHOSPHORUS  SMOKE  PARTICLE  SIZE 
DISTRIBUTION  PARAMETERS 


Teat 

Muni  cion 

MMD  („m) 

J 

_i 

RH  (%) 

DPI-G02-T34 

153  mm 

3.39 

2.13 

35.0 

0PI-C02--T16 

4.2  ' 

1.14 

1.76 

34.2 

OPI-0Q5-I16 

130 

3.38 

2.0 

32.6 

DP1-0Q2-T24 

155  mm 

2.65 

2.08 

36.2 

0PI-CQ5-T12 

122  mm 

1.795 

1.33 

31.8 

to  statistical  variation.  Further,  the  observed  field  daca  show  time 
dependence  on  size  distribution  where  settling  or  coagulation  processes 
are  insignificant  if  not  nonexistent;  see  Fig.  7. 

Hence,  it  the  aerosol  distributions  for  w?  are  time  dependent  and 
munition  related,  then  the  use  of  c ransmis some t era  to  establish  concentra¬ 
tion  histories  is  subject  to  further  validation.  *e  now  discuss  this 
important  topic. 

The  determination  of  the  aerosol  concentration  in  an  obscurant  cloud 
has  been  aade  by  aerosol  ^photometric)  samplers  and  craasmissoseters .  There 
is  an  extensive  literature  on  sampling  techniques  and  their  limitations  and. 
therefore,  samplers  arc  not  discussed  here  except  to  note  chat  the  system 
noise  spectrum  for  the  devices  used  in  the  field  by  test  personnel  would 
be  extremely  valuable.  The  reason  is  the  desirability  of  separating  the 
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contribution  of  atmospheric  fluctuations  in  the  temporal  variation  of  the 

o 

C  *  L  curve  from  system  noise,  where  C  =  concentration  (g/nT)  and  L  =  path 
length  in  cloud  (m) ,  This  separation  is  required,  as  the  data  clearly 
indicate  noise  fluctuations.  The  use  of  transmission  measurements  in 
order  to  derive  particle  concentration  requires  some  discussion. 

First,  there  is  r  question  about  the  current  value  of  the  spectral 

2 

extinction  coefficient.  The  parameter  C  ♦  L  (g/xa  )  is  found  from  the 
relation 


C  •  L  «*  in  T/a,  (1) 

where  =  spectral  mass  extinction  coefficent 

T  a  fractional  transmission  through  the  cloud 

The  standard  procedure  is  to  utilize  laboratory  measurements.  These  may 
not  simulate  free  space  conditions  as  discussed  earlier.  We  will  now 
attempt  to  show  that  this  may  not  be  serious  at  some  wavelengths.  For  this 
we  will  follow  the  treatment  of  Chylek.^1 

A  relation  between  the  volume  extinction  coefficient,  k,  and  the 
mass  concentration,  M,  can  be  written  as 


♦  n(r)  dr _ _ 

(m,  r,  X)n(r)  dr 


k 


(2) 


where  p  a  bulk  density 

r  "  particle  radius 

n ( r)  «  number  of  concentrations  per  unit  radius  change 

Q  =  extinction  efficiency  factor 
ext  J 


^Chylek,  P.,  et  al.,  "Infrared  Extinction  and  the  Mass  Concentration  of 
Atmospheric  Aerosols,"  Atmospheric  Environment,  13,  1979,  pp.  169-173. 
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Since  p  and  m  (the  refractive  index)  are  not  significantly  different  in  time 
and  position,  the  relation  (Eq.  2)  strongly  depends  on  n(r).  On  the  other 
hand,  if  the  bracketed  expression  is  equal  to  a  constant,  h,  for  a  specified 
n(r)  then  it  follows  chat 

M  -  k(A)4p/3h  (3) 

A  sufficient  condition  for  this  situation  is 


Qexc(r,  i)  -  hr 


(4) 


or,  the  extinction  efficiency  factor  is  proportional  to  particle  size. 
Let  x  -  2nr/X  be  the  size  parameter;  then 


.  A 

Q  „  =  h  •  *  X  C  •  X 

ext  2tt 


(5) 


and,  therefore. 


M  =  k(\)  •  (6) 

where  c  is  the  slope  of  Qext  vs.  x.  Clearly  if  c  is  a  constant  or  slowly 
varying  for  a  given  range  of  particle  sizes,  then  the  mass  concentration 
is  a  constant  if  k ( A)  *  X  The  latter  occurs  for  moderate-size  particles, 
<1.0  pm.  In  Ref.  11,  Mie  calculations  for  a  realistic  particle  show  that 
Q  vs.  r  (pm)  is  linear  at  A  =  0.55  pm  for  particle  radii  £  0.5  pm  and  for 
r  <  5  pm  at  A  »  11  pm. 

The  range  of  applicability  is  defined  by  the  ratio 


i»00 

R(r  ,  a)  =  3/  r2Qn(r)  dr/4  /  r3n(r)  dr  (7) 

a  •'o  Jo 
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where  r  is  the  geometric  mean  radius  and  a  is  the  geometric  standard 
a. 

deviation;  that  is,  if  R  *  constant  for  a  range  of  r  ,  then  k(,\)  and  mass 

a 

concentration  are  independent  of  size  distribution.  (Note  for  a  fixed  a, 

R  =  3  c/2  E  constant.) 

For  particle  distributions  of  interest  (smoke,  aerosols)  where  r 

Si 

is  of  the  order  of  0.2  pm,  then  at  a  =  11  pm,  R  =  constant;  however,  for 
A  *  0.55  pm,  the  variation  in  R  is  ~100,  hence,  mass  concentration  at 
A  *  11  pm  is  known  to  high  accuracy  whereas,  in  the  visible  and  near  in¬ 
frared  (A  <  2  pm)  the  error  may  be  as  large  as  a  factor  —100.  It  should 
be  noted  that  for  particle  radii  <  3  pm,  the  error  at  A  cc  4  pm  is  similar 
to  A  =  11  pm,  that  is,  relatively  small. 

In  view  of  the  above,  the  use  of  transmission  measurements  to  establish 
mass  concentration  must  be  exercised  with  caution  if  accurate  values  are  to 
be  derived;  that  is,  the  size  distribution  must  be  known  and  the  appropriate 
wavelength  used  for  the  transmissometer . 

With  the  above  background  and  data,  let  us  examine  the  differences 
that  may  exist  between  the  field  test  derived  extinction  coefficients  and 
the  laboratory  values. 

In  Table  4  we  indicate  estimates  of  the  MMD  and  a  for  several  events 

g 

where  bulk  WP  was  used.  The  variation  in  MMD  was  from  1.14  to  3.39  and 

1.3 

correspondingly  1.7  tc  2.13  for  the  value  of  ar  -  Milham  et  al.  examined 

g 

the  variation  in  the  spectral  extinction  coefficient  for  O-phosphoric  acid 
for  two  values  of  MMD  and  0 ,  and  their  results  are  shown  in  Fig.  8  along 
with  the  value  for  WP  smoke  (35%  relative  humidity)  in  the  7-14  pm  band. 

For  the  3-5  pm  band  they  conclude  that  WP  smoke  is  no  different  than  0- 
phosphoric  acid  and,  therefore,  we  can  conclude  that  size  distribution 
properties  can  have  a  significant  effect  on  the  integrated  or  spectral 

12  " 

Milham,  M.  E.,  et  al.,  New  Findings  on  the  Nature  of  WP/RP  Smokes,  \ 

CSL,  ARMCOM  Report  ARCSL-TR-77067*.  July  197  7.  ■  '> 

•(  i 
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Figure  8 


Extinction  Coefficient  vs.  Wet/eiength  ( - bulk  WP)  :  (a) 

MMD  -  1  um,  cr  -1.5;  (b)  MMD  “3.3,  2  urn,  jg  =  2 
(Milham,  M.  E?,  et  al.,  Wav  Findings  on  the  Nature  of  WP/RP 
Smokes,  CSL,  ARMCOM  Report  ARCSL-TR-77067 ,  July  1977.) 
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extinction  coefficient.  The  field  value  at  3.4  um  is  given  as  0.22  t 

0.04  (1.t)  (m“/g),  which  does  net  agree  with  Fig.  8  in  absolute  value  or 

2  1 

variance.  At  9.75  urn  the  field  trial  result  is— 0.32  m  /g,  or  0.27  a"/g. 
which  again  is  lower  than  the  laboratory  values.  The  effect  of  size 
distribution  on  O-phosphoric  acid,  which  should  be  a  good  approximation 
to  WP  in  the  9-11  un  band,  shows  a  marked  difference  in  the  9-11  'an 
spectral  region. 

There  is  a  real  need  to  reduce  the  field  data  on  a  systematic  basis; 
that  is,  to  calculate  the  effective  extinction  coefficient  as  a  function  of 
time  using  the  time  variation  in  particle  size  distribution  properties 
and  C  •  L  and  then  compare  the  results  with  laboratory  observations. 

2.3.2  Field  Data 

There  are  many  probiv.ms  related  to  the  published  data  and  are  best 

illustrated  by  the  differences  between  Trials  9  and  10  of  the  Dugway 
13 

tests  called  Inventory  Smoke  Munition  Test  (Phase  Ha) .  These  are  tests 
of  the  4.2"  WP  mortar. 

In  Trial  DP1-002-T-9  a  single  4.2"  mortar  was  detonated  at  76  m 

* 

from  the  aerosol  sampling  line  under  the  following  conditions: 
u  (m/sec)  »  3.8 
Azimuth  wind  angle,  z  =  21.6* 

Pasquill  Category  *  C 

The  published  cloud  dimensions  as  a  function  of  time  are  given  as 
follows 


13 

US  Army  D.P.G.,  Inventory  Smoke  Munition  Test  (Phase  Ila).  Final  Test 
Report,  Report  0PG-TP-77-315 ,  June  1978. 

One  should  note  that  locations  of  the  detonation  lines  from  the 
sampler  line  are  uncertain  far  some  tests  in  spite  of  the  fact  that 
the  test  array  is  a  fixed  geometry. 
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c  (sec) 

Length  (m)  Width  (m) 

Height 

_ 

Measured 

Calculated 

Measured  Calculated 

0 

3  9 

— 

5 

— 

10 

17  23 

— 

30 

— 

20 

23  33 

27 

58 

43 

30 

27  43 

37 

87 

30 

40 

Plume  aloft 

The  face 

that  the  samplers  registered 

an  initial 

value  of  C  •  L  =* 

0.3  g/m ~ 

at  cz  50 

sec  does  not  agree  with  the  length  or  the 

"plume  was  aloft 

"  after 

40  sec. 

Similarly,  the  dosage  curve 

indicates  a 

width  of  100  m. 

The  height 

values  agree  reasonably  with  the  calculated  values. 

If  che  reported  azimuthal  variation  is  superimposed  on  the  expected 
value  of  the  cloud  width  at  the  sampling  line,  i.e.,  at  76  m,  then  the 
"sampler  width"  would  indeed  be  100  m.  We  observe  a  cloud  that  is  meander¬ 
ing  rather  strongly  in  the  azimuth  direction,  and,  therefore,  some  delay 
to  reach  the  sampling  line. 

In  Fig.  ?  we  present  the  C  •  L  history.  The  calculated  values  prior 
to  50  sec  are  greater  since  the  estimated  mean  arrival  time  =  76/3.3  ■  20 
sec.  The  possible  explanation  is  chat  the  cloud  did  rise  and  then  fall 
after  '*0  sec;  however,  the  plume  would  have  to  have  been  significantly 
longer  than  the  reported  observations .  On  the  other  hand.  Trial  10  is 
relatively  well  behaved  and  agrees  with  calculations. 

The  non-zero  values  prior  to  50  sec  arid  beyond  -  150  sec  are  to 
be  considered  as  either  noise  or  other  contamination  such  as  dust. 

Finally,  it  should  be  noted  chat  the  C  •  L  results  derived  from 
the  transaissoaecer  measurements  are  not  meaningful,  particularly  in 
light  of  agreement  between  saapler  observations  and  model  estimates  shown 
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2.4  AEROSOL  DIFFUSION 

The  diffusion  of  the  aerosol  cloud  is  of  primary  concern  here  since 
it  determines  the  spatial-temporal  history  of  the  airborne  material.  Ex¬ 
tensive  study  of  this  subject,  theoretical  and  experimental,  has  been  in 
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progress  for  the  past  fifty  years.  ’  Our  basic  interest  here  is  to 
adequately  describe  the  obscurant  aerosol  generated  by  a  smoke  munition, 
high  explosive,  or  smoke  generator.  With  the  exception  of  fog  oil,  the 
other  sources  release  large  quantities  of  thermal  energy  which  dominate 
the  early  history  of  the  diffusion  process;  that  is,  we  are  interested 
in  both  buoyant  plumes  (clouds)  and  turbulent  (momentum  exchange)  diffused 
clouds . 

The  general  behavior  of  the  clouds  depends  upon  the  atmospheric 
environment  (temperature  gradient,  wind  velocities  and  gradient,  surface 
roughness,  etc.),  or  the  stability  of  the  atmosphere.  Characterization 
of  the  atmospheric  stability  has  been  studied  by  many  investigators  and 
for  use  here  we  adopt  the  Pasquill  stability  definitions.  The  reason 
is  that  it  depends  primarily  upon  insolation  and  mean  wind  speed  at  10  m. 
However,  it  suffers  from  a  deficiency,  i.e.,  absence  of  surface  roughness. 
We  will  first  give  the  general  definitions  and  then  correct  for  surface 
roughness . 

The  three  basic  stability  classes  are  (1)  unstable  -  Pasquill  Cate¬ 
gories  A,  B,  and  C;  (2)  neutral  -  Pasquill  Category  D;  and  (3)  stable  - 
Pasquill  Categories  E,  F,  and  C,  where  G  is  the  extremely  stable  situation. 
A  simple  summary  is  given  in  Table  5,  and  a  more  precise  method  to  deter¬ 
mine  the  class  is  given  in  Appendix  B. 
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TABLE  5 

PASQUILL  STABILITY  CATEGORIES 


Night 


Surface  wind 
speed  at  10  m 
(m  sec  -1-) 

Insolation 

Thinly  overcast 
or  >4/8  low 
cloud 

<3/8  cloud 

Strong 

Moderate 

Slight 

<2 

A 

A-B 

3 

_ 

_ _ 

2-3 

A-B 

B 

C 

E 

F 

3-5 

B 

B-C 

C 

D 

E 

5-6 

C 

C-D 

D 

D 

D 

6 

C 

D 

D 

D 

D 

Classically  che  stability  definitions  are  given  by  the  environmental 
lapse  rata,  unstable  =  lapse,  neutral  =  neutral,  and  stable  =  inversion; 
that  is,  if  is  defined  as  the  potential  temperature  and  T  the  adiabatic 
lapse  race  (0.0984C/10  m)  ,  then 


'  2 


(3) 


where  3T/3z  =  T  (10  m)  -  T  (0.5  ra) .  Hence,  if  9  <  0  we  have  an  unscable 
condition.  Similarly,  5  ^  -  0,  .'0  we  have  a  neucrai  or  inversion  condition. 

For  convenience,  we  introduce  here  a  parameter  needed  later,  the 
stability  parameter,  s,  where 


s 


l 

T  3z 


(9) 


where  g  and  T  are  gravity  acceleration  and  ambient  ground  temperature  (K) , 
respectively.  The  parameter  can  be  interpreted  as  the  restoring  accelera¬ 
tion  per  unit  vertical  displacement  for  adiabatic  motion  in  a  stratified 
atmosphere. 
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A  somewhac  more  elaborate  scheme  of  Pasquiil  categories  is  reproduced 
in  table  6,  where  negacive  values  are  equivalent  to  instability,  etc. 

In  Fig.  10,  the  seasonal  variations  of  the  Pasquili  categories  for 
the  Netherlands  are  depicted,  and  one  observes  that  Category  A  does  not 
occur  frequently. 

to  show  the  diurnal  variation  of  the  atmospheric  stability,  we  present 
in  Table  7  data  from  the  Great  Plains,  Nebraska,  reported  in  Ref.  2.  We 
note  that  the  frequency  of  Category  A  is  low  and,  in  general,  unstable 
conditions  occur  primarily  during  the  day. 


The  stability  definition  of  Pasquili  (see  Appendix  3)  depends,  as  staced 
earlier,  on  wind  speed  and  insolation  and  is  independent  of  surface  roughness 
and  ground  cover  height.  This  can  be  accomplished  by  introducing  the 
Monin-Obukhov  "mixing"  length,  L  (m) ,  defined  as 


su /it 


kU/T0)«,/V> 


,  !Z 


where  =  friction  velocity  (m/sec) 

g  =  acceleration  of  gravity  (.a/sec') 


0 


a  scandara  temperature 


k  a  Karman  constant 
u  mean  wind  speed 


T  »  mean  temperature 
q  -  vertical  turbulent  heat  current 
e  ■  -»  air  specific  heat  capacity,  density 

r 


(10) 
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KEY  USED  TO  ESTIMATE  THE  CONTINUOUS  STABILITY  CLASSIFICATION 
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Flours  10.  Prevalence  ai  Paaquill  Stability  Categories 
Over  th«  Netherlands* 


S-eproduced  froa  C.  H.  Stroa.  "Atmospheric  dispersion  oi  Stse!<  Efiluents 
?.  13?.  Air  Pollution  (A.  C.  Stern.  14.)  lad  £d.  Academic  Press,  1963. 


TABLE  7 


DIURNAL  VARIATION  OF  ATMOSPHERIC  STABILITY 
GREAT  PLAINS,  NEBRASKA* 


Month  Hour  of  Mean  Local  Time 


&  Day 

00 

02 

04 

06  ! 

08 

10 

12 

la 

16 

18 

|  20 

*S  A 

3  Aug 

1 

, 

D 

H 

G 

9  Aug 

G 

E 

D 

;  d 

C 

C 

c 

D 

E  ! 

j? 
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13  Aug 

- 

- 

F 

£ 

i  3 

B 

a 

A 

D 

G  i 

H 

G 

1A  Aug 

H 

G 

- 

- 

i 

i 

1 

1 

, 

18  Aug 

- 

- 

- 

- 

i 

i 

- 

- 

- 

- 

H 

!  K 

1 

H 

19  Aug 

H 

H 

G 

G 

1 

.*» 

A 

A 

A 

B 

H 

!  H 

H 

23  Aug 

- 

- 

G 

£ 

i 

!  G 

A 

A 

3 

- 

- 

!  _ 

i 

- 

24  Aug 

- 

- 

- 

- 

< 

- 

a 

C 

D 

F 

F 

D 

25  Aug 

E 

D 

Q 

0 

:  D 

C 

c 

C 

D 

£ 

G 

- 

31  Aug 

- 

- 

E 

D 

•'  C 

A 

A 

B 

C 

D 

,  F 

j 

E 

1  Sept 

c 

F 

£ 

D 

1 

1 

8 

- 

- 

- 

- 

- 

- 

- 

7  Sepc 

- 

- 

H 

G 
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3 

3 

3 

C 
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£ 

F 

3  Sepc 
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A 
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Source:  Lettau  and  Davidson,  Ref.  2. 


layer ; 


The  value  of  L  is  characteristic  or  the  "dynamic 
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Tf'z)  »  constant.  we  have  neutral  strat  1 f  icat  Ion  L  *  «;  arid  i:  q  <  0,  •» 
we  have  L  •  0,  <0,  respectively.  One  can  show  that  the  rise,  1 z.Lx.  or 
smoke  cloud  is  -L;  that  is,  the  ratio  of  the  standard  deviations  of 
vertical  and  horirontii  wind  speed  fluctuations  *  and  r  respectively, 
~  5(2/1.).  Hence,  this  quantity  has  si gni: Icance ;  however,  the  reader  i 


referred  to  Csanady  for  a  complete  discussion. 


Colder*^  establi.-.hed  a 


and  Pasquill  stability: 


relation  between  the  Mon in-Obukov  length. 


in 


t  L 


'10 


:(s) 


and 


f(s)  -  -a.  vi  +•  b  jsj  -) 


where  s  *  stability  class  defined  in  Taoie  6 
a  »  *.  b  •  1.3.  e  •  0.33,  d  -  0.1165* 
2^  a  the  surface  roughness  in  meters 


The  surface  roughness,  2,,  vac.es  from  -10  J  a  1  sand)  to  —0.10  m  >.  tal 

0 

vegetation).  Xung  and  Letcau  have  found  an  empirical  relation  between 
2,,  and  plant  height,  : 


log  e .»  »  -l.ld  *  l .  I  *  log  h 
O  c 

Several  values  of  jsq.  h^  are  given  in  Table  <s. 

The  effect  of  surface  roughness  is  to  increase  or  reduce  in  some 
instances  the  stability  class;  see  Fig.  11. 


*  Csanady ,  C.  T.  ■  Turbulent  3  if  fusion  m 
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2.5  CLOUD  HISTORY 


The  clouds  formed  from  various  sources  rise  because  of  their  buoyancy 
and/or  by  momentum  exchange.  Since  our  initial  work,  numerous  validations 
have  provided  confidence  in  the  basic  formulations  for  thermal  sources 
(e.g.,  WP)  and  nonthermal  sources  (fog  oil).  Hence,  we  will  discuss  these 
formulations  and  some  modifications;  we  first  consider  the  nonthermal. 

2.5.1  Non-Exo thermic  Munitions 

Although  the  number  of  pure  non-exo thermic  munitions,  e.g.,  fog  oil, 
are  few,  an  HC  mixture  with  2a  or  ZnO  is  a  quasi-exo thermic  material  because 
of  its  long  burning  time  and,  therefore,  the  cloud  rise  behavior  after  the 
short  buoyancy  phase  is  described  as  momentum  exchange.  Also,  the  lateral 
diffusion  appears  to  be  unaffected  by  any  thermal  effects. 

The  cloud  size  parameters  determined  from  experimental  observations 
are  given  in  Table  9. 


TABLE  9 

CLOUD  SIZE  PARAMETERS  FOR  NON-EXOTHERMIC  MUNITIONS 


Stability  Category 
Pasquill  Turner 

Width, 

y  (m) 

Height, 

.  z  (m) 

Remarks 

A 

-3 

9.1 

+ 

0. 419x 

2.73 

j. 

0. 137x 

Sunny  day 

B 

-2 

9.1 

+ 

0. 328x 

2.73 

+ 

O.llx 

Day,  broken  cloud 

C 

-1 

9.1 

+ 

0. 238x 

2.  73 

+ 

0.073x 

Overcast  day/night 

D 

0 

9.1 

+ 

X 

o 

C'J 

o 

2.73 

+ 

Q.06x 

Neutral 

E 

1 

9.1 

+ 

0. 18x 

2.73 

+ 

0.055x 

Evening/ early  am 

F 

2 

9.1 

+ 

0. 146x 

2.73 

+ 

0. 046x 

Evening/early  am 

The  general  shape  of  the  cloud  is  a  semi-cone,  and  the  slope  is  given 
by  z/x.  Soviet  measurements  for  a  smoke  pot  (AT  —  2t'C),  terrain  roughness 


2q  equal  to  0.4  cm,  gave  the  following  cloud  slopes  in  terms  of  the  Monin- 
ObuRov  length: 


tan  a  =  -0.85/L  +  0.07  (unstable) 

=  0.07  (neutral/stable) 


(13) 


When  z„  =  0.4  cm  and  the  Pasquill  stability  is  A  -  B,  the  value  of 
-1  J  -1 

L  <  -0.07  (m  '*■)  and  the  cloud  slope  —  0.11.  Since  the  from  Table  8 
may  differ  somewhat,  the  value  of  0.11  is  a  good  average  value  (0.09-0.17) 
for  Category  B.  This  indicates  a  good  correlation  between  the  Pasquill 
Turner  and  Monin-Obuxov  stability  parameters. 


One  should  recognize  that  the  discrete  values  given  in  Table  10  for 
the  constants  multiplying  "x"  are  continuous  for  0.5  <  u  <  5  m/sec;  see 
Appendix  C.  The  low-speed  values  are  omitted  since  the  behavior  at  or 


near  u  =  0  should  be  defined  by  a  different  model,  i.e.,  for  still  air. 
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Let  us  examine  the  expressions  in  Table  10.  Since  x  -  ut,  then  for 
t  =  0,  y  =  9.1  m  and  z  =2.73  m.  This  implies  that  the  source  is  defined 
oy  a  radius  of  —4.5  and  height  2.73  m  values  not  inconsistent  with 
observations  (Fig.  1) 

are  represented  by  the  linear  expressions 


In  Ref.  17,  the  standard  deviations,  0  and  o  , 

y  z 


i  x 

y 


i  x 
z 


(14) 


where  i^  ^  are  the  gustiness  intensities,  in  which  i^  is  the  ratio  of  the 
RMS  lateral  speed  fluctuations  to  mean  time  lateral  speed;  similarly 
The  experimental  values  for  i  and  i  are  given  in  Table  10. 

y  z 


TABLE  10 

GUSTINESS  INTENSITIES  FOR  DIFFERENT  STABILITY  CONDITIONS 


Extremely  Unstable 
Moderately  Unstable 
Near  Neutral 
Moderately  Stable 
Extremely  Stable 


i 

v 

-  JL  - 

0.40-0.55 

0.25-0.40 

0.10-0.25 

0.08-0.25 

0.03-0.25 


i 

2 

0.15-0.55 

0.10-0.15 

0.05-0.08 

0.03-0.07 

0.00-0.03 


We  note  that  the  constants  in  Table  10  fall  within  the  range  given  in 
Table  9  and,  therefore,  since  we  use  the  mean  speed  to  define  x,  we 
assert  that  the  factors  x  and  y  and  z  of  Table  9  are  the  mean  values,  and 
the  derived  volume  and  concentration  are  the  mean  values  of  these 
parameters. 

Above  in  defining  the  Pasquill  stability  categories  and  other  con¬ 
cepts,  we  have  used  mean  values.  Clearly,  fluctuations  of  v'nd  components 
play  an  important  role  in  any  spatial-temporal  history;  e.g. ,  the 
horizontal  fluctuations  will  affect  the  concentration  and  plume  rise.  The 
ability  to  relate  the  fluctuations  to  the  Pasquill  category  would  provide 
a  potential  means  to  find  the  bounds  of  the  mean  value.  However,  the 
data  are  limited.  We  present  some  data  in  Table  11  from  Ref.  2,  where 
a(u)  is  the  standard  deviation  of  u,  the  wind  speed  component  in  the 
direction  of  the  wind,  w,  and  the  vertical  wind  speed,  u  and  u1  are  the 
mean  value  and  fluctuations  of  u,  and  (  u ; ^  ^  /  u  is  called  gustiness. 

We  note  the  following  facts: 

1.  The  standard  deviation  a(u)  decreases  in  magnitude  with 
increasing  stability  since  u  increases  from  stability 
category  A  to  H. 

2.  The  ratio  a(w)/a(u)  i  constant. 


TABLE  11 


WIND  FLUCTUATION  PARAMETERS 
FROM  FAST  RESPONSE  DATA 


Stability  Class 

A 

B 

C 

J(u)m/sec 

1.13 

1.20 

1.61 

C (w) /a(u) 

0.38 

0.40 

0.37 

r2  /- 

V  U  /  u  @ 

1 

3  1.5 

m 

0.214 

0.204 

0.220 

<3 

Z 

=  3.0 

m 

0.198 

0.184 

0.200 

V  w  /  U  (3 

Z 

=  1.5 

m 

0.071 

0.074 

0.085 

■a 

Z 

-  3.0 

m 

0.069 

0.068 

0.065 

D 

E 

F 

G 

u 

1.33 

1.03 

0.89 

0.51 

0.30 

0.36 

0.41 

0.41 

0.45 

0.32 

0.219 

0.199 

0. 193 

0.163 

0.155 

0.198 

0.179 

0.174 

0.136 

0.106 

0.079 

0.079 

0.080 

0.061 

0.035 

0.065 

0.067 

0.065 

0.06 

0.029 

3.  The  horizontal  gustiness  decreases  with  altitude. 

For  our  purposes  we  may  utilize  the  gustiness  value  for  ' u  component  at 
~ 3  m  to  establish  limits. 


From  observations  of  the  time  history  of  the  concentration  at  a  point 
we  find  that  the  concentration  varies  from  near  zero  values  to  peak  values 
far  in  excess  of  the  mean;  see  sketch.  These  observations  are  made  per¬ 
pendicular  or  approximately  so  to  u.  This  effect  of  inceraittency  has 


CONCENTRATION 
(mass/vol ) 
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practical  import  s;.nce  it  is  the  period  of  time  there  is  a  near-zero 
concentration  in  direction  of  u  and  this  may  be  as  low  as  351  of  the  time 
even  at  the  center  of  the  cloud.  Although  we  can  specify  that  length  uf 
time  "A'  is  about  20-40  sec,  corresponding  to  the  period  of  gust i ness , 
we  cannot,  at  the  present,  precisely  define  this  effect  in  the  lateral 
(y)  direction.  Therefore,  the  actual  concentration  as  a  function  of 
time  along  the  line  of  sight  is  not  defined  beyond  the  mean  vaiue  ar.d 

A 

standard  deviation.  This  is  a  serious  practical  shortcoming. 

Csanady  (Ref.  17)  defines,  on  the  basis  of  a  theoretical  development,  the 


ratio  of  A'-/N  as  vexp  (j£)  -  L,  where  is  the  standard  deviation  of 
the  concentration,  N.  The  parameter  has  the  values  of  0.2  (stable, 
smooth  terrain),  0.35  (neutral,  smooth  terrain),  and  0.70  (rough).  The 
peak  value  of  the  concentration  can  be  defined  by  the  following  expression: 


^  -  exp  [J  (2.326  -  0.5c  ) 


n- —  a  exp  l:  (2.326  -  0.5-.;  ]  /  ^  exp  Uj) 
•RMS  1  1 


?o;  an  average  condition,  *  0.5,  N  *  2.3  *  N  .  which  appears  to  agree 

1  P 

with  reported  data. 

Before  we  present  the  overall  diffusion  model  we  must  discuss  the 
buoyant  plume. 


2.5.2  The  Buoyant  Plume 

We  consider  here  the  buoyant  piuse  which  is  formed  as  a  result  of 
exe-thermie  chemical  reactions,  i.e.,  oxidation,  deliquesence.  The 


A  series  of  experiments  is  needed  where  samples  and  cransaissececers 
are  placed  at  several  angles  with  u.  One  of  the  parasecers  should  be 
Pasquili  category. 
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■m. 


excess  temperature  reduces  the  density  below  ambient,  resulting  in  a 
rising  buoyant  volume. 

Generally  speaking,  the  sequence  of  events  for  WP  munitions  appears 
to  be  as  follows:  (1)  the  creation  of  a  hoc  volume  due  to  the  detonation 
and  initial  phosphorus  oxidation,  (2)  an  obscuring  cloud  in  the  vertical 
for  scill  air  or  a  bent  cloud  due  to  atmospheric  winds,  and  (3)  the  tur- 
bulenc  diffusion  of  the  cloud  downwind. 

From  limiced  chermal  imagery  the  spread  of  bulk  material  in  the 
earth-plane  does  noc  appear  to  be  extensive  and  probably  depends  on  the 
shell  trajectory.  Hence,  we  assume  the  source  to  be  limited  and  be  ap¬ 
proximated  by  a  quasi-point  source  for  artillery  and  mortars.  When  wicks 
or  wedges  are  used  the  source  is  normally  a  continuous  linear  source.  We 
now  present  the  general  behavior  of  buoyant  plumes. 

To  estimate  the  effect  of  che  initial  vertical  momentum  which  occurs 

20 

within  a  fraction  of  a  second,  we  uciiize  a  result  from  Sriggs:  che 

vertical  momentum  is  less  important  chan  buoyancy  at  a  distance  x^  downwind; 

X1  -  utf  /*(.Hg/T0)  (17) 

where  W  *  initial  vertical  velocity  (m/sec) 

»  initial  temperature  excess 

Tq  a  ambient  ceaperacure  (K) 

g  »  gravitational  acceleration  (a/sec4) 

Assume  u  »  3  m/sec,  w  *  20  a/ sec.  and  IT  -  500  and  substitute  into  Ed.  17. 

4 

We  obtain  for 

'Briggs,  g.  A..  Flume  Rise.  ESSA.  Oak  Ridge.  Tennessee,  1969  (TID-25075) . 
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x1  «  5.20/9.3  (500/300)  =  6.4  m 


or,  a  very  short  distance  from  the  origin.  For  still  or  near  still  air 
(u  <  0.5  m/sec)  the  initial  momentum  should  be  treated  as  part  of  the 
overall  cloud  rise. 


In  Ref.  1,  we  adapted  the  recommended  formulations  of  Briggs 
(Ref.  20)  for  plume  rise;  that  is,  for  unstable  and  neutral  we  have 


.  ..1/3  2/3  17, 

4h  .  " —  ..*■■■  <F1/6) 


(18) 


Ah 


1.6F1/2x2/3 


,2  r 


0.4  +  0.64  —  +  2.2  — 


1  +  0.8  — 

x* 


-2 


(19) 


4  3  -5 

respectively,  where  F  =  buoyancy  flux  (m  /sec  )  =■  3.7  x  10  (cal  r 

0  4 

burn  time),  and  x^  =  transition  value  for  x  =  0.52F  ‘  (English  units). 

It  should  be  noted  that  the  addition  of  F  is  to  correct  for  cloud  top 
instead  of  centroid.  The  appropriate  values  for  the  heat  release  is  given 
in  the  model  structures. 


Venter'  refining  Briggs'  work,  reports  that  Eq.  18  represents  the 
initial  rise  until  x  »  x^  and  Eq.  19  represents  the  trajectory  after  x#, 
where 


x*  »  14F3/^  for  F  <  55  m'Vs3 

0  4  4  3 

xA  *  34F  for  F  >  55  m  /s 


(20) 


■Venter,  G.  P.  N. ,  "A  Comparison  of  Observed  Plume  Trajectories  with 
Those  Predicted  by  Two  Models,”  Atmospheric  Environment,  11,  1977, 
pp.  421-426. 


53 
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.A 

■i 


where  s  =  stability  parameter  and  the  final  rise  is  given  approximately  by 


1/2  ■>/!- 

Ah  »1.6F  '  (3.5xjk)",'7u  (21) 

for  x  >  3.5xif. 

The  constant  (1.6)  appearing  in  the  above  equations  can  be  shown  to 
be  related  to  the  entrainment  factor  3  with  a  value  of  0.6.  3  represents 

the  cloud  growth  rate  with  vertical  rise;  that  is,  3  ~  cloud  radius/Z. 


For  stable  plume  rise,  Ref.  1  shows  chat  Eq.  18  describes  the  rise 
_  1/2 

until  x  =  2.4u/s  .  The  maximum  rise  was  given  as 

Ah  -  2.9  (F/us)1/3  (22) 

max 

where  again  Che  constant  2.9  relates  to  an  entrainment  factor  •  0.5.  Venter 
gives  an  alternative  form  by  Briggs  for  the  complete  rise  as 


Ah  a  2.9  (,F/us)^3[l  -  cos  (xs  *3/u)  1  ll  3 


(23) 


.  5  * 

which  has  a  maximum  at  x  =  ttu/s  and  is  equal  to  Eq.  22. 


The  different  values  for  entrainment  factor  are  not  unusual  since 
the  magnitude  is  one  of  the  major  uncertainties  in  analytical  treatments 

and  the  entrainment  value  is  established  when  the  predicted  results  are 

■vi 

compared  to  experimental  results.  Shwartz  and  Tulin"“  state  chat  2  can 
vary  between  0.3  and  1.0  and  for  their  work  the  values  ranged  between 
0.42  and  0.60. 


From  experimental  results  one  C3n  estimate  for  the  non-buoyant  plume  that 
the  maximum  rise  at  x  *  100,  500,  1000  m  under  neutral  and  stable  con¬ 
ditions  is  60/u,  105/u,  and  150/u,  respectively. 

"~Shwart2,  I.  and  M.  P.  Tulin,  "Chimney  Plumes  in  Neutral  and  Stable 
Surroundings,"  Atmospheric  Environment,  1972,  pp.  19-35.  (An  interest¬ 
ing  paper  on  the  vortex  formation  of  the  plume  in  a  stable  atmosphere.) 
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Equations  (22)  and  (23)  apply  u  >  0,  i.e.,  the  bene  plume  or  cloud, 
or  x  =  ut  >  0.  Whan  u  *  0,  i.e.,  for  still  air,  Morton  ec  al.  (Ref.  14), 
in  their  classic  paper  give  the  maximum  rise  of  a  buoyant  plume  in  the 
following  expressions  for  che  instantaneous  and  maintained  sources.  The 
expressions  are  given  below: 

Instantaneous 


Hma:<  =  <■»)  (24) 

Q  *  energy  release  (joules).  For  example,  a  cloud  produced  by  100  lb  of 
TNT  will  reach  a  height  230  m. 


Maintained 


,0.25 


max 


31 


(1  +  a) 


3/8 


(m) 


(25) 


where  Q  »  energy/unic  time  (kilowatts),  and  n  ■  racio  of  (actual  lapse/ 
adiabatic  lapse). 


In  conjunction  with  the  above,  we  assume  an  entrainment  factor  of 

0.5  so  chat  che  cloud  radius  is  —0.52.  Hence,  we  can  define  che  approximate 

maximum  diameter  of  che  plume;  that  is,  the  diameter  *  H  .  The  diamecer 

max 

at  any  other  height  ~  2z  or  the  plume  half-angle  is  2b5. 


Briggs  indicates  that  for  stable,  calm  conditions  the  maximum  rise 
is  given  accurately  by 


£h 


5.0F0,'S/i 


3/3 


(26) 


where  s  is  the  average  potential  temperature  gradient  traversed  by  the 
cloud.  (Note  che  analogous  forms  of  Eqs.  22  and  2b.) 
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When  a  ground  inversion  exists  then  the  cloud  will  penetrate  cha 

inversion  if  ( lh)  defined  bv  £03.  id  or  26  exceeds  che  height  of  the 
max 

inversion.  When  cha  inversion  layer  is  elevated,  i.e.,  begins  at  Z  0, 

^‘4  0  6 

then  the  plume  will  penetrate  if  the  inversion  top  z  S  4FJ'  /b^'  for 
calm  conditions  (u  =  0)  and  z^  <  2  (F,ub^)  for  u  >  0.  The  parameter 
b  is  defined  by  the  expression 


bi  "  s“Ti'T 


(27) 


where  -T  *  temperature  difference  of  the  inversion  layer 

a. 

T  •  absolute  ambient  temperature  (’!<) 

7 

g  •  gravitational  acceleration  (a/sec") 


Alternatively,  if  the  plume  characteristic  temperature  excess  exceeds  che 
inversion  temperature  2.1^  then  the  plume  will  penetrate;  that  is 


f  \  \ 


plume 


T  •  F  . 

>  -i 


(28) 


Briggs  (Ref.  19)  suggests  chat  c  =  0.5  when  the  plume  is  rising  through 
neuerai  air  (category  D). 


In  Re:.  1,  we  developed  an  expression  for  plume  (2.T)  as  function  of 
concentration  and  relative  humidity  was  developed.  The  result  is  given 
in  Ftp.  12. 


We  now  return  to  the  exothermic  chemical  munition  and  determine  the 
dimensions  of  che  resultant  aerosol  cloud.  The  overall  cloud  rise,  0.  is 
given  by  the  expression 

2*1*2  h  (29.! 

where  0  *  momentum  rise  (see  Table  10) 
rise  due  to  heat  addition 


it 


-Ih  * 


These  parameters  have  been  defined  above.  The  value  of  the  lateral  spread, 
y,  is  given  in  Table  9.  Since  2  >  y  the  aerosol  cloud  shape  is  assumed  to 
be  defined  as  a  quarter  of  an  ellipsoid,  i.e.,  the  volume  for  a  point 
source  is  given  by 

V(c)  *  -xyZ/6  (30) 

2.6  MULTIPLE  AND  ARRAY  SOURCES 

For  a  linear  array  equally  spaced  instantaneous  poinc  sources,  e.g., 
n-munitions  spaced  "d"  neters  apart  where  d  >  y  (c  ■  5)  they  can  be  created 
as  individual  sources;  see  Fig.  13.  From  Fig.  13  we  note  the  overlap 
occurs  at  time  c  *  (d  -  9.1) /Au  where  A  is  Che  constant  given  in  Table  9. 
y  *9.1  +  Ax.  3eyond  this,  the  lateral  path  lengch  for  n  *  2  will  decrease 
oy  the  amount. 

1  »  d  •  (x,  -  •x1)/'xl  (31) 

where  x ^  is  the  value  of  x  at  time  of  intersection  and  x.,  the  value  of  x 
at  che  poinc  of  interest.  For  n  munitions  we  can  write  immediately 

2  *  (n  -  l)  (d)  (x.,  -  x , )  / x ,  (32) 

-  J,  l 

When  che  wind  speed  is  at  an  angle,  3,  appropriate  adjustment  is  required. 

The  overlap  rones  are  significant  for  radiation  transport  (optical 

properties)  laser  propagation  and  some  cases  of  backscatter .  However,  for 

C  •  L  computation  one  may  treat  che  sources  as  additives,  T)  C  L ( . 

i  1  * 

For  a  linear  continuous  source  of  length.  A,  where  Q  grams  of  smoke 
material  of  yield  Y  are  released,  the  quantity  of  ssoke  produced  per  unit 
length  per  sec  is  (Q  •  Y)/U  •  t)  gas/asee. 


The  concentration  to  a  first  approximation  is  defined  by  the  relation 


Q  *  Y 


3d 


C(x.y.3»t) 


(.x.c)  (1  *  k> 


(33) 


FUure  U.  Mulii-ft -.iiian  Scr^aatic 


where  u  =  mean  speed  (m/see) 

Z  =>  smoke  cloud  height  @  "x"  (z  or  Z) 
k  =  incremental  length  of  smoke  width  @  "x" 

For  the  case  when  wind  speed  is  normal  to  the  smoke  line,  k  *  y/l,  and 
when  u  is  at  angle,  9,  to  the  smoke  k  =  y  sin  9/2..  The  non-uniform  spatial 
distributions  in  the  x  and  2  may  have  a  form  similar  to  Eqs.  3.18  and  3.19. 

For  a  distributed  source,  a  combination  of  point  and/or  linear  arrays 

may  be  applied  with  appropriate  accounting  of  the  time  coordinates;  that 

is,  if  we  have  linear  arrays  at  x_^,  and  x^  then  any  x  is  (x  -  x^) , 

(x  -  x  )  and  (x  -  x,  )  and  the  relative  time  is  t  =  (x  -  x  )/u,  etc. 

j  K  1  1 

A  program  for  distributed  steady  sources  like  smoke  pots  and  generators 
is  given  in  Appendix  D,  Ref.  1. 

When  a  moving  ground  generator  is  utilized  a  vectorial  representation 
Is  applied;  that  is,  the  generator  velocity,  v,  is  added  vectorially  to  the 
wind  velocity,  u,  hence  for  a  headwind,  a  line  source  is  produced  with  a 
resultant  mean  wind  speed  of  u  +  v.  The  value  of  x  is  found  by 


applying  Che  expression,  x  =  [|u|  +  |vj]t.  The  expansion  in  y  and  z 
directions  are  defined  by  Table  9. 

When  the  generator  moves  at  an  angle  to  the  mean  wind  then  the 
resultant  speed  would  be  found  again  by  vector  addition. 

2.7  VALIDATION  TESTS 

In  Ref.  1  we  present  numerous  examples  of  the  model's  capability  to 
predict  a  priori  the  mean  concentration,  C,  and  concentration  x  path  length 
(C  •  L) .  Here  we  will  again  illustrate  the  validity  of  the  proposed 
smoke-obscurant  aerosol  model  and  discuss  any  deviations. 

1.  Test  DP1-002-T-4  (HC  Smoke) 

The  initial  conditions  provided  by  the  test  organization  are  as 
follows : 


u,  mean  wind  speed  -  4.0  m/sec 
Pasquill  Category  -  C 
Relative  Humidity  -  36% 

Munitions  -  36,  M84A.1  canisters 

Array  -  linear 

With  these  data  the  resultant  calculated  and  measured  values  are  shown 
in  Fig.  14.  Although  the  general  agreement  is  good,  better  agreement  could 
have  been  achieved  if  the  non-linear  burning  time  had  been  applied, 
particularly  in  the  60-100  sec  time  period. 

2.  Trial  2  (DP1-002)  (HC  Smoke) 

The  initial  data  provided  by  the  test  organization  are  as  follows: 

u  -8.8  m/sec 

Pasquill  Category  -  D 

Relative  Humidity  -  25% 

Munitions  -  8,  M84A1  canisters 

With  these  data  the  resultant  calculated  values  were  compared  to  both  the 
aerosol  Dhotometev  results  (see  Fig.  15)  and  the  C  •  L  derived  by  the 
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Figure  14.  Trial  4  (DP1-002) ,  19  November  1977,  12:58:00,  HC  Smoke 

( -  C  •  L  Values  Computed  from  Aerosol  Photometers; 

•  Model  Calculations) 
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adjacent  transmissometers  (see  Fig.  16).  First  we  note  the  difference  in 
the  two  C  •  L  plots,  particularly  at  times  greater  than  60  sec.  Second, 
the  absolute  magnitudes  differ.  Hence,  correlation  would  be  difficult 
between  the  two  results  except  for  the  fact  that  they  are  not  collocated; 
that  is,  the  transmissometers  are  displaced  10  m  and  this  can  account  for 
differences  as  noted  in  Sec.  2.3. 

The  agreement  of  the  calculated  points  is  better  for  the  trans- 
missometer  with  the  following  exceptions:  (1)  at  early  times  the  non¬ 
linear  burning  would  reduce  the  first  point  to  about  the  correct  value  and 
also  raise  the  second;  and  (2)  the  observed  dip  at  about  80  sec  is  due  to 
atmospheric  fluctuations. 

3.  Trial  17  (DP1-002)  -  WP  Smoke 

The  initial  conditions  are  as  follows: 

u  -  3.7  m/sec 

Pasquill  Category  -  C 

Relative  Humidity  -  36% 

Munition  -  1,  81-mm  mortar 

The  comparative  results  are  shown  in  Fig.  17  where  the  photometer  C  •  L 
curve  is  used.  The  initial  difference  between  the  two  is  probably  due  to 
atmospheric  fluctuations. 

In  Fig.  18  we  have  plotted  the  C  ♦  L  values  (the  outside  one  applies 
to  "1")  for  the  three  transmissometer  lines  which  are  spaced  60  m  apart, 
i.e.,  r60  m  on  either  side  of  the  one  near  the  sampler  line.  We  note  that 
the  single  mortar  behaves  more  like  a  puff  than  a  quasi-point  source, 
i.e.,  like  multi-mortar  arrays  or  artillery  shells. 

Lastly,  the  small  C  •  L  values  or  the  "flare-up"  of  (3)  at  late 
times  are  not  afterburning  but  concentration  fluctuations  due  to  atmospheri 
instabilities. 
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TRIAL  n  (DPI-002) 

Date:  27  October  1977 
Smoke :  HC 

8  Function  Time  12:24:0:0 
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Figure  16.  C  •  L  Values  Versus  Time  for  Row  0  ( -  Calculated 

Using  Transmittance  and  Extinction  Coefficient, 

• - •  Model  Results) 


Trial  10,  D? 1-002;  WP  Smoke 


Input  data: 

u  -  7.9  m/sec 

Pasquill  Category  -  D 
Relative  Humidity  -  52% 

Munition  -  4,  4.2"  mortars 

The  comparative  results  are  shown  in  Fig.  19,  and  the  agreement  is  excellent. 

5.  Trial  21,  DPl-002,  WP  Smoke 
Input  data: 

u  ~  4.6  m/sec 

Pasquill  Category  -  C 
Relative  Humidity  -  25% 

Munition  -  1,  155-mm  shell 

The  comparative  resulcs  are  shown  in  Fig.  20,  and  the  agreement  is 
excellent.  Here  the  effect  of  fluctuation  at  different  values  are  evident 
and  the  mouel  prediccs  adequately  the  tail  particularly  when  the  fluctua¬ 
tions  are  considered. 


Figure  19.  C  •  L  Tine  History;  Trial  10,  DP1-002,  4-42" 
Mortars  (•  Model  Predictions) 
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Figure  20.  C  •  L  Time  History;  Trial  DP1-002,  1-155  mm 
WP  (•  Model  Predictions) 
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3  TRl,  v_T3,RiL 


MODEL 

The  models  described  here  are  written  in  a  self-contained  format  for 
each  ty pe  of  smoke  and  obscuring  aerosol.  Hence,  there  ia  considerable 
repetition  o:  expressions  that  could  be  compiled  in  several  subroutines 
a  general  computer  program  were  generated. 


3.1  COMPUTATION  SCHEME  FOR  CLOUD  CONCENTRATION  OF  SMOKES  OBSCURING  AEROSOLS 
In  this  section,  a  procedure  is  outlined  to  permit  the  determination 
of  the  concentration,  C  (g/m^) ,  path  Length  in  the  cloud,  L  Cm),  and,  the 

i 

Important  parameter,  C  •  L  (g/m“).  The  units  used  here  are  defined;  however, 
we  utilize  primariiy  grams,  meters,  and  seconds.  «e  initiate  the  discussion 
with  the  requisite  input  data  and  then  proceed  to  the  procedure  proper. 


3.1.1  Inpuc  Data 

A.  Scenario 


Local  target  and  observer  locations 
Munition  deployment  plan 
Munition  race 

Sensor  type  and  operational  wavelength  band 


3.  Smoke  Munitions 


Type 

i.Mk  No . , 

Caliber. 

<2  CO  .  1 

."ill 

weight  ( 

d) 

Fill 

material 

CVP.  R? , 

HC>  Oil*  etc.) 

Burn 

£  iac  and 

rate 

C .  Meteorological  Conditions 
Time  of  day 
Cloud  cover 


}  " 


establish  stability  category 


l 


Mean  wind  speed  (m/sec)  -  surface  to  10  a 

Temperature,  T  (®C)  -  §  10  a  ar.d  0.5  a  above  ground  level 

Wind  direction 

Relative  humidity  (RH  T) 

Elevated  Inversion;  specify  height  above  ground  level  and 
temperature  gradient 

D.  Optical  Properties 

Scattering,  absorption,  and  extinction  coefficients 
Visibility,  if  Eq.  28  and  Tables  9  and  10  of  Ref.  1  are  used. 

3,1.2  Determine  Concentration  for  Sin2le  or  Seoarate  Sulk  V?  Munitions 


u  "  0 

C  (x,y,2,t)  »  CW  •  c  •  Y)  ?  fx  •  y  •  Z/h) 

u  »  0 

C  (x.y.Z.t)  *  (W  •  "  •  Y)  f  <  tZ  3 / 1 2 ) 

(Note:  Cloud  shape  is  a  cone  with  apex  -“8  2  below  the 
ground  level.) 

where  *'  »  fill  weight  (gj 
“  »  fill  efficiency 
Y  *  yield 

x  *  length  in  direction  of  Z  (a) 
y  ■»  lateral  spread  —  to  x  (a) 

2  »  buoyant  rise  (ih)  *  sotasntua  rise  (5)  -  a 
u  *  mean  time  wind  speed  fa/suo 

3 . 1 . 2 . 1  Ha a a  Calculation 
H  *  ii  •  *  •  X 


%  9 

/  • 


1 

V  *  fill  weight  is  specif. for  given  munition 
-  *  1.0  for  ordinary  terrain  surfaces 

*  0.35  tor  soft,  plowed  field  surfaces 
-  0.50  for  muddy  terrain 

Y  -  2.9  *  0.003  CRH  -  iO*)1’*' 

3 .  1.2.3  Calculation  for  y,  Z,  and  Ih;  u  •  0 
f'a)  Establish  Stability  Category 

It  -  T  ( 10  m)  -  T  <  0 . 5 ) 

0  -  IT  +  0.093 

z 

if  -i  •'  0,  unstable 
2 

»  0,  neutral 


■  -T-  -  <V/?  a" j ufa&w-tz*. 


1 


(c)  Cloud  Dimensions 


1.  x  =  ut 

2.  Given  Pasquill  Category,  find  y  in  table  below. 


9.1  +•  G.419x 

Pasquill  Category 

9.1  +  0.328x 

Pasquill  Category 

9.1  +  0. 238x 

Pasquill  Category 

9.1  +  0. 20x 

Pasquill  Category 

9.1  +  0. 18x 

Pasquill  Category 

9.1  +  0. 146x 

Pasquill  Category 

9.1  +  0. 

Pasquill  Category 

3.  Given  Pasquill  Category,  find  z  in  table  below. 

Note:  This  z  value  represents  momentum  exchange;  see  use  when 
Z  =  Ah  +  z  is  presented. 


2.73 

+ 

0.39x  / 

1  - 

u  =  0. 5 

m/sec 

Cat 

- 

A 

2.73 

+ 

0. 25x  ' 

l  =  1.0 

m/s  ec 

Cat 

- 

A 

2.73 

+ 

0.18x  ( 

'  =1.5 

ra/sec 

Cat 

- 

A 

2.73 

+ 

°.15x 

)  =  2.0 

m/ s  ec 

Cat 

- 

A 

2.73 

+ 

0.137x 

2.0 

m/sec 

Cat 

- 

A 

2.73 

+ 

O.llx 

Cat 

- 

B 

2.73 

+ 

0.073x 

Cat 

- 

C 

2.73 

+ 

Q.06x 

Cat 

- 

D 

2.73 

+ 

0.055:: 

Cat 

- 

E 

2.73 

•r 

0.046x 

Cat 

- 

F 

2.73 

+ 

0.03x 

Cat 

•• 

G 
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4.  Buoyant  Plume  Rise,  Ah 


-5 

(a)  Compute  F  =  3.7  *  10  [800  *  n  •  W  f  bumtime] 

j  ^  3 

j  (m  / sec )  using  values  from  Sec.  3. 1.2. 2  and  a 

burntime  of  1  sec  for  artillery  and  1.5  sec  for  mortars. 

i 

J 

(b)  For  Unstable  and  Neutral,  CAT  -  A,  B,  C,  and  D, 

‘  1/2  2/3 

!  Ah(t)  =  1.6F  '  x  f  u,  CAT  A,  B,  and  C 

:  1/2  2/3  2 

=  1.6F' x '  [0.4  +  0.64x/x*  +  2.2(x/x*)  ]  ^ 

■  u  [1  +  0.8CX/X*)]2 

j  where  x^  *  14F3^  F  <  55  mVsec3 

,,.-,0.4  c-xcc  4,  3 

=  34F  F  >  55  m  /sec 

! 

] 

Maximum  Buoyant  Rise  (Cloud  Top) 

;  (Ah)  =  1. 6F1/2  (3.5x*)2/3  t  u 

[  max  * 


n 

I 


Alternative  Rule 

3 

(Ah)  @  x  when  concentration » 0. lg/m  or 
max  _ 

exact  expression  for  AT  in  Fig.  12. 

WARNING:  Check  for  elevated  inversion;  if  none, 

Z  =  Ah  +  z,  otherwise  follow  the  steps  below. 

1.  Determine  height  of  inversion  and  AT  *  T(Z2)  -  T(Z^) 

2.  Calculate  bi  =  gAT^/T  (see  Eq.  27,  Part  1) 

3.  If  inversion  top  defined  by  Z^,  then  plume  pene¬ 
trates  if  ZL  <  2  (F/ub  ,)0,5 

4.  Clearly  if  (AT)  >  (AT^)  the  cloud  will  pene¬ 

trate 
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(c)  For  Stable  Conditions:  CAT  E,  F,  etc. 

Compute  s  =  (g/T)  9  (see  Sec.  3. 1.2. 3) 

Determine  =  2.4  u/s^’3 
1/2  2^3 

Ah  =  1.6F  x  v  u  for  x  <  x^ 

(Ah)  =2.9  (F/us)1/3 
max 

Note:  If  Ah  >  H.,  ground  inversion  height,  plume  will  pene- 
max  l 

trate. 

(Ah)  £  IL  ,  then  is  maximum 

Z  =  z  +  Ah 

The  substitutionof  x,  y,  and  Z  into  volume  expression  gives 
the  concentration  when  u  >  0. 

3. 1 . 2 . 4  Still  Air,  u  =  0 

Concentration  =  (W  •  n  *  Y)  -r  Volume 
(W  •  n  •  Y)  defined  in  Sec.  3. 1.2.1. 

Z  =  (Ah)  -  5.0F°'2Vs3/8 
max 

F,  s  defined  above 

Volume  =  (tt/3)Z3  •  B2,  8  =  0.5 

if  inversion  top  <  Z,  the  cloud  will  penetrate  or, 

.  .-0.4  ,  -0.6 
zi  <  4F  bi 

Note:  Plume  is  a  cone  with  apex  —8  m  below  ground  level. 
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2  2  1/2 

6  =  r/z  =  (XZ  +  y  ;  f  Z 


C(Zmax>  =  (4Z/§) 


1/2 


3.2  METHOD  FOR  WR/RP  -  WICKS,  WEDGES,  PLASTICIZED  MATERIAL 


Municions  Characteristics  Bumtime  (min) 

2.75"  -  Wick,  Wedge  4.3 

3"  -  Wick,  Wedge  7.8 

6"  -  Wick,  Wedge  6.5 

Wavy,  81  mm  -  Wedge  (RP)  4.3 

Navy,  155  mm  -  Wedge  (RP)  6.4 

PWP/PRP  2.0 


We  assume  chat  the  source  is  a  linear -array,  continuous  emitter.  The  material 
is  dispersed  in  a  radius  of  ~25  m,  and  we  treat  the  source  as  three  lines 
25  m  apart  where  the  center  line  contains  50%  munitions  weight  and  the  other 
two  25%  each. 


1.  Basic  Equation  for  Concentration 


C(x,  y,  Z,  t) 


(Q  *  Y  •  n)/(t 


l) 


[12/ttu 


Z(x,  t)  (1  +  k)  ] 


(g/m3 


2.  For  each  line,  Eq.  1  applies;  however,  we  have  at  each  value  of  t 
the  values  of  x  are  (ut  c  25).  Hence,  at  the  position  x  the  cloud 
arrives  at  times,  t  =>  x/u  t  25/u. 


where 


Q  *  fill  weight 
Y  *  yield  for  WP/RP 
H  °  efficiency  factor  of  fill 
t  =  burntime  (sec) 
l  =  length  of  array  (m) 
k  “  y(x)/£  when  u  1  ) 
u  ■  mean  wind  speed  (m) 
y  ■  mean  cloud  width  at  x 
Z  -  cloud  height,  z  +  Ah  at  x 
x  *  u  .  t 
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If  u  is  at  an  angle  9  Co  l,  then  k  =  y(x)  sin  d/l 

Y  =  2.9  +  0.003  (RH  %  -  hJ%)1,67 

where  n  =  0.5  for  muddy  terrain 

0  =  0.85  for  soft,  plowed  terrain 
0=1  for  general  conditions 
l  =  determine  from  operational  plan 
t  =  bumtiae  look-up  value 

3.  From  input  data,  determine  Pasquill  Category 


PASQUILL  STABILITY  CATEGORIES 


Surface  Wind 
Speed  ac  10  m 

(m/sec“T) 


Night 


_ Insolation _ 

Strong  Moderate  Slight 


Thinly  Overcast 
or  >4/8 
Low  Cloud 


<3/8  Cloud 


2 

A 

A-B 

B 

— 

— 

2-3 

A-B 

B 

C 

E 

F 

3-5 

3 

B-C 

C 

D 

£ 

5-6 

C 

C-D 

D 

D 

D 

>6 

C 

D 

D 

D 

n 

4.  Expression  for  y 


y  =  9. 1  +  0. 419x 
=  9.1  +  0. 32 8x 
=  9.1  +  0.238x 
=  9.1  +  0 . 20x 
=  9.1  +  0. 18x 
=9.1+  0. 146x 


Pasquill  Category  -  A 
Pasquill  Category  -  B 
Pasquill  Category  -  C 
Pasquill  Category  -  D 
Pasquill  Category  -  E 
Pasquill  Category  -  F 
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5.  Calculate  k  =  y  (x)/Z  and  (1  +  k) 


6.  Procedure  for  Z  =  z(x)  +  Ih 


(a)  For  Cat  A,  B,  and  C 


.  c  1/2  2/3  .  - 
uft  a  1.6F  x  r  u 


=  3.6  *  10~5 [ 800  *  n  * 


[800  •  n  *  Fill  Weight  (g)  1  burntime  (sec)] 


(b)  For  Cat  D 


dh  *  1.6F1/2x2/3  [0.4  +  0.64  x/x*  +  2.2(x/x*)2] 


u[l  +  0.8  x/x.  ] 


2 


(c)  (Ah)  for  Cat  A,  B,  C,  and  D 
max 


1  /°  2  M 

(ih)max=  (3*5x*)U  *  - 


where  x*  =  14  F 


5/8 


.  4  .  2 

F  <  o5  m  /sec 


34  F 


0.4 


4  2 

F  >  55  in  /sec 


The  values  for  z  for  different  categories  are  given  below. 


+  0. 39x 

u  =  0.5  a/ see 

+  0. 25x  ' 

l  =  1.0  m/sec 

+  0 . 18x  < 

'  =  1.5  m/sec 

+  0. 15x 

|  =  2.0  a/sec 

+  0. 137x 

2.0  m/sec 

+  O.llx 

\ 

=  2.73  +  0 . 0  7  3x 
=  2.73  +  Q.OGx 
-  2.73  +  Q.055x 
«  2.73  +  0.046x 


Cat 

-  A 

Cat 

-  A 

Cat 

-  A 

Cat 

-  A 

Cat 

-  A 

Cat 

-  8 

Cat 

-  C 

Cat 

-  D 

Cat 

-  E 

Cat 

-  F 

(See  notes  for  Bulk  UP  on  inversion  procedure.) 


(d)  Cat  £,  F,  and  G 


Compute  s  =  (g/T)5. 


x^  ®  2.4u  t  s 


Ah  *  1.6Fi/4x2//3  f  u  for  x  £  x 


(Ah)  =  2.9(F/us)L/3 
max 


If  (“^)max  >  H,  ground  inversion  height  will  penetrate, 

(Ah)  <  H,  Z  -  H. 
max  max 


(e)  Calculate  z  for  specific  Pasquiil  Category 


=  2.73  +  0. 39x 
»  2.73  +  0.25x 
=2.73+  0.18x 
=  2.73  +  0. 15x 
=2.73+  0. 137x 
=  2. 73  +  O.lix 
=  2.73  +  0. 073x 
=  2.73  +  0. 06x 
=  2.73  +  0.055x 
=  2. 73  +  0.046x 


u  =  0.5  m/sec 
*  1.0  a/sec 
=  1.5  m/'sec 
=  2.0  a/sec 
2.0  a/ sec 


Cat  -  A 
Cat  -  A 
Cat  -  A 
Cat  -  A 
Cat  -  A 
Cat  -  3 
Cat  -  C 
Cat  -  D 
Cat  -  £ 
Cat  -  F 
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7.  Calculate  concentration 

Substitute  values  for  Q,  y,  ),,  t,  u,  2,  and  k,  compute 
C(x,  y,  z,  t)  for  each  line,  and  add  value  at  each  x  or  time, 
t,  as  needed. 

8.  Find  L,  the  path  length  through  the  cloud,  and  then  compute  C  •  L. 

9.  For  u  *•  0,  use  algorithms  described  for  bulk  W?  in  Sec.  3.2. 

3.3  METHOD  FOR  HC  SMOKE  CONCENTRATION  TIMS  HISTORY 

3.3.1  Munitions  Characteristics 

105mm  Ml  Canister;  Fill  Wt  =  1.65  lb;  Bumtime  =  180  sec  ;  M84  shell 

155mm  *M1  Canister;  Fill  Wt  =  5.8  lb  ;  Bumtime  =  120-140  sec;  M11631  (3/shell) 

*M2  Canister;  Fill  Wt  a  3.0  lb  ;  Bumtime  a  60-240  sec;  Ogive  (l/shell) 

Basic  Equation  for  Concentration 

C(x,  y,  z,  t)  ■  — 

t  ♦ 

Q  *  quantity  of  fill  (g) 

Y  =  '?*■.  i-'  t<  HC 
t  *  bumtime  (sec) 

1.  •  length  of  canister  array  (a) 
k  *  y(x)/2  when  u  is  j_  to  l 
u  «  time  mean  wind  speed  (a) 
y  a  cloud  width  increase  at  x 


12 


•  u  •  Z(x,  ;)(1  »  t) 


(g/n3) 


A 

New  fill  weights  are  5.46  and  2.69  lb,  respectively. 
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Note: 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


(1)  We  assume  a  near-linear  array  for  the  artillery  delivered 
canisters;  i.e.,  the  starting  point  (x  =  0,  t  =  0)  will 
be  the  mean  distance. 

(2)  If  u  is  at  9  to  2,  k 


Y  -  1  +  0.051(RHZ  -  5)°*85 

y  =  9.1  +  0 .419x 
=  9.1  +  0 . 328x 
-  9.1  +  0. 238x 
=  9.1  +  0. 20x 
=  9.1  +  O.lSx 
=  9.1  +  0. 146x 

Z  =  2(x)  +  Ah 


y(x)  sin  9/2 


Pasquill  Category  -  A 
Pasquill  Category  -  B 
Pasquill  Category  -  C 
Pasquill  Category  -  D 
Pasquill  Category  -  E 
Pasquill  Category  -  F 


Ah  correction  for  exothermic  effect 
For  Cat  A,  B,  C 


,  iV'V'3  1.6F1/2c2/3 

;h - - -  or  ~ - 

u  u 


F  «  3.7  •  10  8  [500  cal/g  •  Fill  Wt  (g)  /Bumtime  (sec)] 


t  <,  lOsec 


For  Cat  0 


1.6F  x  •> 

Ah  *  -  [0.4  +  0.64  (x/x*)  +  2.2  (x/x*)“]  i 

u 

(1+0.8  (x/x*)]2 

(”h)»M  *  x.)2/i  *  a 

X,  -  14  F5/8  F  <  55  a4/s5 

**  •  34  F<3,i  F  i  55  a4/s3 
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(7)  For  Cat  E,  F,  etc. 


,,  .  ,  1/2  2/3  .  -  ,  ,  .  ,  -  .  1/2 

Ah  =  1.6  F  x  -  u  for  u  <_  x^  (x^  a  2.4u  r  s  ) 


(Ah)  °  2.9  (F/CIi> 1/3 
max 


a  •  |  e 

T  2 


0  3  Potential  Temp  (2  *  10  m,  2  a  0.5  m) 


(8) 


If  (Ah) 

> 

max 

H 

,  grouns 

and , 

if 

<H, 

2 

*  H. 

max 

z(x) 

M 

2.73 

+ 

Q.39x 

m 

2.73 

+ 

Q.25x 

a 

2.73 

+ 

0. 18x 

m 

2.73 

+ 

0.15x 

sa 

2.73 

+ 

0.137x 

• 

2.73 

+ 

O.llx 

a* 

2.73 

+ 

0.073x 

* 

2.73 

X 

0 

0 

■ 

2.73 

+ 

0 

0 

On 

X 

m 

2.73 

0.046x 

-  0.5  m.'sec 
»  1.0  u/sec 
■  1.5  m/sec 
»  2.0  m/sec 
>2.0  m/sec 


Cat  — 
Cat  - 
Cat  - 
Cat  - 
Cat  - 
Cat  - 
Cat  - 
Cat  - 
Cat  - 
Cat  - 


For  u  a  0  and  inversion  layer  conditions,  see  algorithm  in  bulk  VP 
Sec.  3.1. 

3.3.2  Inputs 

I  (0.5  m)  Munition  Type  -  105  sss  or  155  a 

T  (10  a)  Munition  Quantity  -  '/Mi  and  M2 

u  Length  of  Array 

Pssquill  Category  Relative  Humidity 


max 
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A .  Define  0,  Burntlae 


3 .  Compute  Yield  use  Eg.  1  using  Input  relative  humidity. 

C .  Compute  lh 

(1)  Stipulate  Pasquill  Category 

(2)  Calculate  F,  x*,  or  s 

(3)  Unstable  -  use  Eq.  4 
Neutral  -  use  Eq.  6 
Stable  -  use  Eq.  7 

D.  Compile  2(t)  *  z(t)  4-  lh 

For  given  Pasquill  Category  and  u,  select  appropriate  2(c)  from 
Eq.  8  and  appropriate  lh  Eqs.  4,  6,  or  7.  Calculate  Z(c). 

E.  Compile  y(t) 

For  given  Pasquill  Category,  select  appropriate  y(t)  from  Eq,  2. 
Calculate  y (t) . 

F.  Calculate  "k"  and  (1  +  Sc) 


G.  Compute  C(s,  v ,  2,  t) 

Input  A,  3,  U,  and  F,  u,  and  t  as  defined  under  munitions  into  Eq. 

1.  This  is  mean,  homogeneous  concentration. 

H.  To  find  C  •  L 

Multiply  appropriate  C(c)  value  by  line-of-sight  path  length,  L, 
in  meters. 

I .  Cloud  Width  at  x  »  y(x)  »  i 

Note:  The  persistence  of  the  screen  is  terminated  at  the  end  of  burnt ise. 
For  u  •  0,  use  algorithms  -described  for  bulk  li?  (Sec.  3.1). 
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METHOD  FOR  FOG  OIL  GENERATOR 


3.4 


3.4.1  General  Cons Idarac Iona 

t 

The  M3  fog  oil  generator*  is  a  pulse  Jet  that  vaporizes  15  to  50 
gal/hr  of  oil.  The  weight  is  3.36  kg/ gal.  In  addition,  the  engine  consumes 
3  gal  of  gasoLine/hr  and  exhausts  the  combustion  products.  For  the  present, 
we  neglect  the  latter  although  it  may  have  a  near  field  effect  on  Infrared 
transmission. 

There  is  a  thermo-mechanical  effect  in  the  near  field  of  the  generator 
that  is ,  the  oil  vapor  is  hoc  and  leaves  the  nozzle  at  near  local  sonic 
speed.  It  is  estimated  chat  about  5  kcai/sec  is  released  due  to  condensa¬ 
tion  and  cooling.  Also,  the  jet  would  come  to  near  ambient  air  speed  at 
—3m  from  the  nozzle  exit  plane.  Again  these  effects  are  omitted;  however, 
the  source  is  now  located  -3  at  from  the  generator. 

Hence,  the  source  is  defined  as  a  static,  non-exo thermic ,  point- 
emitting  Q  grams  of  oil  droplets  per  second. 

We  assume  the  sooke  cloud  to  be  defined  by  a  semi-cone  resting  on 
the  ground;  chat  is,  the  x,  y  axes  ace  in  the  ground  plane  and  the  z  axis 
is  _  to  tne  'X,  y)  plane. 

3.-.1  Cloud  Dimension 

3.*. 1.1  Cloud  Length 

x  «  uc  (x  is  in  direction  of  u) 

where  u  »  • ime  mean  wind  speed  at  1-m  height 
:  *  time  (see) 


'When  the  air  temperature  is  -*  60'*?.  the  cloud  dissipates  rapidly. 


3. 4. 2. 2  Cloud  Width,  v 


'/(*)  -  9.1  +  0.4i9x 

-  9.1  +  0.323x 

-  9.1  +  0.233x 

-  9.1  r  0.20x 

-  9.1  +  0.18* 
■9.1+  0. 14Sx 

3.4.2. 3  Cloud  Height,  z 

z(x)  -  2.73  +  0. 39x 

-  2.73  +  0 . 25x 

■  2.73  •>  0. 18x 

■  2.73  +  0.13x 

-  2.  73  +  0 . 1 3 7 x 

-  2.73  +  0. llx 

•  2.73  +  Q.073x 
3  2.73  ▼  0.066x 

•  2.73  +  0.055x 

•  2.73 


Pasquill  Category  ■  A 
Pasquill  Category  3  3 
Pasquill  Category  *  C 
Pasquill  Category  3  D 
Pasquill  Category  3  Z 
Pasquill  Category  3  F 


u  ■  0.5  m/sec 
■  1.0  m/sec 
m/sec 
m/sac 
a/  sec 


Category  A  ■  1.3 

■  2.0 
>  2.0 

Category  3 
Category  C 
Category  0 
Category  E 
+  Q.Q46x  Category  F 


3.4.3  Concentration 

To  Had  the  concentration  during  the  calsaion  phase,  ve  assuse  a 
series  or  discrete  instantaneous  sources  where  the  concentration  re r  an 
instantaneous  source  is  given  by 


C  -  (g/a3)  (31) 

rr*  •  f 

where  Q  »  grass  released  in  1  sec  and  s  and  y  as  defined  in  Secs.  3. -*.2. 
and  3. 4. 2. 2  respectively. 

Hie  conversion  rroa  gal/hr  to  g;' sec  is  0.93. 
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Solution  (31)  can  now  be  write 


en  as  a  function  of  tiuae  since  u 


constant , 


C(t)  - 


T-  «  Cl  9.1  i-  (Au)  •  c]' 


M30  1 _ 

u  C ( 9 . 1  +  Au  •  t)* 


A  13  3  *  *•  specific  PisqclU  Category  and  q 

is  in  gal/hr. 


:or  1-sec  intervals,  we  have  a  series  o ; 


equations  s  .ch  as 


. Lm 


u  [  ±  v  ✓ .  1  +  Au ) ' 


CO  .  LIB  _ 1 _ _ 

u  [C(9.l  .  A  •  u  *  C) ' 


LIB  f 


u  I a(9. 1  ♦a  *  u  *  n ' '  ‘ 


where  n  *  burnt lae  in  seeoaus,  t. 


Now  since  x  is  as  isp licit  function  of  t,  we  have  45 
tiae  1  sec),  etc.: 


x  *  u  (.or  *t 


cfxuij  *  14^  "  £’b[— ■  * 

u  i*l  *  Au  •  t)\ 


cuu>i  .14^  E  [-  ;  -1 

^  V  Au  •  t)’| 


2.X  m'A w, 


For  t  =  t,  ,  then  x  =>  ut,  ,  and 
o  b 


C[x(tfe) ] 


Note:  The  above  simple  procedure  is  suggested  for  ease  of  computation. 
Also,  the  concentration  will  fall  rapidly  and,  therefore,  the  concentration 
calculation  will  be  terminated  at  x  values  <<  x  =  ut^. 


For  t  >  t^  ,  it  is  suggested  as  preliminary  estimate  to  find  the 
mass  near  the  origin  and  treat  as  a  new  instantaneous  source,  at  x  =  0, 
retaining  the  initial  concentration  arrival.  Hence, 


(9.1  +  Au)2 


+  C  ( 5 ) 


•n=- 


5(9.1  +  Au5 ) 2 


3.4.4  C  •  L  Parameter 

With  the  C  given  above  and  the  specified  line  of  sight,  find  C  •  L 
(see  Appendix  F  of  CR-231,  Ref.  1).  It  should  be  noted  that,  near  the 
source,  the  cloud  is  highly  non-uniform. 
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PART  2 

EXPLOSIVE  OUST  OBSCURATION  MODEL 


1 


INTRODUCTION 

In  this  part  of  the  report  we  consider  an  explosive-generated  dust 
model.  The  basic  data  are  examined.  The  early  work  chat  pertains  to  the 
subsurface  events  (performed  by  the  PM/Smoke-Obscurants  during  Smoke  Week  I 
and  at  other  field  trials)  is  described.  In  addition,  several  validations 
of  the  live  and  static  firings  performed  at  Smoke  Week  I,  Fort  Sill,  and 
Grafenwoht  II  are  presented. 

As  will  become  evident,  there  are  difficulties  associated  here  with 
the  validations.  The  number  of  single  events  are  few,  and  multi-events 
are  extremely  complex  phenomena.  Also,  the  available  data  on  the  initial 
conditions  are  limited;  that  is,  the  soil  properties  are  not  defined 
adequately,  the  extinction  coefficients  appear  variable,  cloud  geometries 
are  semi-quantitative,  and/or  meteorological  data  are  limited.  Nevertheless, 
the  geometrical  histories  of  the  dust  clouds  may  be  calculated;  however, 

Che  ability  to  calculate  the  transmission  may  be  limited  by  the  afore¬ 
mentioned  shortcomings. 
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2 


BASIC  CONSIDERATIONS 


2.1  THE  SOIL  MEDIUM 

Since  we  are  Interested  in  the  dust  clouds  generated  by  munitions 

exploded  near  the  earth's  surface  (under  or  on),  the  soil  characteristics 

of  the  first  meter  or  so  of  the  earth's  crust  are  important.  This  includes 

the  physical  and  mechanical  characteristics  and  the  constituents  for  d if- 

* 

ferent  degrees  of  saturation  where  the  latter  is  due  to  precipitation. 

To  avoid  some  of  the  difficulties  such  as  mechanical  properties,  we  will 
depend  upon  experimental  cratering  results  whenever  available  for  the 
different  soils. 

Table  1  lists  the  densities  of  several  soil  types.  Note  the  profound 
effects  of  vegetation  and  moisture  on  the  density.  These  variations  affect 
the  mass  of  dust  injected  into  the  atmosphere  as  a  result  of  an  explosion. 
In  the  literature  we  find  an  occasional  loose  statement  that  the  average 
soil  density  is  2.5  g/cc. 


TABLE  1 

SOIL  MATERIAL  DENSITY 


Material 

Density 

Type 

g/cc 

Ib/ic3 

Quarts 

2.67 

167 

Clay 

2 

124.8 

Loose  loam,  dry 

l 

62.4 

Loam,  wet 

1.8 

112.3 

Loam,  sandy 

1.5 

93.6 

Fine  sand,  loose 

0.7 

44 

Coapact  sand,  wet 

1.8 

112 

Light  soil-grass  roots 

0. 1-0.3 

<31 

Sock  aaterials 

-2.4 

130 

* 

Or  water  table  level. 


At  Grafenwohr  the  results  of  soil  analysis  show,  primarily  after  the 
explosion,  that  near  the  surface  (0—5  cm),  the  density  was  0.95  g/cc,  which 
included  "28.?  vegetacion;  Che  moisture  concenc  was  28%.  In  another 
location  the  value  was  —1.13  g/cc,  and  che  moiscure  concent  was  27%. 

In  general,  che  soil  type  up  to  depth  of  —0.55  m  was  loamy  sand  or 
sandy  loam  wich  an  occasional  gravelly  loamy  sand.  The  moiscure  content  in 
Zone  4  was  typically  — 2 IS— 2 72  and  the  specific  gravity  —2.62  ^ 163  ib/ft3). 

Hence,  che  surface  layer  will  have  a  low  density;  however,  the  pre¬ 
dominant  crater  material  will  have  a  near  average  value  o:  specific  gravity 
for  a  wet  sandy  soil. 

The  depth  of  che  vegetated  soil  (sod)  is  important  since  che  density 
is  lower  chan  che  subsoil,  and  the  likelihood  that  che  sod  does  not 
participate  ia  the  airborne  dust  cloud  due  to  che  weight  of  the  clumps. 

An  analysis  of  the  subsurface  soil  at  Fort  Polk,  Louisiana,  a  silt- 
stone,  is  listed  in  Table  2.  Noce  that  the  depth  is  in  feet. 


TABLE  2 


SUBSURFACE  SOIL 

ANALYSIS  OF 

FORT  POLK,  LOUISIANA 

Moisture 

Degree  of 

Wet  Density 

Dry  Density 

Depth  (ft) 

Content  (X) 

Saturation 

(?)  (Ib/ft3) 

(lb/ft3) 

2-3 

21.1 

42.6 

35.7 

70.8 

3. 3-4. 3 

23.0 

50.2 

91.2 

74.4 

4. 5-5 . 3 

20.2 

36.8 

80.0 

68.9 

)  _  \ 

24.2 

84.8 

113.1 

91.2 

3. 3-4.3 

24.0 

30.2 

111.9 

90.0 

4.5-5 . 3 

22.4 

94.7 

121.9 

100.0 

i-3 

23.5 

79.5 

nrs 

91.2 

3. 3-4.3 

30.0 

80.1 

106.9 

81.9 

4. 5-5. 3 

24.9 

70.8 

106.2 

85.0 

2-3 

29.0 

100.0 

127.8 

107.3 

3. 3-4.3 

26.4 

96.5 

126.2 

105.0 

4.5-5. 3 

25.3 

100. 0 

106.2 

98.1 

2-3 

29.2 

100.0 

121.9 

94.4 

3. 3-4. 3 

24.0 

100.0 

125.0 

100.6 

4.5-5. 3 

29.2 

69.8 

113.1 

37.6 
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We  note  chat  the  moisture  content  does  not  correlate  with  degree  of 
saturation  or  density  at  different  locales  in  the  same  general  vicinity; 
however,  the  percent  of  saturation  correlates  with  the  wet  density.  The 
dry  density  is  obtained  by  dehydration  of  th.e  wet  (actual)  densicy. 


Alluvial  moist  sandy  soil  (sandstone  subsoil)  in  Colorado  has  a 
density  of  about  120  lb./ft^. 

From  Tables  l  and  2  we  observe  two  Interesting  facts:  (1)  the 
quancity  of  silica,  clay,  and  loam  will  determine  the  soil's  basic  density; 
and  (2)  in  Table  2,  sample  1  has  a  greater  specific  gravity  and  a  lower 
dr/  density.  For  sample  1  there  is  probably  somewhat  more  quartz  then  in 
the  other  samples.  Further,  at  half  the  degree  of  saturation  the  percent 
weight  increase  is  the  same;  that  is,  sample  1  is  less  compact  as  evident 
from  the  low  dry  density.  At  100%  saturation  the  wet  weight  would  be 

3 

comparable  far  all  samples  m:  120  Ib/ft  . 

The  chemical  composition  of  soils  daper.cs  on  location  and  the  nature 
of  terrain  (wooded,  desert,  etc.'.  In  Ref.  1  data  for  several  locales  are 
given  and  are  reproduced  here  in  Table  3.  One  must  recognize  that  soil 
composition  is  highly  locale  dependent;  that  is,  to  specify  West  Germany 
is  not  to  imply  that  it  is  applicable  to  ail  of  West  Germany!  Further,  the 
existence  of  only  5%  quartz  in  che  Sinai  is  highly  questionable.  Samples 
from  the  periphery  of  the  Sahara  and  other  deserts  were  show  to  have 
quartz  ss  their  major  mineral  constituent. 

Flanigan  and  Delong’  round  from  their  sampling  effort  at  several 
European  locations  chat  silica  (quarts)  concentrations  were  high  and  that 
little  or  no  CaCO^  was  found  in  6s  different  locations  around  the  globe. 

The  latter  is  eosaon  where  shellfish  are  or  were  present. 


4Gs6U.  G.,  Obscuration  Due  to  Dust  of  a  laser  Sea a  in  a  Gun  Fir  in 
sent.  t‘5A&hC0S.~ Stock  Island  Arsenal.  Report  AC-Ta-75005.' October 

i 

'Flanigan  and  Delong.  Spectral  Absorption  Characteristics  of  the 
Components  of  Dust  Clouds.  EAHUOO.  September  1910. 
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TABLE  3 

CLASSIFICATION  OF  COMPARATIVE  SOIL  BY  WEIGHT 


Mineral 

Chemical 

Aberdeen 

Mo \ ave 

Germany 

Sinai 

Montmor-_ 

iilonice 

(Mg,Ca)0-Al,03 

5SiO,-nH.,Q 

35 

25 

45 

25 

Kao  Ionite 

Ai4S>i4°10(OH)8 

20 

20 

19 

17 

Illite' 

iClAl4SV1l°20(0H)4 

20 

20 

10 

20 

Calcite 

CaC03 

15 

15 

20 

20 

Quarc2 

3iO, 

5 

5 

* 

5 

Hemacice 

F*2°3 

■k 

5 

5 

5 

CGrundum 

Al,°3 

* 

5 

it 

h 

Thenardice 

NaSO , 

* 

) 

* 

£ 

Halite 

NaCi 

* 

■> 

1 

Otners  MaOv  TiQ,  CuO 

MgO,  K.,0.  Na,Q, 

«  r\  *  • 


la  the  aorsal  state  these  clay  minerals  arc  usually  hydrated. 
Met  measured. 

i 


93 


Chemical  analysis  of  dust  aerosols  over  che  Indian  desert  showed  chat 
on  a  percentage  basis  they  fail  into  the  following  groups:  quartz  -  36 ", 
mica  -  20%,  feldspar  -  22",  (includes  the  clay  minerals),  and  carbonates 
(calcite,  dolomite)  -  22".  The  existence  of  the  carbonates  is  expected, 
as  a  high  proportion  of  shellfish  fragments  are  present  in  the  southern 
part  of  the  desert. 

The  important  of  soil  chemistry  is  chat  the  particulate  scattering 
and  absorption  properties,  particularly  in  the  7  to  14  „a  region,  are 
determined  by  che  mineral  type.  Hence,  to  find  the  effective  extinction 
coefficient  integration  of  che  complex  band  structure  is  required.  More 
important,  if  a  laser  transmitter  or  transmissometer  is  used,  significant 
differences  will  be  obtained,  and  che  resultant  may  not  be  representative 
of  che  integrated  8  to  13  cm  bands.  The  latter  may  be  a  problem  in  seme 
published  dust  data. 

In  Fig.  1,  we  show  che  spectral  extiiccion  coefficients,  a  ,  at  3 
to  13  cm  for  Fort  Banning  and  Fort  McLeiian,  respectively.  The  silica 
content  is  low  at  Fort  3enning  and  of  medium  value  at  Fort  McClellan.  The 
effect  ac  9.1?  dm  should  be  noted  and  also  the  difference  in  che  absolute 
value  of  che  extinction  coefficient,  namely,  the  Tort  3enning  values  ire 
2*  Fort  McLeiian.  The  minimum  near  12. 3  am  is  real. 

2.2  PARTICLE  SIZE  DISTRIBUTION 

Anot ho r  parameter  of  basic  Importance  is  particle  size  distribution 
since  it  enters  into  the  determination  of  che  dust’s  optical  properties. 
Normally,  che  distributions  found  in  the  literature  relate  to  surface 
conditions,  whereas  for  near-surface  explosions  we  are  interested  in  the 
subsurface  conditions.  The  former  is  germane  to  ambient  aerosols,  vehicular 
du3t  clouds,  and  airbursts  without  cratering.  For  complete  coverage,  we 
will  consider  bath  surface  and  subsurface  detonations. 

A  standard  soil  claasi: icac ion  chat  is  not  unreasonable  is  given 
in  Table 
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FT.  BENNING 


WAVELENGTH  -  um 


_ _ I _ 1 _ i _ L'V__ . L 

8  9  10  11  12  13 

WAVELENGTH  -  u  T 


Figure  1.  Extinction  Coefficient  for  Ft.  Benning  and  Ft.  McClellan 
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TABLE  4 

STANDARD  SOIL  CLASSIFICATION 


Soil  Class  Sl2e  (urn)  \  by  Weight 


Clay 

5 

22 

Silt 

5-50 

38 

Sand,  very  fine 

50-100 

15 

Sand ,  f ine 

100-250 

14 

Sand,  medium 

250-500 

5 

Sand,  coarse 

50=1-1000 

4 

Sand,  very  coarse 

1000-2000 

1 

Fine  Gravel 

2000-10,000 

1 

In  Table  5,  we  have  listed  the  mass  size  distribution  for  various 
locations  and  soil  types.  Again  the  predominance  of  the  mass  fraction  is  in 
smaller  particles;  i.e.,  for  a  particle  of  similar  density  the  mass  of  50  urn 

3 

particle  is  10  greater  than  for  a  5-um  particle  or  on  a  number  density 
basis  the  distribution  is  predominantly  of  smaller  size  particles. 

Observations  of  ambient  dust  may  show  few  particles  beyond  10  '-a  in 
size,  whereas  explosives  may  generate  dust  clouds  with  several  orders  of 

4 

magnitude  (  —  10  )  of  100- 'am  particles  or  larger.  For  example,  at  Dugvay 
Proving  Ground,  the  surface  soil  mass  median  diameter  (HMD)  was  254  urn, 
whereas  at  3- in  depth  the  subsoil  had  an  hWD  -  354  u».  The  Smoke  Week  I 
tests  held  at  Dugway  Proving  Ground  showed  chat  only  0.331  of  the  particles 
were  greater  chan  75  van.  On  the  other  hand,  at  Fort  Sill.  Oklahoma,  the 
particle  sizes  terminated  at  <,10  ua. 

In  Ref.  3  measurements  of  the  particle  size  lofted  by  a  tank  main 
gun  muzzle  blast  at  a  dry  lake  bed.  China  Lake,  showed  that  particles  of 
radii  «;100  ua  made  up  501  of  the  dust  cloud  mass.  This  aass  distribution 


measurement  of  Laser  geas  Transmission  and  ftlS  Sight  Obscuration  in  a 
Tank  Gun  Firing  Environment,  Fard-Aeroautronic  Division,  California, 
Report  So.  U-8362,  Karsh  1977. 
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is  characteristic  of  the  surface  layer.  A  similar  mean  mass  diameter  is 
found  for  the  DIRT-I  test  area  at  the  White  Sands  Missile  Range;  see  Fig.  2. 
The  additional  data  points  are  results  from  the  Grafenwohr  tests  (November 
1978):  the  dots  are  from  a  sandy  clay  loam;  the  lower  points  are  from  a 
sand  soil;  and  the  squares  are  the  results  for  a  soil  like  the  one  at 
WSMR.  We  observe  that  the  MMD  varies  from  100  ym  to  ~  2000  ym  (coarse 
sand).  Interestingly,  the  two  silty  sands  are  similar  except  for  the 
variation  in  gravel  content;  however,  the  overall  effect  should  be  small. 

In  Table  6,  we  present  the  percentage  of  particles  below  74  ym  at 
several  US  military  installations. 


TABLE  6 

PERCENTAGE  OF  SURFACE  SOIL  PARTICLES  BELOW  74  ym 


Location 


Particle  Percentage 


Ft.  Hood,  Texas 
Ft.  Bliss,  Texas 
Ft.  Bliss,  Texas 

White  Sands  Missile  Range,  New  Mexico 
Ft.  Huachua,  Arizona 
Yuma  Proving  Ground,  Arizona 
Yuma  Proving  Ground,  Arizona 
Aberdeen  Proving  Ground,  Maryland 


16-38 

5-8 

24.4  (Field) 
>8.3 


1.4-9. 4 
9.1-21.0 
<45 

7.6-30.4 


Hence,  we  should  anticipate  large  fractions  of  particles  in  excess  of  75  ym. 
This  statement  therefore  raises  a  question,  "How  valid  are  the  Smoke  Week 
and  Ft.  Sill  sampler  results?"  What  did  the  samplers  actually  intercept? 

Although  100-ton  detonations  are  not  comparable  to  artillery-like 
detonations  due  to  deeper  craters  and  greater  energy  release,  nevertheless 
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Figure  2.  Typical  Grain  Size  Distribution  Curve  for  DIRT-I  Test  Area 
(Graf enwohr  Data  Supplied  by  Waterways  Exp.  Station,  Depart 
ment  of  Array) 


their  particle  distributions  have  qualitative  interest.  Particle  distributions 
following  explosions  of  charges  on  the  ground  are  presented  below. 

The  in  situ  material  (clay  soil;  water  table  4  ft)  prior  to  detonation 
has  a  distribution  given  by 

dS  5.6  x  104 

df - ~  <l) 

r 

It  should  be  noted  that  the  primary  emphasis  was  on  the  larger  particles, 
radii  >10  urn,  and  therefore,  extrapolation  to  smaller  radii  are  required. 

We  assume  the  distribution  is  valid  to  —  1  ym.  At  92  sec  after  the  detona¬ 
tion  the  sampling  of  the  rising  main  cloud  showed  the  distribution  in  Table  7. 


TABLE  7 

PARTICLE  SIZE  DISTRIBUTION  IN  MAIN  CLOUD 
(100  TON  HE,  COLORACC) 


D  (pm) 

dN/d  (log  D) 

26 

0.212 

48 

0.124 

80 

0.0606 

112 

0.035 

160 

0.021 

224 

0.011 

320 

5.5  *  10~3 

448 

2.37  «  10° 

640 

1.19  *  IQ’3 

896 

6.4  x  io“5 

1280 

0 

This  distribution.  like  others,  has  3  fora  given  by  the  expression 


<1S  ^  A 

d  (log  07  *  ”n-l 


(2) 
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where  A  is  a  "number  density,"  or  dN /dr  *  A/r11.  This  form  is  approximately 
Che  same  as  the  surface  layer;  however,  che  exponent  was  3.3  for  particles 
less  chan  500  ya. 


At  a  later  time,  —  150  sec  sampling  of  che  cloud  base  and  stem  cop 

showed  an  enrichment  of  large  particles  £200  urn  due  to  settling;  that  is, 

-2  34 

the  size  distribution  ~  r  '  for  che  particles  <500  pin.  In  che  middle 
of  the  stem  che  slope  is  r  ^  and  is  compatible  with  the  ejected  dis¬ 
tribution.  Shortly  thereafter,  we  established  from  the  data  a  particle  size 
distribution  function  near  the  cloud  base/stem: 


dN  7.75  x  103 
dr  “  r3.69 


The  major  difference  is  a  variation  in  the  number  of  particles  >10G  urn 
and  their  iocacion.  For  particle  sizes  <60  yn  che  numbers  remain,  approxi¬ 
mately  ;.he  same. 

This  sampling  cook  place  during  high  winds.  When  the  winds  were  mild, 
che  efrect  on  the  large  particles  results  in  earlier  fallout. 

In  another  test,  che  oartic.e  distribution  near  ground  level  was  found 

to  be 

dN  8.5  *  10* 

37  *  — TT”  (4) 

r 

With  respect  to  particle  size  distribution  from  artillery  and  smaller 
changes,  two  particle  distributions  observed  in  DIRT-1,  under  mild  wind 
conditions,  have  been  reported  in  Ref.  4.  The  first  was  obtained  at  120  sec 
after  a  barrage  of  8-155  s  projectiles  where  che  sampler  was  at  10.6  8  above 

’'lindberg,  J.  S..  A  Preliminary  Report  on  Dusty  Infrared  Teat  -  l  (DIRT-1) , 
A3L-TS-0021.  January  1979. 
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ground  level.  We  note  chac  che  particle  sizes/dlaaecers  varied  from  — 1 
to  100  -o.  Ac  this  time  the  cloud  should  be  ac  —  40-50  m  alcicude  and 
retain  most  of  che  lofced  particles  (no  ballistic  trajectories  assumed) ; 
however,  che  large  particles  (>,100  um)  should  have  settled  to  the  —  10-m 
level.  The  second  event  consisted  of  15-lb  bare  charges  totaling  140, 
deconaced  in  a  rectilinear  array.  Here  the  measurement  was  made  ac  150  sec 
after  detonation  and  ac  an  alcicude  of  —26.3  m  above  ground  level.  The 
number  of  small  and  large  particles  were  both  larger  chan  in  the  above  case. 
The  greater  number  of  small  particles  is  understandable  in  view  of  che  later 
time  and  sampling  alcicude.  With  respect  to  the  larger  number  of  large 
particles,  this  may  be  due  to  che  charge  array  and  che  interacting  efface  of 
charges  which  may  have  lofted  a  larger  fraction  of  loose  surface  dust  to  a 
higher  alcicude  chan  normal. 

A  fit  to  che  data  gives  the  expression. 


dN  1.79  <  IQ& 
dD  *  Q2.65  ; 


dicer"1',  uo”1) 


(5) 


or  in  units  of  (cm  \  um 


1 


(6) 


The  exponent  (2.65)  is  characteristic  of  large  particle  enrichment  as  noted 
earlier. 


We  have  assumed  in  our  analysis  that  the  particle  sise  distribution 
is  given  by 


4 £  m  8,5  <  10S 

dr  * 


(7) 


which  is  tare  characteristic  of  standard  aerosol  distributions  and  observed 
functions  of  airborne  dust  generated  by  explosives.  Although  there  is  a 


difference  between  che  two,  che  basic  quescion  is:  What  is  che  cime-space- 
dependenc  particle  si2e  distribution?  The  optical  properties  will  depend 
upon  this  question,  particularly  che  particles  wich  radii  i  5  pm. 

We  will  attempt  to  formulate  the  parcicle  size  distribution  history. 

First,  we  must  establish  che  distribution  of  the  soil  type  that  will 
be  lofted  into  che  acmosphere.  For  this  we  assume  that  che  distribution 
is  given  by 


dN  _A _ 

dr  *  3.6 
r 


(8) 


Obviously,  che  exponent  will  vary  between  —3.3  and  3.7.  Mext  we  assume  che 
important  particle  radii  lie  between  —1  and  200  pm.  Now,  che  mass  density 
(g/cc)  of  this  function  is 


m  • 


(9) 


where  o  *  density  of  material  (g/cc). 

If  we  assume  p  ®  2  then  M  can  be  obtained  directly;  that  is 

a(t)  *  76.7  <  1Q~12A(c)  (10) 

\  /  3 

To  find  "A"  we  consider  a  time,  c’,  »  Q.ivr  ,  where  U  ®  lb  of  HE.  Now  we 
define  a  cratering  efficiency,  n(cc/lb  HE).  Then  the  ejected  mass  for  W  lb 
of  HE  is 


M  •  n  •  u  •  p  (grams)  (11) 

where  p  is  the  sail  density.  Since  the  entire  mass  is  not  lofted  we 
introduce  a  fractional  fallout  mass,  s»q,  to  represent  the  "instantaneous 
fallout."  Hence,  che  airborne  fraction  is 


MU')  •  (1  -  Oq)  •  n  •  u  •  o 
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(12) 


n^rau 


The  factors  Oq  and  o  are  functions  of  soil  cvpe  and  soil  wetness  and  will 
be  given  later.  Thus  the  factor  A  at  t*  for  a  cloud  volume,  V,  is  defined  by 


(1  -  mn)  -!  •  W  •  c 

A(t’)  -  - - - p -  (13) 

76.7  x  lQ~l~V(c') 


If  che  concentration  were  uniform  than  A(t)  would  decrease  as  V(t')/V(t) 
ar.d  similarly  m(t)  .  However,  there  is  settling  of  the  larger  particles 
and,  therefore,  the  vertical  particle  size  distribution  changes  with  time 
as  well  as  the  concentration.  The  lacter  depends  upon  wind  speed  and 
particle  size. 

To  simplify  the  problem  we  note  that  electro-optical  transmission  is 
limited  to  less  than  60  sec,  probably  ~30  sec.  Within  this  time  period 
the  upper  part  of  the  cloud  should  be  free  of  particles  with  radii  >200  urn; 
however,  these  should  have  settled  near  the  base  of  the  main  cloud.  Further 
particles  from  the  stem  will  be  blown  into  the  cloud  base;  however,  only 
10%  of  the  entire  lofted  mass  is  in  the  stem. 


Hence,  we  assume  that  the  particle  size  distributions  are  given  by 
the  following  expressions,  when  a  >0.3  m/sec,  u  <  10  m/sec.* 

For  c  <  5  sec 


dN  ^  A 
dr  "  p3.6 


(14a) 


For  5  <  t  S.  15  3ec 
The  cloud  base  23  m 


dN  A 

«  73 


(14b) 


See  Part  i  for  definition  of  u. 
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The  remainder,  cloud  and  seem 


dN  A 

dr  “  3.7 

r 


For  c  >  15  sec 


dN  A 

dr  *  3.7 

r 


\ 

(14c) 


(14d) 


The  values  of  A  are  found  by  the  procedure  described  above  with  the 
auxiliary  condition  that  the  total  mass  for  t  >  0.3W^3  is  distributed 
accordingly,  90S  in  the  main  cloud  and  10S  in  the  stem. 

Additional  discussion  is  given  later  when  the  concentration 
distribution  is  considered. 


The  above  functions  are  needed  to  establish  the  Hie  values  for  the 
extinction,  scattering,  and  absorption;  chat  is,  the  optical  properties  of 
the  dust.  We  consider  this  topic  in  the  next  section. 

2.3  DUST  OPTICAL  PROPERTIES 

From  the  above  discussion  we  observe  that  to  specify  the  dust 

(explosive)  optical  properties  is,  to  say  the  least,  a  difficult  task. 

In  Table  3  we  have  assembled  several  experimental  values  and  a  theoretical 
* 

one.  The  theoretical  value  for  the  photopic  region  (0.4  to  0.7  ya)  is 
baaed  on  a  water-like  particle.  Me  have  omitted  the  effect  of  diffraction, 
which  will  increase  the  forward  scattering  and  reduce  the  extinction 
coefficient. 


• 

Patterson,  S.  M. ,  "E.ttlnctjon  Between  0.55  ya  and  10. 6  urn  due  to  Sojl-Ber iued 
Aerosols,"  Applied  Optics,  1977,  pp.  24)4-2418.  In  this  reference  Patterson"' 
reviews  and  discusses  his  results  on  the  extinction  coefficient  induced  by 
varying  atmospheric  dust  aerosol  site  distributions.  The  variation  due  ;g 
distribution  is  *-2  orders  of  magnitude.  For  each  spectral  distribution,  the 
photopic  value  is  -l*  the  10-11  „o  value  and  the  variation  within  the  3-5  ya 
region  is  small. 
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T ABLE  3 


SPECTRAL  EXTINCTION  COEFFICIENTS  FOR  OUST"' 

(m2/g) 


Wavelength  (urn) 


Duac  Origin 

. 4- . 74 

.67  3-5 

3.3 

1.06 

3-13 

Source 

Edgewood 

.55  .4 

.30.4 

CSL 

Fort  Sill 

.06 

.04 

.05 

.03 

DPGFR-78-313 

Fort  Sill 
(Surface) 

•  24 

.16 

.19 

.13 

DPGFR-73-313 

Dugway 

.32 

.27 

.35 

.21 

Smoke  Week  I 

Wyoming 

Clay3 

.07- 

.125 

Aerodyne  T.nc. 

Theoretical 

.  4-.  3 

.13 

Eglin  iVFB 

.2? 

.25 

.26 

.23 

Smoke  Week  II 

See  Part  3  for  addle  tonal  data  on  the  values  for  the  photopic  region. 

A  single  major  deviation  is  the  Fore  Sill  photopic  extinction  coefficient; 
that  is.  The  lowest  value  observed  for  field  data. 

^Crushed  to  optimize  size  to  4  urn  particle  diameter,  i.e..  comparable  to 
surface.  Also  a  minimus  should  exist  at  ~  3.3  urn  for  the  3-5  ua  band. 


The  variations  in  Table  8  are  primarily  related  to  composition  and 
particle  sice  distribution  (location}.  To  appreciate  variation  in  absolute 
magnitude  see  also  Fig.  1. 

The  theoretical  calculation  at  the  3  to  13  ;a  region  was  jerforsed 
by  &.  Frickei.  C31.  and  the  input  data  were  provided  by  this  writer  as 
representative  o;  Crafenwohr  (vet  site).  The  inputs  were  mineral  content, 
index  of  refraction  and  particle  distribution.  Our  estimate  for  the  mineral 
content  of  che  soil  is  as  follows;  quarts  (40t);  feldspar  (28S).  mica  (30*). 
and  carbonates  (?i) .  From  this,  th«  complex  index  of  refraction  was  estab¬ 
lished  with  che  effort  focused  on  the  3-13  ^  region;  see  Table  ?.  -ith 
these  values  and  the  particle  sise  distribution  dS/dr  *  %  «  10  /r  ‘  .  the 
mass  extinction  coefficients  were  computed  for  spherical  particles  whose 
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TABLE  9 

WAVELENGTH  VARIATION  OF  REFRACTIVE  INOEX 


Wavelength  (um) 


Real  Index 
1.130 
1.300 
1.400 
1.700 
1.730 


Imaginary  Index 


1.620 

1.470 


upper  and  lower  size  limits  were  varied.  Some  results  are  given  in  Table 
10.  and  a  variation  of  the  coefficient  with  particle  size  is  shown  in 
Fig.  3.  For  the  present  scudy  we  used  a  lower  limit  for  r  ■  2  ua  and  an 
upper  limit  of  200  ua.  From  Fig.  3  we  note  that  the  important  region  for 

this  wavelength  band  is  r  <  10  ym.  Hence,  a  good  representation  of  the 

* 

particle  distribution  is  significant.  For  shorter  wavelengths  the  impor¬ 
tant  region  shifts  inversely  with  wavelength  since  the  particle  size 
parameter,  x,  is  defined  by  2"rA  and  x  <  5  is  the  significant  range. 

In  the  above  discussion  we  did  not  examine  an  effect  of  the  explosion; 
that  is.  the  possible  dehydration  of  some  ejected  mass  {particularly  the 


Although  this  statement  is  self-evident,  we  make  this  apparent  from  the 
analytical  expression  for  the  extinction  coefficient. 


Ua"1)  »  ^Qe<a.  x^  'O  dr  d< 


iXfe  * 


Q  Cm,  s..  dr  dV 
c  i 


where  Q  »  efficiency  factor  extinction.  C  »  cloud  concentration,  and 

S  m  palticle  density.  ’'<a' 

? 


TABLE  10 

SPECTRAL  EXTINCTION  COEFFICIENT 


.  (urn) 

3  (m“/.i 

a. 2 

.059 

8.5 

.089 

8.7 

.101 

9.0 

.128 

9.5 

.133 

10.0 

.140 

11.0 

.142 

13.0 

*-* 

V>> 

00 

3 


a 

w 


Z 

Ui 


8 


s 

M 


**A«riae  oivcrf*  Cum:? 


Figure  3.  SsjiaecieH  Ce^f'r'iii^as  aa  4  Fuae«4«a  if  fzt tieU  Diis^Sc: 

(waveUagch  *  9 . S  us;  1.730,  0.315) 


minerals)  since  the  fireball  will  be  >600*  K  for  1  sec  detonation.  Also 

soil  moisture  is  —  20*  of  the  total  mass  ejected.  The  moisture  should 

* 

appear  either  in  liquid  or  vapor  phase.  A  possibility  is  that  the 
moisture  recondenses  on  the  particles  in  the  dust  cloud  and  therefore, 
the  cloud  particles  are  similar  to  the  ejecta.  For  the  present  we  neglect 
chis  effect;  see  brief  dicussion  in  3ec.  3.3.1. 

Also,  there  is  produced  abouc  0.2  g  of  carbon  per  gram  of  TNT.  This 
aerosol  should  be  included  with  the  dust.  A  preliminary  estimate  of  the 
affect  would  indicate  that  the  transmission  can  vary  by  —25X. 

Thus,  we  see  that  the  optical  properties  are  varied,  and  a  quantita¬ 
tive  understanding  of  the  soil  type,  etc.,  is  required  to  obtain  a  reason¬ 
able  estimate  of  the  transmittance. 

2.4  EXPLOSIVE  EFFECTS 

The  effects  of  explosive  type  can  have  some  importance  either  for 
charge  weight  equivalent  or  possibly  cratering  phenomena.  We  follow 
convention  and  use  TNT  as  the  reference  explosive  in  Table  11.  The  explosive 


TABLE  11 

EXPLOSIVE  TYPES  AND  PROPERTIES 

cal/g  TNT  Equivalent  Detonation  Velocity  (as/usee) 
TNT  934  1.0  6.93 

Cos?  3  1119  1.19  7.99 

(51DX/TNT/UAX)  1033  1.19 

Coop  C4  1200  1.1  3.04 

Explosive  0  761  0.36 


This  is  not  true  of  the  Fort  Sill  dess. 

*  3 

Consider  the  crater  efficiency  of  1  ft  /lb  then  the  lofted  material 
-10  ib/lb  of  HI  or,  SO*  the  amount  of  carbon.  If  the  mass  extinction 
iot  carbon  is  10*  dust,  shea  the  contribution  to  attenuation  is  -  AQI. 
The  aseuat-  of  carbon  will  depend  upon  the  type  of  explosive  used! 


Corap  3  is  also  an  RDX/TUT  macoriai.  the  racio  of  RDX/TNT  differs  from  che 
mix cure  lisced  diraccly  below.  With  respect  to  cratering  a  lower  detona¬ 
tion  velocity  explosive  will  be  more  effective.  We  assume  the  effect  is 
linear  with  velocity. 

Above  we  used  the  value  of  984  cal/g  for  TNT,  a  quantity  whose 
exact  value  is  uncertain.  Some  investigators  give  a  value  of  —740  cal/g. 
Experiments  performed  on  encased  charges  indicate  chat  the  effect  of  total 
shell  weight  should  be  included.  The  expression  established  experimentally 
is  given  as 


Heat  Released  (cal/g)  •  740  (total  weight /HE  weight)'"  (15) 
Thus,  for  a  155-ram  shell  ve  would  have  - 1000  cal/g. 

Since  che  high-explosive  filler  is  not  a  sphere,  che  lengch/diaraeter 
racio  of  the  filler  and  Initiator  location  should  produce  strong  effects  on 
cratering,  etc.  These  effects  are  omitted  at  this  time;  however,  they 
should  be  examined  In  an  upgraded  model.  The  reader  is  referred  to  Ref.  5. 

To  account  for  ambient  conditions,  i.e. ,  pressure  altitudes,  we 
propose  che  use  of  the  Sach’s  scaling  lav.  Hence,  che  effective  W  of 
explosive  is 


V  ,  . 
^5 


Jistsjaai) 

^  'altitude  ' 


(lb) 


“or  later  use.  we  include  Fig.  «.  which  gives  the  significant  blast 
parameters  for  an  air  blast  in  sealed  form.  These  include  pressure, 
arrival  tine  of  the  shack  wave,  shock  duration,  and  material  velocity, 
for  a  surface  burst  the  equivalent  charge  weight  is  -!  and  therefore  at 

*  /  J 

a  given  range  the  peak  pressure  will  be  equivalent  ee  a  range/-*  .  for 


isetsWi.  J. .  and  «.  U.  Snyder,  Characteristics  of  blast  Waves  Obtained 
from  Cylindrical  K5  Charges.  Denver  Research  Institute,  DRID-dd.  19eS, 


Hi 


f  I  ILK  <V  ASUtlVdi.  KKU  K>3|TjtE-»"!tA&f  UUK4TIOP  (•••k/ib 


PEAK  OVERPRESSURE  (PS. ) 


Figure  4  (Cont.).  Significant  Blast  Parameters.  Relation  of 

ideal  blast  wave  characteristics  at  the  shock 
front  to  peak  overpressure. 


115 


PEAK  DYNAMIC  PRESSURE  (PSI) 


underground  deconations,  the  equivalent  free  air  charge  weight  is  given 
by 


W(air)  =  W(ground)  exp  (-3oX.)  (17) 

d 

where  p  =  soil  specific  gravity 

Ad  -  charge  depth/W1/3  (ft/lb1/3) 

The  latter  is  valid  for  X  <  2,  the  range  of  interest  here. 
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CLOUD  DEVELOPMENT  MODEL 


Here  we  expand  the  model  presented  in  Ref.  6  to  include  detonations 
occurring  at  various  depths  and  heights  relative  to  the  surface,  and  an 
upgrade  of  the  thermal  effects.  Although  we  discuss  the  cratering  effects 
produced  by  spherical  bare  charges,  we  recognize  that  differences  can 
exist  with  artillery  shells  and  ocher  munitions. 

3.1  CRATERING  EFFECTS 

The  explosive  energy  coupled  into  the  ground  to  produce  cratering 
is  a  function  of  burial  depth.  For  a  half  buried  charge,  i.e.,  a  surface 
burst  with  the  charge  center  at  surface,  the  coupling  is  abouc  202  and  will 
rise  to  1002  when  the  depth  (ft)  *  where  W  is  the  charge  weight 

(lb).  For  near-surface  bursts,  i.e.,  above  ground  level,  the  coupling  drops 
rapidly.  Hence,  the  fuzing  effect  can  have  a  profound  effect  on  the 
crater  si2e. 

Although  the  mechanical  interaction  between  the  coupling  energy  and 
terrain  is  significant,  the  topic  is  beyond  the  scope  of  this  study.  For 
our  model,  we  will  relay  upon  experimental  data  to  formulate  the  cratering 
effects  in  different  soils  and  conditions  as  a  function  of  charge  location 
relative  to  the  ground  surface. 

We  will  now  describe  the  results  of  several  experiments  performed 
under  varying  conditions  for  the  express  purpose  of  deriving  a  basis  for 
scaling  laws.  Again,  we  note  chat  we  have  adopted  the  established  procedure 
of  relating  spatial  coordinates  and  times  to  (charge  weight) . 

Since  we  are  primarily  interested  in  near-surface  bursts,  discussion 
will  emphasize  the  range  where  the  charge  center  location  is  :1  charge 
radius  above  or  below  ground  level  (see  Fig.  5)  where  ST  and  HB  are  used  to 
denote  surface  tangent  and  half  buried,  respectively. 


Zirkind,  R.,  A  Preliminary  Description  of  an  Explosive  Dust  Cloud  Model, 
General  Research  Corporation  Report  TM-235,  1978. 
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HOB  -  1  (ST) 


HOB  -  0  (HB) 


DOB  =»  1 


Figure  5.  Definitions  of  Charge  Placement 


For  clarification  we  note  that  a  1000-lb  spherical  TUT  charge  has 
a  radius  of  1.33  ft.  Thus,  at  an  HOB  »  1  (geometrical  height  to  charge 
radius  ratio  *  1)  the  center  of  the  charge  is  located  1.33  ft  above  ground 
level.  Similarly,  for  HOB  »  2,  the  center  of  a  1000-lb  charge  is  located 
2.66  ft  above  the  ground  level.  For  a  1-lb  charge  the  radius  is  0.133  ft. 

In  each  of  these  cases  a  crater  will  be  formed  whose  siae  depends 
upon  cnarge  weight  and  soil  conditions.  In  Fig.  6  we  provide  some  defini¬ 
tions  applicable  to  the  crater  where  "apparent”  means  the  created  physical 
void.  It  is  assumed  here  that  the  volume  of  the  apparent  crater,  V  , 
is  initially  discharged  into  the  atmosphere  and  creates  the  primary  and 
secondary  dust  clouds.  The  latter  may  be  caused  by  stones,  etc.,  (salta¬ 
tion  process)  but  for  the  present  is  neglected.  Our  major  interest  will  be 
on  the  following: 

c  or  r  •  crater  radius,  apparent 
r  a 

c^  or  d^  »  crater  depth,  apparent 

One  should  note  that  the  lip  contains  some  ejecta  and  compact  soil  upthrusced 
by  the  explosive  force.  The  lip  ejecta  is  fall-back  material  and  is  non¬ 
trivial.  Let  us  now  examine  the  test  data  and  derive  therefrom  a  suggested 
set  of  rules. 

We  initiate  the  discussion  with  a  sec  of  data  obtained  at  Fort  Polk, 
Lousiana:  1000-lb  TUT  spheres  were  detonated  in  or  on  silt  stone.  The 
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soil  was  wee  (>  35%  saturation)  except  for  one  event,  >H,  where  the  satura¬ 
tion  was  only  about  50%.  The  pertinent  results  are  given  in  Table  12. 


TABLE  12 

APPARENT  CRATER  DIMENSIONS* 


Event  No. 

1 

2 

3 

4 

5 

H03  (charge  radius) 

-1 

0 

+0.5 

+1 

-0.5 

r  (ft);(r  /W1/3) 

a 

11.4 

(1.14) 

10.4 

(1.04) 

7.9 

(.79) 

6.9 

(.69) 

10.2 

(1.02) 

da(ft);(d/Wi/3) 

5.7 

(.57) 

4.3 

(.43) 

2.7 

(.27) 

3.0 

(.30) 

5.7 

(.57) 

V  (ft3)  -  parabolic 

359.3 

703.1 

254.6 

114.4 

729.4 

Soil  density  (0/ft3) 

85 

112 

108 

120 

120 

Crater  Efficiency 
(ft  3/ con) 

(:c3/lb) 

1719 

.56 

1406 

.703 

509 

.255 

229 

.115 

1459 

.73 

Harvey,  W.  T. .  et  al.,  Near-Surface  Cratering  Experiments,  Fort  Polk, 
Louisiana,  AFVL  TR-74-351.  November  1975. 


ye  ebserve  that  the  crater  raaius  a^u-rs  ?«  increase  with  decreasing 
height  of  burst  siaee  the  coupling  sr.ergy  into  the  ground  increases. 

In  Fig.  7,  we  see  the  effeec  of  some  soil  types  on  the  crater  radius 
and  depth  as  a  function  of  scaled  burial  depth.  Clearly  the  two  excretes, 
narine  muck  and  basalt,  behave  33  expected.  Further,  as  the  burial  depth 
increases  the  crater  depth  increases  initially  and  then  decreases. 
Eventually  no  crater  is  f orated.  The  data  free  Ft.  Folk  casts  agree  with 
the  sandstone  curve. 
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In  Figs.  8  and  9  we  present  data  on  the  cratering  ability  of  > 1000-lb 
and  larger  spherical  charges  in  Colorado  soils  and  different  percentages 
of  soil  saturation.  These  were  scaled-down  versions  of  DNA  tescs,  Middle 
Gusts,  and  Mixed  Company.  ORES  is  the  Defence  Research  Establishment, 

Canada.  The  effects  of  wet  and  dry  soil  are  evident.  Additional  specific 
data  are  given  in  Appendix  D. 

The  shape  factors  for  che  various  craters  are  shown  in  Fig.  10, 
where  the  crater  shape  is  normally  parabolic  or  ellipsoidal.  For  artillery 
shells,  the  elliptical  crater  has  been  suggested. 

Since  we  are  interested  in  cratering  efficiency,  we  show  this  parameter 
as  a  function  of  burst  height  (in  charge  radii)  in  Fig.  11.  Clearly,  for 
most  situations  where  artillery  or  mortars  are  used,  the  height  of  burst 
lies  between  0  and  1,  since  the  charge  radii  are  only  several  inches. 

Since  we  are  interested  in  scaling  charge  weights  4,100  lb  then  ve  should 
establish  the  range  of  validity  for  scaling.  In  Ref.  7,  experiments  were 
performed  with  0.08  lb  to  0.247  lb  of  HE  in  simulated  alluvial  soil.  For 
HOB  =■  0,  che  scaled  crater  radius  and  depth  are  0.69  and  0.62,  respectively. 
This  result  agrees  reasonably  well  with  the  data  obtained  with  fraction 
to  multi-ton  explosive  charges.  Similarly,  for  other  conditions  we  can 
assert  the  scaling  is  applicable. 

In  Table  13  we  have  estimated  che  parameters  needed  in  this  model 
for  several  soil  types  where  Item  1  is  the  values  used  for  the  Grafenwohr 
tests  in  November  1978.  The  values  for  VS MR  are  limited  since  the  moisture 
content  is  not  available  at  this  time.  We  caution  the  reader  these  are 
current  values  and  additional  data  can  modify  them  significantly. 


Seifert,  K.,  et  al..  Experimental  Report  of  Cratering  Displacement  and 
Ejecta  Processes,  Physics  International  Report  DNA  4482T,  August  1974. 
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TABLE  13 

ESTIMATED  PARAMETERS  FOR  SHELLS  AND  WARS 


Itaqt 

Soil  Type;  Condition 

c  ,..1/3 
u  r ,  ■  w 

cd/Wl/3 

-  f  r  c  3 /  1 

a 

hj.  jl 

l 

* 

Sod  it  soil;  vet 

.  7V  • 0  3 

.7W-04 

.55 

90% 

Sod  it  soil;  dr/ 

.7:w04 

i.ow"*03 

.33 

75% 

2 

Clay;  wee 

1. 44 

.72 

2.0 

90% 

Clay;  dry 

.  5 

.3 

.72 

752 

3 

Sandy  clay;  wee 

1.6 

.  56 

2.5 

90% 

Sandy  clay;  dry 

.55 

.32 

.75 

752 

4 

Silc  scone;  wee 

■  .04 

.43 

.55 

90% 

Silt  scone;  dry 

.9 

_  2 

.32 

752 

3 

WSMR  soil;  wee 

1.3 

.50 

1.67 

902 

W3MR  soil;  dr/ 

.96 

.45 

.37" 

752 

6 

Alluvial  soil;  moi3t 

.71 

.31 

.25 

50% 

* 

Cr3te? 

nay  be  larger  ae  the 

sure see.  we 

assume  here  an 

er  receive 

crater. 

135-bub 

live  shell. 

3.: 


thermal  effects 


3.2.1  General  Discussion 

the  detonation  of  a  high-explosive  charge  generates  a  high— temperature 
environment  including  shock-heated  air.  The  environment  is  a  high- 
temperature,  low-pressure  hutiole  (fireball)  which  expands  and  rises.  It 
is  responsible  for  the  later  stages  of  the  dust  clrud  behavior  and  the 
radiant  power  history.  Hera  we  will  discuss  the  •  'mieracuce  sources  and 
the  emitted  thermal  radiation  and  attempt  to  estimate  the  lace,  low 
temperature,  characteristics  of  the  cloud  chat  are  important  to  thermal 
electro-optical  systems. 

First  let  us  consider  the  shock  wave  effects  on  the  ambient  air  chat 
1/3 

occur  at  times  <10*^  msec.  Several  significant  variables  are  given  in 
Table  14. 


TABLE  H 

SPATIAL  DISTRIBUTION  OF  SHOCK-HEATED  AIR  rQft 
ATMOSPHERIC  DETONATION 


Seduced 

Radial 

Distance 

< 

?aak  shock 
Overpressure 
fMt) 

Peak  Shock 
Temperature 

CO 

Residual 

Temperature 

CC) 

:,93 

250 

911 

245 

4.32 

iOO 

399 

101 

5.52 

so 

225 

53 

6.3e 

20 

153 

37 

7.59 

20 

114 

30 

10.0 

10 

7  4 

2e 

The  shock-heated  air  will,  be  engulfed  by  the  dust  cloud  at  later  times. 

1/3  1/3 

As  will  be  shown  later,  the  cloud  radius  is  — 5.5W  at  t  =  Q.1W  sec 

1/3  1/3 

and  >  10W  (ft)  at  t  — 0.35W“  sec.  Since  we  are  concerned  about  longer 
1/3 

time  periods  (— 4W  sec)  ,  the  rireball  etfect  or  heated  cloud  is  of 

greater  significance. 

The  luminous  cloud  has  a  definitive  thermal  pulse — a  first  peak, 
second  peak,  and  a  later  third  peak  resulting  from  the  rising  TNT  debris 
afterburning.  We  can  write  for  the  peaks  the  following  scaling  laws: 

1/3 

t  =  I2WX/J  (psec) 
t2  =  150W1/3  (usee) 
t^  =  2.5W3//3  (msec) 

The  radiated  power  (for  an  HOB  =  1)  normalized  to  t2  is  shown  in 
Fig.  12,  and  the  integrated  value  in  Fig.  13.  Although  the  latter  implies 
that  the  normal  radiation  is  —  10%  of  the  blast  yield,  the  probable  value 
is  closer  to  30%  since  the  sensors  measure  primarily  visible  and  near- 

g 

infrared  radiation. 

Let  us  now  consider  surface  heating.  For  a  1000-Ib  charge  the  fire¬ 
ball  radius,  at  1U  msec,  is  —3.6  m.  For  a  500-ton  charge  we  observe  a 
value  of  —39  ro  at  100  msec.  Thus  the  scaling  agrees  with  the  observed 

results.  These  times  correspond  approximately  to  times  of  their  respective 
second  peak.  Hence,  if  scaling  holds,  we  have  a  1-lb  charge,  a  sphere 
(fireball)  whose  radius  is  0.36  m.  At  this  time  scale  —10%  of  the  total 
energy  is  radiated. 

For  a  155-mrn  shell,  W  =  7  kg,  the  fireball  radius  is  0.88  m  at  2.26 
msec.  Hence,  the  energy  radiated  is 
* 

It  should  be  noted  that  it  takes  about  6  usee  for  the  detonation  wave  to 
reach  the  surface  of  a  1-lb  spherical  charge. 

g 

See  discussion  in  F,  B.  Porzel,  Introduction  to  a  Unified  Theory  of 
Explosions,  NOL  TR-72-209,  September  1972. 
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(18) 


Q  *  (0.1) (0.1)W  x  106  cal/kg 
=  7  x  104  cal 

where  the  thermal  power  is  10%  of  the  explosive  yield.  The  ground  ir- 
radiance  is  therefore 

H  3  — cal/cm^  «  —  *■— r  (19) 

Trr  4tt(88)Z 

2 

Consider  this  input  to  1  cm  of  surface  of  thickness,  t;  then  the  temperature 
rise  is 


2 

0.72  cal/cm  =  oc  AT 
P 


(20) 


where  p  =  soil  density 
0^  *  specific  heat 

or,  for  t  =  1  cm 

AT  =  0. 72/(2. 6) (0.22)  =  1.26°C  (21) 

Since  the  irradiation  occurs  at  shorter  times  and  radii,  the  increase  may 
be  <,3°C  which  under  most  conditions  is  unimportant. 

Thus  we  suggest  that  the  surface  heating  for  explosive  events  of 
interest  to  this  study  (artillery/mortar  munitions)  is  unimportant  in 
calculating  the  overall  cloud  thermal  effects.  There  is  semi-quantitative 
evidence  that  the  crater  is  warm.  Long-term  observations  (42  min)  after  a 
multi-ton  event  showed  that  in  the  near  field,  i.e.,  2*  crater  radius,  the 
surface  temperature  was  characteristic  of  the  subsurface  ambient  temperature 
and  then  dropped  to  the  normal  ambient  surface  temperatures  at  5*  crater 
rad ius .  This  implies  that  background  temperatures  can  be  disturbed  by  the 
explosion  for  long  periods,  and  furthermore,  may  depend  upon  the  season 
of  the  year. 
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The  work  done  by  che  explosive  on  che  soil  and  the  hoc  cloud  fallout 
into  che  crater  can  contribute  to  che  warm  crater  temperature. 

3.2.2  Cloud  Temperature 

Above,  we  have  shown  that  the  ground  heacing  and  shocked  heated  air 
is  not  significant  since  che  fireball  heacing  of  the  dusc  cloud,  which 
includes  the  latter,  is  che  dominant  mechanism.  The  temperature  history  is 
important  for  two  reasons:  (1)  residual  heat  to  generate  a  buoyant  cloud; 
and  (2)  target  obscuration. 

Quantitative  measurements  of  thermal  radiation  produced  by  nigh- 
explosive  detonations  are  limited  and  these,  until  recently,  to  relatively 
large  events.  In  Fig.  14  we  present  the  power-time  results  from  a  500-ton 
event  in  che  spectral  bands  a  4  2.5  um  and  \  <_  4.5  urn.  The  effect  at  —4 

aSr 

sec  is  due  to  either  afterburning  or  uncovering  of  low  temperature  debris. 

Correspondingly,  the  mean  temperature  of  cloud  vs.  time  is  shown  in  Fig.  15. 

The  extrapolation  to  lower  temperatures  is  a  complex  problem;  however, 

this  has  been  done  analytically  (see  Fig.  16)  and  can  be  used  for  other 

1/3 

charge  weights  by  using  the  scaling  factor  (y/15.6)A 

The  cooling  of  che  hot  -*  warm  dusc  cloud  is  given  by  a  relationship 
of  the  form 


2.T  *  2T0  exp  [-(3h/sc  r)c]  (22) 

where  h  «  thermal  transfer  coefficient  (cal/cm"  see  4C) 

o  »  particle  density,  g/cc 

c  -  particle  specif ie  heat,  cal/g  ’C 
P 

r  «  particle  radius,  cm 


We  note  here  thn  aven  at  t  »  30  sec,  th-are  is  an  amount  of  thermal 
radiation  emittea  at  wavelengths  of  interest  to  some  thermal  systems. 
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The  exponent  i 3  -(3  *  1.356  x  10^/(2.6)r  =  6.51/r,  where  r  is  in 
micrometers,  and  therefore,  the  (AT^/AT^)  for  r  =  10  um  is  —0.521/1.49  x 
10  ^  —  1/350  and  1/671  for  AT  ^  ^  /  AT  ^ .  These  values  do  not  disagree  with 
observed  values  for  the  155-mm  event,  and  the  curve  scaled  from  large  TNT 
charges;  see  Fig.  16.  For  r  different  than  10  um,  the  times  may  be  faster 
or  slower  if  les3  or  greater  than  10  pm,  respectively. 

For  buoyant  plume  calculations  we  assume  that  300  cal/g  of  TNT  is 
the  appropriate  quantity  and  the  duration  is  t  (sec)  4.1W  '  J,  This 
quantity  does  not  scale  precisely  since  for  larger  events  the  cloud  pene¬ 
trates  to  greater  altitudes,  i.e.,  to  lower  temperatures,  giving  greater 
buoyancy.  For  example,  a  100-ton  event  was  reported  to  stabilize  at  an 
altitude  of  —3000  ft  after  300-480  sec,  whereas  the  above  expression 
would  suggest  —240  sec. 

Alternatively  we  can  use  Fig.  16  in  the  following  fashion:  if  the  mean 
cloud  temperature  >  AT  of  the  ambient  air  temperature  corresponding  to 
cloud  Az  (temperature  at  cloud  top  minus  temperature  at  cloud  base) ,  it  is 
buoyant. 

Although  the  above  is  directly  applicable  to  an  explosion  on  the 
surface,  extension  to  near-surface  bursts  will  introduce  smaller  errors 
than  other  uncertainties.  For  example  in  the  HOB  =  -1  case,  the  high- 
temperature  bubble  will  be  non-spherical  due  to  the  initial  venting  process. 
Also  the  crater  will  be  hotter  chan  for  HOB  «  +1. 

For  events  like  Fort  Sill,  i.e.,  small  particles,  the  cooling  should 
be  faster. 

3.3  CLOUT)  DISPERSAL 

In  Ret.  6  the  generation  of  the  dust  cloud  and  dispersion  of  this 
cloud  are  described  in  detail  for  the  case  of  HOB  ®  1.  The  results  are 
given  in  Sec.  3.4;  that  is,  for  the  lofted  crater  mass. 
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3.4  INITIAL  CLOUD  VOLUME 


The  sec.cnd  stage  is  defined  when  an  initial  volume  is  lofted  into 
che  atmosphere,  forming  a  parabolic  dust  cloud.  This  stage  occurs  within 
a  scaled  time,  c  £  0.  1W1'3.  Hence,  the  volume  ia 


V  . 

C  4 


(C) 


a  Dct 


K. . (c) 


(21) 


where  D  H  ,  are  specified  in  Fig.  17. 

cl  cl 

For  example,  at  t  -  0.1W1/3,  D/Wi/3  •  10.3  f c/ (lb) L/ 3  and  H/W'l/3  - 
6.6  ft/{lb)3//3j  therefore 

V  (0.1W1'*)  -  32GW  (24) 


Hence,  for  an  explosive  charge  of  weight  W  »  15  lb  for  wet  terrain,  the 

3  3  * 

concentration  C  is  ~-3.5  *  10  g/m  . 

3.5  GROWTH  OF  MAIN  CLOUT)  AND  STEM 

For  the  time  period  0.1W^3<  c  -  Q.3W^3,  we  treat  che  growth  of  che 
stem  and  the  main  cloud  separately. 

The  stem  diameter  is  considered  to  oe  cs  2  crater  diameters.  We  con¬ 
sider,  for  computational  reasons,  that  the  stem  is  a  cylinder. 

The  main  cloud  is  also  defined  as  a  cylinder  whose  volume  is  given  by 
che  expression 


T  3,  .'CO 

4  ;  t 


1(0 


(25) 


where  T(c)  is  che  cloud  thickness;  see  Fig.  17.  The  stem  value  is  defined  as 


* 

Sate 


fcAib) 


3 


»  1.6  -  10* 


lid 


(26) 


V  (c)  -  4t  C2  •  H  (c) 

3  r  3 

where  H  (t)  is  the  seem  height  and  given  in  Fig.  13.  The  overall  cloud 
height  is  H  (t)  +  T(t). 

At  this  time  we  assume  that  the  cloud  has  stabilized  at  a  height 
defined  by  W  only;  chat  is,  we  include  the  effecc  of  lapse  rate  as  part  of 
the  plume. 

The  buoyant  plume  in  still  air  is  a  special  case,  see  discussion  in 
Part  1.  On  che  basis  of  limited  data  and  analysis  one  may  estimate  that 
the  height  of  che  cloud  is  given  by  che  value 

H  -  3.7FL/4  (nO  (27) 

-5  3 

where  F  -  3.6  *10  Q  (cal/sec)  and  Q  •  (10  cal/gm  TNT)  •  454W  f 
0.3U2,/2  or  Q  ■  (1.5  *  10^  cal)W“^2,  where  W  is  in  pounds. 

3.6  (.LOUD  DISPERSION 

For  times  greacet  than  c  a  Q.3W2^2  we  assume  che  dust  cloud  is 
dispersed  like  a  low-order  thermal,  stationary  source. 

For  the  present  we  limit  the  discussion  for  dust  cloud  dimensions 
at  or  near  stabilization  when  HOB  *  0  and  -l.  These  are  listed  iu  Table  15. 
Like  the  dispersion  of  the  HQ3  a  i  case,  the  dust  clouds  for  HOB  »  0  and 
-l  are  dispersed  by  the  methodology  described  in  Part  1.  This  includes 
the  thermal  effect  described  in  See.  3.2,  namely  300  cal/g  of  TNT  is 
available  for  a  time  t  *  4.1W f  sec. 

Again  we  assume  to  a  first  approximation  the  mass  concentration  is 
uniform.  In  the  next  subsection  we  consider  qualitatively  the  subject  of 
non-uniform  distribution. 

: 

Sirkind.  R..  An  Obscuring  Aerosol  Dispersion  Model.  Cenerai  Research 
Corporation  Report  CR-231,  December  1^78.  ” 
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TA3U  >5 

CLOUD  DIMENSIONS  FOR  HOC  =  0,  -1 


HOB  --  0 

HOB  -  • 

Cloud  Heighv  (fc/Wi/3) 

7 

6.2 

Cloud  DUmecer  (tc/W^'^) 

6.84 

5.3 

item  Diameter  (fc/W^^) 

1.14 

1.0 

’  /"» 

Stem  Height  ffc/W*  ) 

2.28 

2.0 

3.7  MASS  DISTRIBUTION 

One  of  che  uncertainties  in  explosive- related  phenomena  is  the 
spatial-temporal  history  of  che  lofced  dusc.  Although  che  earlv  history 
is  responsible  for  che  flov  field  internal  to  che  dusc  cloud  and  che  dusc 
distribution,  our  inceresc  is  in  che  later  stages.  Experimental  data  are 
limited  particularly  for  artillery /rnorcar  munitions;  however,  che  general 
behavior  is  reasonably  understood.  It  should  be  noted  chat  che  data  from 
Urge  KE  detonations  are  not  directly  applicable  in  view  of  che  heights 
reached,  thousands  of  feet  above  ground  level,  to  which  che  dusc  cloud  rises 
the  vertical  atmospheric  wind  3nd  temperature  structure  are  different. 
Nevertheless,  che  creation  of  a  cloud  stem  and  a  rising  vortex  in  the  main 
cloud  are  also  characteristic  of  near-surface  explosions. 

Observations  for  Urge  detonations  after  the  initial  rise  indicate 
that  the  horisontai  distribution  in  the  cloud  Has  a  fora  shown  below;  that 


I 


is,  the  primary  mass  is  located  in  the  vicinity  of  the  vortex  center,  and 
a  relatively  small  amount  of  dust  is  concentrated  around  the  axis.  The 
stem  distribution  will  have  a  form  similar  to  the  one  shown  below:  the 
stem  extremities  (x)  are  increased 


due  to  limited  quantities  of  upswept  ground  dust.  The  vertical  distribution 
is  a  result  of  the  vertical  motion  and  the  sedimentation. 

3.8  MASS  CONCENTRATION 

The  swept-up  dust  is  part  of  the  stem  and  increases  the  outer  radius 
somewhat . 

From  the  above  sampling  results  the  vortex  center  appears  to  exist 
equidistant  from  the  axis,  and  we  may  assume  the  vertical  location  is  the 
center  of  the  cloud.  In  time,  turbulence  will  vertically  mix  the  particulates 
in  the  cloud.  On  the  other  hand,  sedimentation  (fallout)  plays  an  important 
role  in  altering  the  mass  ana  size  distribution  in  the  cloud  and  stem;  that 
is,  the  larger  particles  settle  out  by  inertial  forces.  This  is  illustrated 
in  Fig.  19  for  a  100-ton  event,  HOB  =  1.  The  mass  values  given  are  for 
a  horizontal  column  of  unit  area  at  that  altitude  and  time. 

On  the  basis  of  the  above,  for  charge  weights  of  battlefield  munitions, 
we  consider  chat  for  normal  wet  soil  about  90%  of  the  soil  ejected  will 
settle  early  and  therefore,  10%  will  constitute  the  dust  cloud  and  stem. 

In  fact  most  particles  in  excess  of  100  ym  will  fall  out  in  the  first  10 
sec.  On  the  other  hand,  if  the  soil  in  the  dust  cloud  is  made  of  small 
particles,  as  exemplified  by  Fort  Sill,  the  large  fraction  of  the  ejected 
mass  may  remain  aloft. 

If  there  is  a  need  to  establish  a  non-uniform  concentration  then  for 
the  main  cloud,  above  the  stem  height  plui  3  m,  we  suggest  the  vertical 
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ALTITUDE  (:<10  ft) 


9 


•  -  MASS  (g  •  10'H) 


o 

31 


40 


Figure  19.  Cloud  History  for  Middle  Gust  III;  (100  T) ,  HOB  =  1 
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mass  distribution  is  of  Gaussian  form  with  the  halfwidth  about  the  cloud 
toroid  center  equal  to  toroid  thickness/2.  Similarly,  for  the  horizontal 
distribution,  we  may  assume  the  shape  as  shown  above  with  the  same  half¬ 
width  as  the  vertical  distribution.  Since  turbulent  vertical  and  horizontal 
mixing  will  exist,  uniform  concentration  occurs  after  a  time  —20-30  sec. 

For  the  model  presented  here  we  assume  uniform  concentration  in  the 
cloud  and  stem;  however,  the  mass  differences  will  yield  different  concentra¬ 
tions  for  both. 

Since  the  height  of  the  clouds  attain  values  of  30-40  m,  we  disperse 
the  cloud  in  the  following  manner:  (1)  stem,  (2)  top  layer,  and  (3)  cloud 
center.  Obviously  the  cloud  height  is  established  from  the  top-layer 
behavior,  the  overall  width  from  the  stem  and  cloud  widths.  The  dispersion 
algorithms  require  the  mean  wind  speeds  at  the  respective  heights.  Thus 
we  use  an  average  u  for  the  stem  at  —2m.  To  find  the  values  at  the  dif¬ 
ferent  altitudes  we  can  suggest  the  following: 

V 

Stem:  u(z)  *  In  (z/z^)  (28) 

-2  -2 

where  V  =  frictional  velocity  =  (1/2)C  ,u.  (C  ,  —  4  x  10  ) 

x  d  i.  u 

Zq  =  surface  roughness;  see  Table  8,  Fart  1 
u  *  mean  speed  at  0.5  m 

For  heights  >4  m,  employ  the  power  law: 

u  *  u1(z/z1)P  (29) 

where  z^  is  4  m  and  is  the  mean  speed  at  4  m.  The  values  for  p  are  the 
following:  neutral  -  0.142;  unstable  -  0.1;  stability  G  -  0.8.  Values  of 
0.3  and  0.5  are  suggested  for  stability  classes  E  and  F. 

Although  a  surface  surge  is  associated  with  the  shock  wave,  we  will 
neglect  this  effect  for  several  reasons: 


1.  The  height  of  the  surge  cloud  can  be  expressed  as 

H  (ft)  =  2*AP  •  t  (sec) 

and,  since  we  are  interested  in  APs  ~  0.1  psi  and 
t  —  .003  sec,  H  (ft)  =  2.  (0.1) (0.003)  =:  0.  Further, 
the  shock  strength  would  be  too  weak  to  raise  signi¬ 
ficant  dust  from  the  compacted  soil. 

2.  The  material  velocity  of  the  air  behind  the  shock  wave 
is  only  —  1.7  m/sec. 

3.  The  effect  of  the  asymmetric  jets  should  be  more 
important  on  the  near-surface  obscuration. 

Hence  any  dust  raised  by  the  surface  surge  will  for  modeling  purposes 
be  considered  as  an  integral  part  of  the  initial  cloud. 


3.8.1  Settling  Speed 

The  particulate  material  within  the  cloud  will  settle  out  provided 
the  vertical  downward  velocity  exceeds  the  upward  wind  speed  where  the 
latter  is  approximately  one-fifth  the  mean  horizontal  wind  speed. 


Consider  a  spherical  particle  to  be  acted  on  by  gravity.  The  force 
is  then  defined  by 


where  r  =  particle  radius,  cm 

2 

g  =  981  cm/sec 
pp  =  particle  density,  g/cc 
p  =  air  density,  g/cc 

cL 

For  particles  at  constant  velocity  we  can  apply  Stokes  Law  to  define  the 
resistance  force;  i.e., 


R  =  67m  •  r  *  v 


(31) 


where  r,  =  coefficient  of  viscosity  (1.8  x  10  g/cm  sec) 
v  =  particle  velocity 

At  this  point  we  need  to  examine  the  character  of  Stokes  Law.  The 
drag  can  be  defined  by  the  expression 

R  -  ( 1/ 2 ) p  v2C  A  (32) 

a  D 

where  CQ  3  drag  coefficient 
A  =  particle  area 


or,  for  the  Stokes  case, 


CD  3  24/Re  (33) 

where  Re  is  Reynolds  number  3  pvD/r,.  One  may  show  when  Re  >  1  that  the 
Stokes  value  for  CQ  underestimates  the  drag,  and  a  more  appropriate  expres¬ 
sion  is 


0.06  + 


36 

Re 


(34) 


To  find  the  settling  speed  for  Stokes  particles,  neglecting  0  ,  we 
have  F  -  R  ■  0  or, 


2  a  2 
7  4  0  r 


V  ■  rr  ■**  q 

9  n  P 


(35) 


Substituting  Che  values  for  g,  r,  and  0pr*’»  we  have  a  2  g/cc  particle, 


6  7 

v  •  2.5  x  10  r“ 


(36) 


Hence,  for  a  10- urn  particle,  v  ■  2.5  cm/sec.  For  a  100-um  particle,  v  • 
2.4  m/sec.  Experiments  have  shown  that  the  terminal  velocities  for  100-um. 
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1000-uo,  and  1-cm  particles  are  0.5  m/sec,  6  m/sec,  and  22  m/sec,  respee- 

if 

lively.  For  our  use  here  we  assume  Stokes  Law  for  r  <  50  um  (v  <.  0.5  m/sec 
and  che  experimental  values  for  r  »  50  urn.  It  this  is  the  situation,  then 
in  10  sec  all  particles  with  r  z  50  urn  at  heights  <  5  a  will  settle  out. 

Alvarez^  examined  che  effect  of  relative  humidity  on  precipitation 

«4* 

of  ambient  hygroscopic  particles.  His  data  indicate  chat,  if  che  relative 
humidicy  >  74*,  the  suspended  particles  approach  precipitation  weight  and 
by  91 %  che  particles  are  heavy  enough  to  precipitate.  This  situation  may 
be  exascerbated  in  a  dust  cloud.  Again,  there  is  little  data  on  che  relaciv 
humidity  within  che  dust  cloud;  however,  the  local  value  must  be  elevated 
above  ambient  conditions. 


Hi 

The  general  equation  of  motion  is  defined  as 


y  -  -g  -  (18u/cD')y 
x  »  -(18u/dD4)x 


from  which  we  obtain 

^  u  ,,  -at. 
*  *  *0  *1  (l  "  e  } 


y  *  yQ  ~  +  <Va  +  3/Od 


-at. 
e  ) 


where  u.»  mean  horizontal  speed,  VQ  »  initial  vertical  velocity,  and  a  • 

18  \xlOu  -  18  1/D2. 

^Alvarez,  R. ,  "A  Study  in  the  Changes  on  the  Quantity  of  Aerosol  Particles 
in  the  Ambient  Medium  as  a  Result  of  Fluctuations  in  che  Relative  Humidity 
of  che  Air,"  lav,  Atmospheric  $  Oceanic  Physics,  Vo  1 .  13.  So.  12.  1377. 
p.  911. 

It  should  be  noted  chat  this  effect  on  hygroscopic  smokes  has  not  been 
examined  and  say  affect  the  yield  of  phosphorus-derived  smoke  clouds. 

Also  ambient  partieles  are  <,  10  urn.  whereas  dust  parcieles  are  primarily 
>  10  LB. 
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VALIDATION  TESTS 


As  in  Part  1,  we  present  here  several  tests  or'  the  dust  model  dis¬ 
cussed  earlier.  Unfortunately,  the  many  tests  performed  to  data  have 
limitations;  that  is,  the  data  are  unavailable  (Grafenvohr  1979,  DIRT-1); 
the  data  are  not  complete  (Grafenvohr  1973  ar.d  others);  for  live  firings 
the  exact  detonation  location  relative  to  the  samplers  are  unknown.  Further, 
few  tescs  of  single  or  multi-individual  rounds  are  available  for  analysis. 

If  the  test  consists  of  several  adjacent  rounds,  <.  25  m*,  then  interactive 
affects  occur  that  are  difficult  to  model  accurately.  Notwithstanding 
the  above,  we  will  present  several  cases. 

Test  1  -  Grafenvohr  -  November  1978;  NV  £  EO  LAB 

Input  Conditions: 

u  -  3.35  m/sec 

Pasquill  Category  -  C 

Ground  level  -  960  mb 

Munition  -  One  155-nxm  shell  artillery  delivered 

The  data  provided  were  limited  to  semi-quantitative  cloud  geometry 
and  the  transmission  tv  the  8-12  urn  spectral  region. 

In  Fig.  20  we  show  that  the  cloud  height  is  accurately  predicted 
when  the  buoyancy  (heat  added)  is  taken  into  account.  In  Fig.  21.  we  present 
the  calculated  concentration  and  transmission  and  the  measured  transmission. 
Although  the  practical  implication  between  the  two  transmissions  are  small, 
the  measured  value  is  2*  the  calculated.  There  are  several  possible  ex¬ 
planations:  (l)  the  extinction  coefficient  is  too  large;  (2)  the  con¬ 
centration  is  too  high;  or  (3)  both.  The  underlying  reason  is  that  the 
transais seme ter  is  naar  the  cloud  base  and  consequently,  the  particle  sices 
may  be  different  (larger)  chan  in  the  other  parts  of  the  cloud,  and  the  parts 
of  the  line  of  sight  are  not  filled  with  dust. 

ft 

Half  distance. 


Figure  20.  Cloud  Height  as  a  Function  of  Time 


I 
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T e sc  2  -  Grafanwohr  -  N ova mb e r  197 3 


Input  Data: 

u  -  2.5  m/sec 

Pasquill  Category  -  C 

Visibility  -  2  km  (high  humidity) 

Ground  temperature  -  l’C 

Munitions  -  Two  155-mm  artillery 

shells  delivered;  detonations  25  m  apart 


Results 

Apparently  the  two  events  occurred  simultaneously,  and  therefore,  the 
shock  interacts  at  '-12.5  m  at  97  msec.  A  Mach  shock  will  form  at  an 
altitude  of  —5  m  and  then  propagate  in  both  directions.  The  wave  will 
reach  the  ejecting  mass  at  125  msec  at  a  velocity  of  -40  m/sec.  The  wave 
will  interact  with  the  dust  cloud  ..‘or  0.02  sec  in  the  direction  of  air  flow. 
Vectorial  considerations  yield  an  average  displacement  of  1.7  a  due  to  the 
wind;  i.e.,  the  average  ejecta  velocity  is  -40  m/sec  at  45*  to  the  vertical. 
The  horizontal  displacement  a  and  43  m,  and  the  cloud  diameter  is 
increased  -3.4  a.  The  effect  on  cloud  height  is  small,  4 1  m. 

The  calculated  initial  cloud  dimensions  are  as  follows: 

t  (sec)  S  (m)  H  (a) 

0.031  3.72  1.49 

0.122  5.53  2.79 

0.742  13.2  10. 0 

Subsequent  dimensions  are  determined  by  conventional  expressions  tor 
Pasquill  Category  C. 

In  Fig.  22  the  claud  height,  as  measured  by  a  silicon  TV'  and  3-12  ua 
Ft ISt »  is  presented  as  a  function  of  time.  The  calculated  values  are  pre¬ 
sented  as  well  as  the  settling  location  for  particles  4IOQ  „s.  The 
decrease  of  the  height  at  infrared  wavelengths  is  due  to  the  cooling  of  the 
dust  cloud. 
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Figure*  22.  Cloud  Height  History  in  iiilicon  and  8-12  >iui  Region 


The  cloud  width  comparative  results  are  shown  in  Fig.  23.  These  and 
the  above  show  good  agreement.  One  must  recognize  chat  the  measured  results 
are  subject  to  some  error. 

Finally,  the  calculated  and  measured  transmissions  through  the  cloud 
are  presented  in  Fig.  24 .  The  thermal  transmission  agrees  reasonably  welt; 
however,  the  silicon  region  appears  to  have  a  lower  extinction  coefficient 
than  theoretically  estimated. 

Test  3  -  SWII  -  Trial  29 

The  input  data  are  given  in  Fig.  25,  where  the  second  row  was  detonated 
5  sec  after  row  A.  The  distance  of  S,  C,  and  M  transmissometers  were  40, 

90,  140  m  from  row  A,  respectively. 

The  calculated  and  measured  C  •  L  for  the  three  cransmissoaracers  are 
shown  In  Fig.  26.  It  should  be  noted  that  the  shock  interaction  has  beer; 
taken  into  account;  i.e. ,  intra  and  inter  row. 

The  agreement  between  the  two  C  •  t  curves  is  excellent. 

Test  4  -  Ft.  3111,  DPt-005-n 

In  Fig.  27,  we  present  the  initial  conditions,  the  measured  C  •  I. 
(solid)  and  the  calculated  results.  If  the  arrival  time  is  displaced 
several  seconds  then  the  agreement  is  improved.  The  methodology  used  here 
is  as  described  in  the  text. 


Test  3  -  Ft.  Sill,  DPI-QQ5-T19 

The  available  input  conditions  are  given  below: 


u 

Relative  Suaidlty 

Time 

Sky 

Pasduill  Category 
^initio* 


6  a/ sec 

1141 

Overcast 
iftstisuated)  0 

1.  K-l.  IQS-bs  artillery  shell 
delivered 


15- 


■«**-*: 


r %?**rg*K**&* 


Fifewrc  2'i.  Cloud  Width  lUalury  in  Silicon  and  Thermal  Region  Measurements 


In  Fig.  28  we  present  the  calculated  and  observed  cloud  height  and 
width.  The  results  are  in  fair  agreement  in  spite  of  the  limited  input 
data  to  establish  the  prevailing  meteorological  conditions.  In  Fig.  29 
we  present  the  comparative  results  for  C  •  L,  and  the  results  are  in 
good  agreement. 

On  the  basis  of  the  above  we  conclude  that  the  proposed  model  can 
provide  a  good  estimate  of  the  dust  obscuration  generated  by  explosive 
munitions. 

With  respect  to  the  detailed  explosive  dust  model  presented  here, 
a  form  suitable  for  operational  use  is  described.  For  specific  calculations 
where  the  actual  burst  height,  soil  type,  etc.,  are  known,  the  reader  can 
obtain  the  necessary  cratering  expressions  from  the  basic  text. 


160 


Figure  28.  Comparison  of  Calculated  and  Measured  Cloud  Height  and  Width 


Figure  29.  Comparison  of  Calculated  and  Measured  C 
History  (•  -  model  predictions) 
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MODEL  STRUCTURE  FOR  EXPLOSIVE  DUST 


I 

f 


i 


5.1  GROUND  INPUT  DATA 
A.  Scenario 

Line  of  sight 
Munition  deployment  plan 
Munition  rate 

Sensor  type  and  operational  wavelength 

3 .  Munition  Characteristics 
Type  (caliber,  Mk  No.) 

HE  material  and  weight 
Fuze  cype  and  setting 

C .  Meteorological  Conditions 
Time  of  day 
Cloud  cover 

Mean  wind  speed  (m/sec)  -  surface  to  10  m 
Temperature  T(°C)  -  0.5  a  and  10  m 
Elevated  inversion  layer 
Wind  direction 


0*  Terrain  Characteristics 
Soil  moisture  conditions 
Soil  type  and  gradation 
Vegetation 

E.  Optical  Properties 

(1)  Effective  Charge  Weight  Calculation 


■f  *  ^  (pressure  at  sea  level/pressure  at  detonation  alt) 
w  »  a  (explosive  weight) 
a  »  1  for  TNI 

•  1.19  for  Coop  B 
*1.1  for  C-4 

•  0.86  for  Explosive  D 
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(2)  For  ambient  conditions  and  explosives  used,  compute  V 


and 


1/3 


err 


(3)  Establish  soil  type  and  condition,  enter  Table  16,  and 

find  scaled  parameters,  "  and  m^.  From  (2)  obtain  crater 
radius  and  depth.  The  latter  is  only  for  general  Interest. 


TABLE  16 

ESTIMATED  PARAMETERS  FOR  SHELLS  AND  MORTARS 


I  cam 

Soil  Type;  Condition 

cr/v1'3 

crt/V3 

"  (ft3/'/) 

* 

907. 

I 

Sod  and  soil;  v* t 

0.7V03 

0.71W"'04 

0.55 

Sod  and  soil;  dry 

Q.72W.04 

l.ou’*03 

0.33 

757 

2 

Clay;  vet 

1.44 

0.72 

2.0 

907 

Clay;  dry 

0.50 

0.30 

0.72 

757 

3 

Sandy  clay;  vet 

i.60 

0.56 

2.50 

907 

Sandy  clay;  dry 

0.55 

0.32 

0.75 

757 

4 

Silt  stone;  vet 

1.04 

0.43 

0.55 

907 

Silt  stone;  dry 

0.90 

0.20 

0.  32 

757 

3 

USMR  soil;  vet 

1.30 

0.50 

1.67 

907 

WSJ®  soil;  dry 

0.96 

0.43 

0.37* 

757 

6 

Alluvial  soil;  moist 

0.71 

0.31 

0.25 

507 

Fallout  fraction  (7). 

Crater  nay  be  larger  at  the 

surface,  ve 

assume  here 

an  efi'eet 

ive  crat 

J155-esi  live  shell. 
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(a)  Determine  ejected  weight 
"  -  See  Table  16 


charge  weight  (lb) 


'Lb/fc^)  -  soil  type 


100 

Alluvium  moist 

35 

Sod,  dry 

70 

Sod,  wet 

125 

Clay,  wet 

100 

Clay,  dry 

100 

Sandy  clay,  dry 

125 

Sandy  clay,  wet 

140 

Sandy  clay  and  gravel 

125 

Siltstone,  wet 

90 

Siltstone,  dry 

(b)  Lotted  weight  *  (1  -  m^)  •  ejected  weight 


(t)  Main  cloud  dust  weight  -  (0.9)  •  lofted  weight 
Stem  dust  weight  =  (0.1)  •  lofted  weight 

P-etain  above  for  use  to  find  concentration  and  part  icle 
distribution,  if  desired. 


(*)  Cloud  History  -  initial 


(a)  Eater  Fig.  30.  find  D^u)  and  H  ,  =  ,  for  c  <  0.1  wi/3 


lb)  Eater  Fig.  30.  for  Dc.(t).  T(t)  when  t  >  0.1  II 

and  <  0.3  Vl/i 


i/3 


(e)  Enter  Fig.  31,  for  stem  heigh,  and  ac^aa  diameter,  i.e., 


0g  •  4C 


(i)  ?roa  ««  get  Vo l  »  •  H  (t> 

(e)  ?rom  (b)  and  (t)  we  have  Voi  ^  (^/r)S’(c)  *  T(t) 


V  •  4if  C'  *  H  (t) 

^  r  d 


Jf2.ce:  2  »  *  ♦  a  st  t  •  0.3  W'^3  (sec) 
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T&C." 


Hal 


(5)  Concentration  -  Initial 


„  0.9  Loire  Weight  (?)  ,  3 

C  - - - - 3 - g/m 


t/4  •  D“.  •  I 

Ci. 


0. 1  Loft  Weight  (?)  3 

■  '  a  1  .a  i  rj  ,  m 


g/O 


*tC‘  *  H 
r  s 


We  next  consider  the  dispersion  of  the  main  cloud  and  stem,  realizing  chac 
che  cloud  is  scill  buoyant. 

(6)  Cloud  History  -  Buoyant  Phase 

(a)  Calculate  c.  *  resultant  is  time  for  buovanc 

0  " 

phase  of  rnajor  cloud  only. 

(b)  Calculate  F  -  3.6  -  10"0  (300  •  g  HE  r  c.  ) 

J 

(c)  Pasquill  Category  A.  5,  and  C 


.  ,  r 1/2  2/3 
1.65  x 


u  3t  •*  n 

(d)  Pasquill  Category  0 


(t)’l 


i  *  o.aii 


sa  -  14F 


i/6 


?  <  55  a4/ see 3 


*  3-? 


0 .( 


?  j.  55  s’*/ see 


3 


i  /  ■»  *  n 

(e)  (lh)  *  l.e,»  “s”'  j  u 

sax 

where  x  •  J.Sx*  for  Pasquill  Category  A.  B.  C.  aod  D 


in 


00 


(£)  Cat 


c 


F,  etc. 


i  'T 

Compute  3  =  -J-  7—  +  0.098  ;  X  is  actual  temperature 

T0  “Z 

_  ,.-,0.5 

Compute  x,  ■  1.-*  u/3 

Ah  ■  1.6  F^^x'  ^  r  u  tor  x  <  x^  u  s.  ?  4m 
_  1/2 

Ah  -  :.9  (F/us)  ' 
max 

(g)  u  -  0 


1. 37Q° ‘ 


(o) 

W  (kgj  ;  W 


TNT  equivalent  weight 


Noce:  The  dimensions  here  are  given  by  a  righc  circular 

cone  where  cop  radius  =  0.5  *  H  and  the  cone  hair-angle 

max 

is  16*. 

(7)  Stem:  Momentum  Dispersion  Phase  (u  >  0) 

(a)  From  input  data  and  the  t3ble  determine  Pasquili  Category. 


3;sq'j:li  s'a3:l:ty  CATiX^:is 


Suriaee  Kind 
Speed  at  10  s 
(»  aee  ) 


Insolation 

Strong  Moderate  Slight 


Might 

Thinly  Overcast 
or 

low  Claud  1  3/ 8  Claud 


1 

A 

A-S 

3 

— 

* 

2-1 

A- 8 

8 

c 

*’ 

3-S 

5 

3-C 

c 

3 

5> 

T 

5-a 

C 

e-o 

3 

li 

9 

>6 

3 

3 

3 

9 
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(b)  it  ■  uc,  u  n  a 

(c)  Cloud  Width,  y 

From  (a)  enter  cable  and  corapuci  v. 


Pasquill  Category  -  A 
Pasquill  Category  -  il 
Pisquiil  Category  -  C 
Pasquill  Category  -  D 
Pasquill  Cacegory  -  £ 
Pasquill  Category  -  F 


y  *  9.1  +  0  .  1 9 :< 

•  9.1  +  0.328* 

•  9.1  +  C.23P-X 
-  9.1  +  G.2Cx 
■  9.1  +  0 . ' 3  x 

•  9.1  +  0.  U6x 


(d)  Cloud  Rise,  s 

Note:  The  stem  dust  will  penetrate  chi  .min  cloud  base 
crotn  (a)  enter  cable  and  find  l- 


9 

2-73 

♦ 

0 . 39x  /  u  •  0.5 

a,  itc 

Cat 

- 

A 

9 

2.73 

0 . 23x  \  «  1.0 

k/scc 

Cat 

- 

A 

■9 

2.73 

•v 

o.u*  {  -i.3 

»/  see 

Cat 

- 

A 

2.72 

•F 

0.13*  j  »  2.0 

3/S3C 

Cat 

- 

A 

« 

2.73 

* 

0.137v  I  2.0 

5/  see 

Cat 

- 

A 

» 

2.73 

♦ 

0.U*  ' 

Cat 

w 

2 

- 

2.73 

0.073s 

cm 

- 

i* 

«• 

Hm 

2.73 

♦ 

0.06s 

cm 

- 

5 

9 

2.73 

O.OpJs 

■Cm 

- 

V 

w 

2.73 

♦ 

O.Oiex 

Cat 

- 

V 
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(e)  Cloud  Volume:  Stem 


V ( t)  =  (tt/6)  x  ♦  •  z1 

x  *  uc  +  D  /2 
s 

y"  *  y(t)  +  Ds 

z1  -  z(t)  +  Hs 

(f)  Concentration  (t)  -  0.1  - yjt) — C —  S'm 

(8)  Main  Cloud  Dispersion 

(a)  Compute  u(x) 

u(z)  =  u  (4m)  •  (z/4)^ 

p  =  0.1,  i.  B,  and  C 

p  *  0.142,  Cat  D 
p  =  0.3,  Cat  E 
p  =  0.5,  Cat  F 

p  s  0.8,  Cat  G  (+3,  Table  1) 

(b)  For  unstable  and  neutral,  cloud  sheared  at  base  (4m), 
H  (top)  and  (H^  -  4)  f  2. 

For  stable  categories,  shear  only  top. 

From  input  data,  find  Pasquill  Category,  see  (7a). 
x  =  ut  where  u  is  given  by  above  rules 

y,  see  (7c),  again  use  appropriate  value  of  x. 

z,  see  (?d),  again  use  appropriate  value  of  x. 

(c)  Inversion  layers  can  be  treated  as  defined  for  smoke. 

(9)  Main  Cloud  Dimensions 

x  ■  ut  +  D  , 1 2 
cl 

y1  *  y (t)  +  Dcl 

2  (t)  -  Ah(t)  +  5(c) 
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(10) 


Main  Cloud  Volume 


A  simple  approach  is  to  assume  an  ellipsoidal  shape;  i.e., 
ellipsoid  r  4,  where  x,  y,  and  Z(t)  are  defined  for  the  top 
layer.  This  will  give  a  somewhat  larger  volume.  '  Otherwise, 
we  can  treat  the  volume  in  two  parts,  the  top  layer  and  the 
remainder. 

For  top,  insert,  compute  values  for  x,  y,  Z 
Vx  =  (tt / 6)  xy 1  Z(t) 

For  remainder 

TT  D  ?T  . 

V2  =  +  y  (t)  *  z(c)  *  u  *  t 

where  u  is  at  cloud  center. 

(11)  Main  Cloud  Concentration 

Concentration  ~  0.9  loft  weight  *  (V^  +  V^) 

Note:  The  validations  in  the  main  text  were  performed  using 
the  simple  procedure. 

(12)  Particle  Size  Distribution 

Follow  procedure  described  in  main  text. 

(13)  For  a  precise  calculation  where  the  charge  center  height  is 
known  explicitly,  then  the  crater  dimensions  and  crater 
efficiency  can  be  estimated  from  material  in  text  and  appendix  D. 
The  remainder  of  the  procedure  is  as  given  above. 

The  values  in  Fig.  4,  Part  2,  are  for  free  air  explosions. 

For  the  condition  HOB  <  0  the  reader  must  employ  the  relation 

Wair  =  W(charge)  •  exp  (-3pAd) 

1/3 

where  p  =  soil  specific  gravity  and  »  charge  depth/W 
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(14)  Multi-Munitions 

In  the  above  procedure,  we  assumed  each  explosion  is  mutually 
independent.  This  is  correct  provided  the  charges  are  about 
50  m  apart,  or  the  material  velocity  (see  Fig.  4)  is 
comparable  to  the  ambient  u.  Otherwise,  there  are  shock 
interactions.  Since  operationally  there  is  no  a  priori 
knowledge  of  detonation  locations  and  lapse  time  between  these 
events;  i.e.,  space  and  time  must  be  considered,  the  effects 
cannot  be  modeled  until  a  detailed  study  of  the  statistical 
characteristics  are  made. 

Again  the  basic  data  to  analyze  the  interactive  event 
is  contained  in  the  text,  and  the  reader  may  analyze  any 
particular  event. 


/ 
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PART  3 

VEHICULAR  DUST  OBSCURANT  MODEL* 


i 


Note:  The  basic  work  was  performed  under  sponsorship  of  the  Defense 
Advanced  Research  Projects  Agency  and  Is  Incorporated  herein 
for  completeness. 
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INTRODUCTION 


Experience  has  shown  that  military  tracked  vehicles  and  wheeled 
vehicles  raise  considerable  dust  as  evidenced  at  F roving  Grounds  and  in 
combat  areas.  Under  some  conditions,  the  generated  obscuration  creates  a 
serious  problem  for  offensive  systems.  We  will  consider  this  problem  of 
vehicle  dust  obscuration. 

Ground  vehicles  can  raise  loose  surface  material  (dust,  water,  etc.). 
The  atmospheric  loading  of  this  materi',1  is  dependent  on  the  type  of  vehicle, 
vehicle  speeds,  fender  arrangement,  etc.  For  dust,  we  must  add  the  soil 
type,  soil  cohesiveness,  grain  size,  and  related  factors.  The  launching 
speed  and  angle  in  combination  with  atmospheric  conditions,  wind  speed,  and 
stability,  will  affect  the  dust  cloud  behavior.  For  example,  particles  of 
~  7  pm  launched  at  a  minimum  speed  of  ~7  m/sec  will  be  transported  with  air- 

-ft 

flow  motion;  larger  particles  are  propagated  by  a  process  called  saltation. 
The  latter  is  a  looping  action  where  the  particle  impacts  the  soil  and 
launches  additional  soil  particles. 

The  quantity  and  quality  of  data  on  dust  created  by  vehicles  are 
extremely  limited;  that  is,  there  are  no  real  controlled  experiments  nor 
sufficient  statistical  data  to  generate  a  comprehensive  model.  Further,  the 
soil  types  are  limited,  namely,  to  those  of  Army  Proving  Grounds.  The 
general  configuration  of  the  vehicle  will  affect  the  dust  cloud  geometry. 

For  the  available  data,  this  parameter  has  not  been  specifically  defined 
except  in  a  general  sense,  i.e.,  light  or  heavy  wheeled,  heavy  track.  Never¬ 
theless,  an  attempt  has  been  made  here  to  provide  several  algorithms  to 
permit  an  estimate  of  vehicle  dust. 


See  Bagnold's  Blown  Sand  &  Desert  Dunes  (1934). 


2  TECHNICAL  DISCUSSION 

First,  to  generate  a  dust  cloud  relatively  loose,  dry  soil  is  required - 
Second,  the  soil  particle  size  distribution  affects  not  only  the  attenuation 
coefficient  but  also  the  dust  concentration  lofted  into  air  and  staying 
aloft.  In  Fig.  1  is  presented  the  relation  between  concentration  and  p, 
percent  by  weight  of  particles  with  d  <  74  ym,  where  the  areas  define  the 
cloud  visibility  at  —25  ft  from  the  source.  We  see  if  p  <  10%,  there  is 
unrestricted  visibility  except  for  high  concentration.  On  the  other  hand, 
if  p  >,  24%,  the  visibility  reduces  to  s  0  ft.  This  is  understandable  since 
the  more  loose  powder  in  the  surface  and  subsurface  layers,  the  denser  the 
dust  cloud.  Thus  we  suggest  Fig.  1  is  a  general  plot  which  may,  however, 
differ  somewhat  for  a  specific  situation. 

In  Part  2,  Table  6,  we  presented  the  particle  size  fractions  at 
several  US  Army  facilities.  The  fraction  of  particles  £74  um  falls  into 
the  category  of  dust  cloud  potential.  Similarly,  we  have  shown  that  at 
Grafenwohr  the  soil  has  been  sufficiently  pulverized  at  a  number  of  sites 
to  allow  dust  cloud  generation.  Of  course,  the  soil  should  be  relatively 
dry  to  be  conducive  for  dust  making.  It  should  be  noted  that  clay  soils 
tend  to  have  large  particle  sizes  (HMD  —  650  um)  and  high  cohesiveness 
when  wet.  Therefore,  the  likelihood  for  dust  generation  differs  radically 
depending  on  locale. 

From  Fig  1  we  can  estimate  that  the  photopic  extinction  coefficient 

2 

(in  m  /g)  will  vary  significantly;  see  Table  1. 


TABLE  1 

EXTINCTION  COEFFICIENT  AS  A  FUNCTION  p 


1  <  74  yo) 

»e(®  /g) 

>21 

>0.93 

9-21 

<0.093-0.93 

<9 

>0.093 
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Figure  i.  Visibility  Zones  as  a  Function  or  Concentration  and  p 
(Estimates  or  Extinction  Coefficients  Are  Provided) 


Ic  should  b-i  noted  chat  for  typical  soil  p  —  15%  and,  therefore,  a  cr  .46  m"/g 
is  a  valid  number.  The  spread  in  Table  1  is  very  interesting;  it  shows  an 
order-of-magnitude  difference  in  extinction  coefficient;  see  Part  2. 

Next,  we  analyzed  the  concentration  as  a  funccion  of  distance  for 
heavy  wheeled  and  heavy-tracked  vehicles,  where  the  vehicle  speeds  were 
between  15  and  35  mph  and  heavy  dust  conditions  existed.  The  expressions 
are  given  below  where  the  concentration  is  constant  from  vehicle  axis  to 
edge  (halfwidth) . 

Wheeled  Vehicle: 

C  (g/m3)  -  6. 9489/d3’ 32  (1) 

Heavy  Track: 

C  (g/m3)  -  525.3/d2’367  (2) 

where  d  is  the  distance  in  meters  from  the  vehicle  velocity  vector  and  is  >. 
vehicle  width/2.  When  the  surface  dust  is  light,  the  values  are  reduced  by 
an  order  of  magnitude. 

Lastly,  we  have  reduced  a  set  of  data  for  heavy-wheeled  vehicles  in 
light  soil  that  related  the  number  of  vehicles  in  a  convoy,  vehi>.  '  e  speed, 
duration  of  dust,  and  distance  from  vehicle  route,  i.e.,  N»  u  (mph,  >  c  (sec), 
and  d  (ft),  respectively. 

The  general  form  of  the  expression  is 

N/t  -  (N/t)asy  (1  -  e’3u/dj  (3/ 

where  3  is  an  empirical  determined  constant.  The  preliminary  values  of  (N/t) 
asymptotic,  i.e.,  limiting  value,  and  3  are 
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Vehicle  Type 

Soil  Type 

(N/C) 

asy 

3 

Heavy  Wheel 

Light/Moderate 

0.2 

1.8 

Light  Track 

Moderate 

0.24 

2.86 

Light  Wheel 

Light 

0.14 

2.0 

Heavy  Track 

Heavy 

0.28 

5.0 

The  soil  cypes  are  defined  as  follows: 

Light  p  <  52 

Moderate  p  <  92 

Heavy  p  >  102 

One  should  note  that  c  and  d  are  positive  (non-sero) ,  Lastly,  the 

effect  of  wind  is  minor  and  included  in  the  above,  i.e.,  for  winds  <  5  m/sec. 

The  several  expressions  given  above  provides  an  initial  approach  to 
the  modeling  of  vehicle  dust.  Basically,  the  visual  observation  produced 
by  generated  dust  has  a  duration  comparable  to  the  time  tne  vehicle  passes 
through  the  line  of  sight.  Secondly,  the  extrapolation  to  another  wavelength, 
\,  the  required  concentration,  would  have  to  be  increased  (decreased) 
directly  proportional  to  the  ratio  of  the  extinction  coefficients.  Finally, 
the  above  algorithms  are  valid  to  cloud  heights  of  3  m. 
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MODEL  STRUCTURE 

To  apply  the  above  expressions,  we  suggesc  Che  following  procedure. 
Inputs : 

Vehicle  Speed  -  u  (mph) 

Vehicle  Type  -  wh«#l.<*d  or  cracked 

Vehicle  Size  -  light,  heavy 

Soil  Condition  -  light,  heavy 

Soil  Type  -  p,  %  wt  of  particles  (d  £  74  ym) 

Number  of  Vehicles  -  N 

Steps : 

1.  If  the  soil  is  either  wet  or  u  <  10  mph  or  p  <  902,  then 
the  concentration  »  0. 

2.  Establish  value  of  p,  soil  condition,  location  from  source 
edge  d,  type  of  vehicle. 

3 

Enter  Eqs.  1  and  2  and  find  concentration.  (Note:  g/m  » 

35  g/ft3.) 

-5  3 

If  concentration  <  10  g/ft  and  d  >  25  ft,  terminate. 

-5  3 

If  concentration  >  10  g/ft  ,  find  C  vs.  d.  C  •  L  is 
determined  by  geometry. 

3.  To  find  transmission: 

Multiply  C  •  L  by  appropriate  extinction  coefficient.  If 
the  region  of  interest  is  the  photopic,  then  enter  Fig.  I 
and  find  visibility  range. 

4.  To  find  the  duration: 

Enter  Eq.  3  with  N,  u,  and  d  and  (N/t)  and  §  and 

asy 

compute  t. 
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MUNITIONS  CHARACTERISTICS 


Table  A2 

Some  Smoke  Generators  and  Pots 


1 


Smoke 

i 

1  Type 

i  Fin 

'  Fill  Wt/Flow  Pate 

1  Smoxe  Generator 

M3 

|  Fog 

Oil 

1 

25-50  gal /hr 

Smoke  Pot  ; 

AN-M7 

Fog 

Oil 

(SGF)  ; 

— 

1  j 

AN-M7A1 

i  Fog 

Oil 

(SGF)  ; 

— 

Smoke  Pot  I 

Ml 

j  HC 

1 

1 

i 

10.0 

! 

ABC -MS 

j  HC 

i 

i 

■ 

3.0 

!  Air  Smoke  Pot 

M4A2 

! 

1 

1 

Table  A3 

UK  Grenades 

Smoke 

i 

Burn 

! 

i 

! 

Screen 

Range  ; 

Time 

1 

1 

Time 

Type . . . 

(m) 

(sec) 

Source  ; 

Smoke  1 

(sec) 

:  Fi  1 1 

L5A1 ,2,3,A 

70  (grd)  ! 

55-75 

Point  1 

Whi  ta  j 

3 

i 

i 

i 

ground  i 

1 

! 

! 

l7A1 

70  (grd)  ! 

30-50 

Point  | 

Green  | 

i 

i 

t 

| 

1 

I 

1 

ground  | 

j 

.<17-2 

120  (grd)  ! 

30-50 

Point  I 

1 

Green  j 

i 

* 

i 

ground  i 

| 

XL361 

30  i 

60-90 

3  m  air,  ! 

•  ; 

2 

j  RP» 

'Certain ) 

1 

ground 

j 

,  butyl  rucoer 

|  L3A1 

| 

i 

. 

i 

i 

XL13EI 

i 

120  ; 

1 0  m  a  i  r ,  ; 

| 

!  nCV 

,  [for  L3A1 

| 

ground  i 

1 

♦ 

t  m 

j  silica 

j  practice) 

j 

l 

} 

XL  136 : 

no  ; 

3-10 

1 

i 

! 

Green  ! 

i 

t 

i 
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APPENDIX  B 

A  STABILITY  CLASSIFICATION  BASED  ON 
HOURLY  AIRPORT  OBSERVATIONS 


m 


A  STABILITY  CIaSSIJICATICM  3ASZ3  Q'J 
HCU3LY  AlBBCBT  CSoir/AlICh'S  a 


This  system  of  classifying  scab iiicy  on  an  hourly  basis  for  research 
in  air  pollution  is  based  upon  work  accomplished  by  3r.  ?. 

Pasquili  of  the  British  Meteorological  Office  (see  reference  1  of 
text) .  Stability  near  the  ground  is  dependent  primarily  upon 
net  radiation  and  wind  speed.  Without  the  Influence  of  clouds, 
insolation  (incoming  radiation)  during  the  day  is  dependant 
upon  solar  altitude,  which  is  a  function  of  time  of  day  and 
time  or  year.  When  clouds  exist,  their  cover  and  thickness 
decrease  incoming  and  outgoing  radiation.  In  this  system 
insolation  is  estimated  by  solar  altitude  and  modified  for 
existing  conditions  of  total  cloud  cover  and  cloud  ceiling 
height.  At  night,  estimates  of  outgoing  radiation  are  made 
by  considering  cloud  cover.  This  stability  classification 
system  has  been  made  completely  objective  so  that  an  electronic 
computer  can  be  used  to  compute  stability  classes.  The  stability 
classes  are  as  follows:  (A)  Bxtremely  Unstable,  (3)  Unstable, 

(0  Slightly  Unstable,  (D)  Neutral,  (a)  Slightly  Stable, 

(?)  Stable,  (G)  Bxtreaeiy  Stable.  Table  8-1  gives  the 
stability  class  as  a  function  of  wind  speed  a ns  net  radiation. 

The  net  ra’istion  index  ranges  from  4,  highest  positive  net 
radiation  (directed  coward  th,e  ground),  to  -2,  highest  negative 
net  radiation  (directed  away  from  the  earth) .  Instability 
occurs  with  high  positive  net  radiation  and  law  wind  speed, 
stability  with  high  negative  net  radiation  and  light  winds, 
and  neutral  conditions  with  cloudy  skies  or  high  wind  speeds. 

The  net  radiation  index  used  with  wind  speed  to  ooca.'n  stability 
class  i3  determined  by  the  following  procedure: 


*  The  following  explanation  of  the  Pas^oill  Stability  c-assirication 
has  been  extracted  from  an  article  by  3.  Bruce  Turner  in  the 
February  1H4  .’oumal  of  Art  lied  Meteorology 


1 1 


(i) 


chan  7CC0 
or  night) 


Xi  she  coca!  cloud  cover  is  10/10  and  che  calling 
£aac,  use  nec  radiacian  index  equal  co  0  (whether 


is  less 
day 


(2)  For  night rise  (night  is  dafinad  as  che  period  iron  one 
hour  before  sunset  co  one  hour  after  sunrise) : 

a.  If  coca!  cloud  cover  <4/10,  use  nec  radiation  index 

equal  to  -2. 

b.  If  total  cloud  cover  >4/10,  use  nec  radiation  index 

equal  to  -1. 

(3)  For  daytime : 

a.  Determine  che  insolation  class  number  as  a  function 
of  solar  altitude  from  Table  3-2. 

b.  If  total  cloud  cover  <5/10,  use  che  nec  radiation  index 
in  Table  A-I  corresponding  to  che  insolation  class  number. 

c.  If  cloud  cover  >5/10,  modify  the  insolation  class  numb  a 
by  following  these  six  steps: 

(1)  Ceiling  < TOGO  ft.,  subtract  2. 

(2)  Ceiling  >7000  ft.  but  <  loGCO  ft.,  subtract  1. 

(3)  Total  cloud  cover  equal  10/10,  subtract  1.  (This 
vill  only  apply  to  ceilings  >7000  ft.  since  cases 

10/10  coverage  below  7000  ft.  are  considered 
in  ices  1  above) . 

(4)  If  insolation  class  number  has  not  been  modified 
by  steps  (1),  (2).  or  (35  above,  assume  codified 
class  number  equal  to  insolation  class  number. 

(5)  If  modified  insolation  class  number  is  less  than 
1,  let  i;  equal  1. 

Co)  Use  the  net  radiation  index  in  Table  3-1 

corresponding  to  the  modified  insolation  cia*3 
number. 


103 


v::ro  spzid 
(SHOTS) 

4 

3 

RAD  14: 

2  T 

o,  : 

A 

A 

3 

C 

2,  3 

A 

3 

3 

c 

■*»  5 

A 

3 

n 

•j 

6 

3 

3 

c 

3 

7 

3 

3 

c 

3 

8,  9 

3 

C 

c 

D 

10 

C 

C 

3 

3 

*  1 

C 

C 

u 

* 

> 

c 

3 

3 

* 

TA3LS  3-3 

insola: 

:o.v 

AS  A 

?: 

SCLA&  XLZZ72Z 


60  *  «■  a 
35*  <  a 
15*  <  a 


a 


Scrar-s 

60  <  y.Q‘i&?zza 
35  i  Sligh: 

V«ait 


<  * 2* 


t/j  o 


APPENDIX  C 


CONTINUOUS  VALUES  FOR  PLUME  DISPERSION  PARAMETERS  PLUME  HEIGHT 


-a  Table  10  of  Parc  1  we  have  defined  che  vertical  rise  as  1.73  - 
A  •  x,  where  A  depends  upon  Pasquill  category.  For  modeling  purposes  we 
have  proposed  a  fixed  value  for  each  category.  Actually  the  parameter 
varies  continuously  with  mean  wind  speed,  u,  for  values  <  5  3/ sec. 
Although  smoke  normally  is  not  discharged  when  u  1  1  m/sec  because  the 
horizontal  coverage  is  low,  f-'r  completeness  we  present  the  value  of  A 
for  0.5  5.  u  -  5  a/ sec.  (Fig.  C.l).  Similarly,  we  provided  an  estimace  of 
the  continuous  values  for  the  lateral  dispersion,  y  “0.1+3  •  x. 


<  0.3 


UNSTABLE 
o  NEUTRAL 


£ 

3  0.2 

< 


STABLE 


-Vi 


u  (m/ $dC) 


Figure  C.l.  Pluae  Parameters  at  Low  Wind  Speeds;  Rise  < A; 
Lateral  Spread  1 3)  “*  •“ 
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ADDITIONAL  CRATERING  DA 


The  US  Army,  Waterways  Experimental  Station,  analyzed  the  available 

data  through  1960  and  reported  the  results  in  Ref.  1.  The  data  included 

results  from  various  explosives  including  shaped  charges,  various  3oils 

and  terrain  conditions,  and  charge  weights  from  <  1  lb  to  10  lb.  The 

results,  in  algebraic  form,  are  given  in  Table  D.l,  where  \  *  scaled 

1/3  C 

charge  center  location  in  units  of  ft/charge  weight  .  It  should  be 
noted  that  the  3a  values  for  the  crater  radii  and  depths  vary  between 
1.5*  to  4x  the  mean  value. 

In  this  appendix  several  nomographs.  Figs.  D.l  and  D.2,  are  reproduced 
here  from  Ref.  1  that  show  the  differences  between  wet  and  dry  soil  types. 

Additional  data  are  provided  in  Tables  D.2  and  D.3  and  Figs.  D.3  and 

D.4. 


‘US  Army,  Waterways  Experiment  Station  (WES) ,  Analysis  of  Crater  Data, 
Technical  Report  2-347,  Report  2,  June  1961. 
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CRATER  DATA 
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NOMOGRAPH  FOR  CRATER 
DEPTH  AND  RADIUS 

apparent  crater  in  dry  clay 


Figure  D.l.  Nomograph  for  Crater  Depth  and  Radius  (US  Amy,  Waterways 

Experiment  Station  (WES) ,  Analysis  of  Crater  Data,  Technical 
Report,  2-547,  Report  2,  June  1961. 
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